楊長(zhǎng)森

（青海新高度房地產(chǎn)投資有限公司，青海 西寧 810000）

在高層建筑中出現(xiàn)大體積混凝土結(jié)構(gòu)的部位主要是基礎(chǔ)底板。這種底板的面積較大，厚度也在80 cm以上，最厚處可達(dá)5～6 m?；A(chǔ)底板澆筑后，在內(nèi)部水化熱和外部環(huán)境溫度變化的影響下，混凝土結(jié)構(gòu)各處溫度時(shí)時(shí)發(fā)生變化，其內(nèi)部最高溫度將達(dá)到50～80 ℃，與外界氣溫相差很大。如果不采取一定的措施，混凝土很容易開(kāi)裂，影響結(jié)構(gòu)的整體性并造成地下水滲漏，進(jìn)一步影響結(jié)構(gòu)的耐久性，給整個(gè)工程帶來(lái)極其嚴(yán)重的危害。

混凝土熱分析不僅是預(yù)測(cè)混凝土溫度變化的必要方法，同時(shí)也是計(jì)算混凝土溫度應(yīng)力的前提條件。由于早期混凝土各項(xiàng)熱性能都時(shí)時(shí)發(fā)生變化，同時(shí)環(huán)境等因素也時(shí)時(shí)在變化，因此，混凝土早期熱分析是一個(gè)較為復(fù)雜的問(wèn)題。ANSYS軟件雖具有專(zhuān)門(mén)的熱分析模塊，適用于計(jì)算一個(gè)系統(tǒng)或部件的溫度分布及其他熱物理參數(shù)，如熱量獲取或損失、溫度梯度、熱流密度等，但當(dāng)運(yùn)用在混凝土早期熱分析這一特定問(wèn)題時(shí)，尚需開(kāi)展一定的研究。［1］

1 軟件二次開(kāi)發(fā)

ANSYS熱分析基于能量守恒原理和熱平衡方程，通過(guò)施加溫度、熱流率、熱流密度、對(duì)流、輻射、絕熱和生熱等邊界條件和初始條件，求解即可得到節(jié)點(diǎn)的溫度、溫度梯度、熱流率等結(jié)果。

ANSYS熱分析過(guò)程可分為3個(gè)步驟，前處理（建模）、加載求解（施加荷載并計(jì)算）和后處理（查看結(jié)果）。在建模階段，主要工作為定義單元類(lèi)型并設(shè)置關(guān)鍵選項(xiàng)，定義單元實(shí)常數(shù)，定義材料熱性能，建立幾何模型，對(duì)幾何模型劃分網(wǎng)格。在加載求解階段，主要工作為定義分析類(lèi)型，獲得瞬態(tài)熱。

分析的初始條件，設(shè)定荷載步選項(xiàng)和存盤(pán)求解。在后處理階段，主要工作為查看計(jì)算結(jié)果，ANSYS提供了通用后處理（POST1）和時(shí)間后處理（POST26）兩種。在POST1模塊中，可以對(duì)整個(gè)模型在某一荷載步的結(jié)果進(jìn)行后處理，對(duì)結(jié)果進(jìn)行彩色云圖顯示、矢量圖顯示、打印列表等，如果設(shè)定的時(shí)間點(diǎn)不在任何一個(gè)子步的時(shí)間點(diǎn)上，ANSYS會(huì)進(jìn)行線(xiàn)性插值。在POST26模塊中，可以繪制計(jì)算結(jié)果隨時(shí)間變化的曲線(xiàn)。

混凝土早期熱分析施加荷載常用命令見(jiàn)表 1，其中施加生熱量較為復(fù)雜。早期混凝土由于水化反應(yīng)產(chǎn)生大量水化熱，在ANSYS中模擬混凝土早期溫度場(chǎng)時(shí)，水化熱量作為荷載輸入，由于每個(gè)時(shí)刻產(chǎn)生的熱量不同，因此首先要將不同時(shí)刻的生熱量做成數(shù)組，在bf命令中使用%tabname%（tabname為數(shù)組名）施加隨時(shí)間變化的生熱量。另外，還可以采用 UPFs進(jìn)行二次開(kāi)發(fā)，修改Usrefl.F程序（主要參數(shù)見(jiàn)表2）。Usrefl.F程序可用于標(biāo)量場(chǎng)（溫度、生熱量、濕度等）荷載的修改和施加，該程序中的key參數(shù)用于確定具體修改哪個(gè)物理量，當(dāng)key＝3時(shí)，輸入的物理量為放熱量，將其寫(xiě)成荷載步的函數(shù)，在每一個(gè)迭代步中都會(huì)調(diào)用該子程序，從而施加水化熱量。［2］

表1 熱分析施加荷載常用命令

表2 USREFL.F參數(shù)說(shuō)明

為了驗(yàn)證 ANSYS有限元計(jì)算結(jié)果的正確性，選取了有理論解的算例與有限元計(jì)算結(jié)果進(jìn)行對(duì)比。［3］

算例：處于絕熱溫升狀態(tài)立方體混凝土塊尺寸為10 m×10 m×10 m，初始溫度為150 ℃，混凝土的絕熱溫升方程為θ＝30（1－e－0.34τ），其中齡期τ單位為 d，導(dǎo)溫系數(shù)α＝0.004 m2/h，導(dǎo)熱系數(shù)λ＝10.76 kJ/（m·h·℃），比熱 c＝1.00 kJ/（kg·℃）。ANSYS求解所得結(jié)果見(jiàn)圖1，兩種計(jì)算結(jié)果比較見(jiàn)表3。從計(jì)算結(jié)果看，ANSYS的計(jì)算結(jié)果與理論值基本相同。

表3 理論解與ANSYS解比較

2 實(shí)例分析

兩幢12層和26層高的辦公樓，設(shè)有一個(gè)單層地下停車(chē)場(chǎng)，總建筑面積約為46 056 m2，地下室底板平均厚度2.5 m，電梯井處最厚達(dá) 3.6 m，屬典型大體積混凝土結(jié)構(gòu)；地下室側(cè)墻厚400 mm，高6 m，周長(zhǎng)400多m，屬典型超長(zhǎng)薄壁結(jié)構(gòu)。根據(jù)大體積混凝土澆筑方案和某商品混凝土公司提供的配合比情況，按施工計(jì)劃地下室基底板共設(shè)三條后澆帶，A樓深基礎(chǔ)及B1－B4，BE－BF部位大底板，自然分成5個(gè)區(qū)域，根據(jù)區(qū)域順序先后及基底施工順序，先深后淺、逐漸按劃分區(qū)域展開(kāi)分別澆筑鹼。其中5號(hào)區(qū)域底板最厚，達(dá)3.4 m。第1區(qū)域混凝土分兩次澆筑，第一次澆1.6 m，第二次澆1.8 m，兩次間隔時(shí)間約4 d。采用ANSYS進(jìn)行有限元分析。首先建立有限元模型見(jiàn)圖2，有限元分析結(jié)果與實(shí)測(cè)結(jié)果比較見(jiàn)圖3～圖5，從圖分析，除混凝土表面溫度用環(huán)境影響變化比較突然外，其他各點(diǎn)有限元分析結(jié)果與測(cè)試結(jié)構(gòu)基本一致。

圖2 底板有限元模型

3 結(jié)束語(yǔ)

ANSYS軟件是國(guó)際著名的通用有限元軟件，是第一個(gè)通過(guò)ISO 9001質(zhì)量認(rèn)證的大型分析設(shè)計(jì)類(lèi)軟件，是美國(guó)機(jī)械工程師協(xié)會(huì)（ASME）、美國(guó)核安全局（NQA）及近二十種專(zhuān)業(yè)技術(shù)協(xié)會(huì)認(rèn)證的標(biāo)準(zhǔn)分析軟件。ANSYS雖然通用性很強(qiáng)，但是在計(jì)算早期混凝土結(jié)構(gòu)溫度場(chǎng)、溫度應(yīng)力場(chǎng)這一特定問(wèn)題時(shí)略顯不足，［4］不能直接求解，需要人為地編程干預(yù)并進(jìn)行開(kāi)發(fā)?；贏NSYS的二次開(kāi)發(fā)功能，編寫(xiě)適用于混凝土早期溫度計(jì)算的模塊，充分利用 ANSYS的前后處理功能，減少軟件開(kāi)發(fā)的工程量，模擬混凝土早期溫度變化過(guò)程，對(duì)大體積混凝土溫度控制具有重要的意義。

圖3 1－1號(hào)測(cè)溫點(diǎn)溫度實(shí)測(cè)與模擬計(jì)算

圖4 1－2號(hào)測(cè)溫點(diǎn)溫度實(shí)測(cè)與模擬計(jì)算

圖5 1－3號(hào)測(cè)溫點(diǎn)溫度實(shí)測(cè)與模擬計(jì)算

1 祝效華.ANSYS高級(jí)工程有限元分析范例精選［M］.北京：電子工業(yè)出版社，2004：26～43

2 李驍春.基于數(shù)值分析和工程經(jīng)驗(yàn)的現(xiàn)代大體積砼裂縫控制實(shí)用技術(shù)［J］.河海大學(xué)，2008

3 李曉春、昊勝興.基于水化度概念的早期混凝土熱分析［J］.科學(xué)技術(shù)與工程，2008（2）：441～445

4 李袋春、吳勝興.基于ALYSYS的混凝土早期徐變應(yīng)力仿真分析［J］.系統(tǒng)仿真學(xué)報(bào)，2008（15）：3944～3947