楊長森
(青海新高度房地產(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 ℃,與外界氣溫相差很大。如果不采取一定的措施,混凝土很容易開裂,影響結(jié)構(gòu)的整體性并造成地下水滲漏,進(jìn)一步影響結(jié)構(gòu)的耐久性,給整個(gè)工程帶來極其嚴(yán)重的危害。
混凝土熱分析不僅是預(yù)測混凝土溫度變化的必要方法,同時(shí)也是計(jì)算混凝土溫度應(yīng)力的前提條件。由于早期混凝土各項(xiàng)熱性能都時(shí)時(shí)發(fā)生變化,同時(shí)環(huán)境等因素也時(shí)時(shí)在變化,因此,混凝土早期熱分析是一個(gè)較為復(fù)雜的問題。ANSYS軟件雖具有專門的熱分析模塊,適用于計(jì)算一個(gè)系統(tǒng)或部件的溫度分布及其他熱物理參數(shù),如熱量獲取或損失、溫度梯度、熱流密度等,但當(dāng)運(yùn)用在混凝土早期熱分析這一特定問題時(shí),尚需開展一定的研究。[1]
ANSYS熱分析基于能量守恒原理和熱平衡方程,通過施加溫度、熱流率、熱流密度、對流、輻射、絕熱和生熱等邊界條件和初始條件,求解即可得到節(jié)點(diǎn)的溫度、溫度梯度、熱流率等結(jié)果。
ANSYS熱分析過程可分為3個(gè)步驟,前處理(建模)、加載求解(施加荷載并計(jì)算)和后處理(查看結(jié)果)。在建模階段,主要工作為定義單元類型并設(shè)置關(guān)鍵選項(xiàng),定義單元實(shí)常數(shù),定義材料熱性能,建立幾何模型,對幾何模型劃分網(wǎng)格。在加載求解階段,主要工作為定義分析類型,獲得瞬態(tài)熱。
分析的初始條件,設(shè)定荷載步選項(xiàng)和存盤求解。在后處理階段,主要工作為查看計(jì)算結(jié)果,ANSYS提供了通用后處理(POST1)和時(shí)間后處理(POST26)兩種。在POST1模塊中,可以對整個(gè)模型在某一荷載步的結(jié)果進(jìn)行后處理,對結(jié)果進(jìn)行彩色云圖顯示、矢量圖顯示、打印列表等,如果設(shè)定的時(shí)間點(diǎn)不在任何一個(gè)子步的時(shí)間點(diǎn)上,ANSYS會進(jìn)行線性插值。在POST26模塊中,可以繪制計(jì)算結(jié)果隨時(shí)間變化的曲線。
混凝土早期熱分析施加荷載常用命令見表 1,其中施加生熱量較為復(fù)雜。早期混凝土由于水化反應(yīng)產(chǎn)生大量水化熱,在ANSYS中模擬混凝土早期溫度場時(shí),水化熱量作為荷載輸入,由于每個(gè)時(shí)刻產(chǎn)生的熱量不同,因此首先要將不同時(shí)刻的生熱量做成數(shù)組,在bf命令中使用%tabname%(tabname為數(shù)組名)施加隨時(shí)間變化的生熱量。另外,還可以采用 UPFs進(jìn)行二次開發(fā),修改Usrefl.F程序(主要參數(shù)見表2)。Usrefl.F程序可用于標(biāo)量場(溫度、生熱量、濕度等)荷載的修改和施加,該程序中的key參數(shù)用于確定具體修改哪個(gè)物理量,當(dāng)key=3時(shí),輸入的物理量為放熱量,將其寫成荷載步的函數(shù),在每一個(gè)迭代步中都會調(diào)用該子程序,從而施加水化熱量。[2]
表1 熱分析施加荷載常用命令
表2 USREFL.F參數(shù)說明
為了驗(yàn)證 ANSYS有限元計(jì)算結(jié)果的正確性,選取了有理論解的算例與有限元計(jì)算結(jié)果進(jìn)行對比。[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é)果見圖1,兩種計(jì)算結(jié)果比較見表3。從計(jì)算結(jié)果看,ANSYS的計(jì)算結(jié)果與理論值基本相同。
表3 理論解與ANSYS解比較
兩幢12層和26層高的辦公樓,設(shè)有一個(gè)單層地下停車場,總建筑面積約為46 056 m2,地下室底板平均厚度2.5 m,電梯井處最厚達(dá) 3.6 m,屬典型大體積混凝土結(jié)構(gòu);地下室側(cè)墻厚400 mm,高6 m,周長400多m,屬典型超長薄壁結(jié)構(gòu)。根據(jù)大體積混凝土澆筑方案和某商品混凝土公司提供的配合比情況,按施工計(jì)劃地下室基底板共設(shè)三條后澆帶,A樓深基礎(chǔ)及B1-B4,BE-BF部位大底板,自然分成5個(gè)區(qū)域,根據(jù)區(qū)域順序先后及基底施工順序,先深后淺、逐漸按劃分區(qū)域展開分別澆筑鹼。其中5號區(qū)域底板最厚,達(dá)3.4 m。第1區(qū)域混凝土分兩次澆筑,第一次澆1.6 m,第二次澆1.8 m,兩次間隔時(shí)間約4 d。采用ANSYS進(jìn)行有限元分析。首先建立有限元模型見圖2,有限元分析結(jié)果與實(shí)測結(jié)果比較見圖3~圖5,從圖分析,除混凝土表面溫度用環(huán)境影響變化比較突然外,其他各點(diǎn)有限元分析結(jié)果與測試結(jié)構(gòu)基本一致。
圖2 底板有限元模型
ANSYS軟件是國際著名的通用有限元軟件,是第一個(gè)通過ISO 9001質(zhì)量認(rèn)證的大型分析設(shè)計(jì)類軟件,是美國機(jī)械工程師協(xié)會(ASME)、美國核安全局(NQA)及近二十種專業(yè)技術(shù)協(xié)會認(rèn)證的標(biāo)準(zhǔn)分析軟件。ANSYS雖然通用性很強(qiáng),但是在計(jì)算早期混凝土結(jié)構(gòu)溫度場、溫度應(yīng)力場這一特定問題時(shí)略顯不足,[4]不能直接求解,需要人為地編程干預(yù)并進(jìn)行開發(fā)?;贏NSYS的二次開發(fā)功能,編寫適用于混凝土早期溫度計(jì)算的模塊,充分利用 ANSYS的前后處理功能,減少軟件開發(fā)的工程量,模擬混凝土早期溫度變化過程,對大體積混凝土溫度控制具有重要的意義。
圖3 1-1號測溫點(diǎn)溫度實(shí)測與模擬計(jì)算
圖4 1-2號測溫點(diǎn)溫度實(shí)測與模擬計(jì)算
圖5 1-3號測溫點(diǎn)溫度實(shí)測與模擬計(jì)算
1 祝效華.ANSYS高級工程有限元分析范例精選[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