高軒能,黃文歡,張惠華
(華僑大學(xué)土木工程學(xué)院,福建泉州 362021)
薄壁U型鋼混凝土梁火災(zāi)-結(jié)構(gòu)耦合的ANSYS分析
高軒能,黃文歡,張惠華
(華僑大學(xué)土木工程學(xué)院,福建泉州 362021)
采用在空間上運(yùn)用有限單元法與在時(shí)間上運(yùn)用有限差分法相結(jié)合的方法,建立組合梁的火-結(jié)構(gòu)耦合ANSYS分析模型;編制組合梁ANSYS熱分析程序和全過程受火非線性分析程序,對組合梁在國際標(biāo)準(zhǔn)ISO 834規(guī)定的標(biāo)準(zhǔn)火災(zāi)的溫度場和撓度-時(shí)間曲線進(jìn)行分析和計(jì)算.數(shù)值計(jì)算結(jié)果表明:在國際標(biāo)準(zhǔn) ISO 834規(guī)定的標(biāo)準(zhǔn)火災(zāi)下,試驗(yàn)梁測點(diǎn)的溫度與ANSYS理論計(jì)算結(jié)果高度吻合,試驗(yàn)梁位移-受火時(shí)間曲線與理論計(jì)算結(jié)果符合良好.
薄壁U型鋼;組合梁;火災(zāi);火-結(jié)構(gòu)耦合;ANSYS模擬
U型鋼-混凝土組合梁由冷彎薄壁U型鋼填充混凝土組成,是一種正在發(fā)展的新型構(gòu)件.U型鋼-混凝土組合梁能夠充分利用薄壁型鋼和混凝土兩種材料,外形整潔、施工快捷、無需模板、易于澆灌混凝土,并且具有降低樓層高度、節(jié)能省材、符合低碳經(jīng)濟(jì)發(fā)展趨勢等優(yōu)點(diǎn),應(yīng)用前景廣闊.國內(nèi)外學(xué)者[1-3]已有不少有關(guān)其在常溫下靜力承載性能的研究,但對于其在火災(zāi)高溫下的溫升特性及耐火性能,迄今研究不多[4-6].然而,火災(zāi)高溫下U型鋼-混凝土組合梁的傳熱,是在非穩(wěn)態(tài)環(huán)境下多介質(zhì)及復(fù)雜邊界條件下的非線性二維傳熱問題,不僅受到靜力荷載的作用,還要受到溫度荷載的作用,理論分析的難度很大.這是因?yàn)榻Y(jié)構(gòu)構(gòu)件在環(huán)境溫度變化的情況下形成了動(dòng)態(tài)的不均勻溫度場,高溫使材料的強(qiáng)度和變形性能嚴(yán)重劣化,又使結(jié)構(gòu)構(gòu)件產(chǎn)生劇烈的內(nèi)(應(yīng))力重分布;而且,溫度和荷載(應(yīng)力)有顯著的耦合作用效應(yīng),使材料的本構(gòu)關(guān)系和構(gòu)件的溫度-荷載路徑有較大變化.本文建立分析組合梁構(gòu)件的火-結(jié)構(gòu)耦合模型,了解高溫下U型鋼-混凝土組合梁的受力性能及其變化規(guī)律.
1.1 火-結(jié)構(gòu)耦合場
為便于計(jì)算,許多研究通常假設(shè)構(gòu)件均勻受熱,構(gòu)件截面溫度線性分布及構(gòu)件處于穩(wěn)態(tài)溫度場中,但這與構(gòu)件受火的實(shí)際工況不符,使得構(gòu)件的耐火性能分析結(jié)果誤差較大.要準(zhǔn)確模擬火災(zāi)條件下構(gòu)件的反應(yīng),就必須考慮溫度場及材料熱工性能和熱邊界條件隨時(shí)間的變化,即對構(gòu)件進(jìn)行瞬態(tài)非線性熱分析.為此,采用間接耦合法建立火-結(jié)構(gòu)耦合分析模型,利用ANSYS通用軟件對構(gòu)件內(nèi)力和變形進(jìn)行分析.即通過傳熱分析得到構(gòu)件隨時(shí)間變化的溫度場,再計(jì)算在此溫度場下構(gòu)件的受力,然后與靜力荷載疊加求解.具體流程如圖1所示.
1.2 傳熱分析及假定
設(shè)U型鋼-混凝土組合梁為室內(nèi)受火,可視為連續(xù)、各向同性、無內(nèi)熱源的無限長梁,為二維溫度場問題,傳熱方式按對流換熱和輻射換熱考慮.采用國際標(biāo)準(zhǔn)ISO 834規(guī)定的標(biāo)準(zhǔn)的溫度-時(shí)間曲線描述構(gòu)件的火災(zāi)環(huán)境[7],其升溫段(t≤th,th為火災(zāi)持續(xù)時(shí)間)表達(dá)式為
式(1)中:t為時(shí)間;θg為t時(shí)刻的空氣溫度;θg(0)為初始空氣溫度.
單位時(shí)間內(nèi),空氣向構(gòu)件對流傳遞的熱量為
圖1 ANSYS熱-結(jié)構(gòu)間接耦合法分析流程圖Fig.1 ANSYS analysis flow chart of thermalstructure indirect coup ling method
式(2)中:αc為對流傳熱系數(shù);θb為構(gòu)件表面(或保護(hù)層)溫度.
單位時(shí)間內(nèi),空氣向構(gòu)件輻射傳遞的凈熱量為
式(3)中:φ為形狀系數(shù);εr為綜合輻射系數(shù),εr=εfεm,εf與著火房間有關(guān),εm與構(gòu)件表面特性有關(guān);σ為Stefan-Boltzmann常數(shù),σ=56.7 nW· (m2·K4)-1.
構(gòu)件內(nèi)部主要以導(dǎo)熱方式傳遞熱量,其截面的導(dǎo)熱微分方程為
式(4)中:ρ為介質(zhì)密度;c為介質(zhì)比熱;l為介質(zhì)的導(dǎo)熱系數(shù);θ為截面(x,y)處在時(shí)刻t的溫度.
為便于分析計(jì)算,分析時(shí)做如下5個(gè)假定[8]:
(1)構(gòu)件截面的溫度場與材料的應(yīng)力、應(yīng)變及是否開裂無關(guān),即不考慮混凝土開裂或表層脫落后的截面局部變化引起的溫度重分布;
(2)構(gòu)件截面的溫度場是二維的,沿構(gòu)件的軸線不變;
(3)忽略鋼材與混凝土之間接觸熱阻對組合梁構(gòu)件傳熱的影響;
(4)忽略鋼筋體積對構(gòu)件截面溫度場的影響;
(5)一般情況下,不考慮水分蒸發(fā)對溫度場的影響.
1.3 參數(shù)選取
表1 熱分析系數(shù)Tab.1 Coefficients of thermal analysis
根據(jù)研究對象特性,在熱分析中將輻射傳熱與對流傳熱分開考慮,并參考已有文獻(xiàn)[9]的邊界條件假定和邊界參數(shù)取值.熱分析參數(shù)的具體取值,如表1所示.表1中:對流換熱系數(shù)為αc,綜合輻射系數(shù)為εres,形狀系數(shù)為φ.
圖2為一榀3層3跨連續(xù)框架,柱腳剛接,層高3 m,柱間距5.5 m,柱頂集中荷載P為151 kN,橫梁上的均布荷載q為25.4 kN·m-1.梁為工字型鋼截面H303.8 mm×165 mm×6.1 mm×10.2 mm,柱為工字型鋼截面H206.2 mm×203.9 mm×8.3 mm×12.5 mm[10].在底層左側(cè)防火間發(fā)生國際標(biāo)準(zhǔn)ISO 834規(guī)定的標(biāo)準(zhǔn)火災(zāi).防火間具有防止火災(zāi)蔓延和隔絕熱量傳遞的作用.
為了準(zhǔn)確模擬截面的非均勻溫度分布和捕捉局部屈曲等現(xiàn)象,防火間的梁和柱采用實(shí)體單元,其余部分采用梁單元,兩者接合處采用約束方程以保證變形協(xié)調(diào).圖3為梁和柱的熱邊界條件.柱內(nèi)部砌有磚墻,僅面向防火間的翼緣受火,梁上部有混凝土樓板,上翼緣上表面不受火.熱量以對流和輻射的形式從熱空氣傳遞到結(jié)構(gòu)表面,又以傳導(dǎo)的形式在結(jié)構(gòu)內(nèi)部傳播.防火間受火框架采用實(shí)體單元Solid 70,防火間外的框架梁柱采用梁單元Beam 188,梁柱節(jié)點(diǎn)自由度耦合及接合處截面約束,如圖4所示.變化曲線(負(fù)值表示向下凸曲),如圖5所示.從圖5中可以看出,鋼框架在受火13 min時(shí),梁跨中位移為118.36 mm;受火14 min時(shí),梁跨中位移急劇增加到190.14 mm.這表明框架在當(dāng)前荷載水平下,耐火極限可確定為14 min.另一方面,根據(jù)文獻(xiàn)[7]可知,在進(jìn)行整體結(jié)構(gòu)分析時(shí),可取跨中位移達(dá)到L/30為183.33 mm(L為框架梁的跨度)作為結(jié)構(gòu)構(gòu)件達(dá)到耐火極限的判據(jù).在計(jì)算中,位移值對應(yīng)的受火時(shí)刻為13.92 m in(約等于14 m in).由此可見,應(yīng)用ANSYS進(jìn)行火-結(jié)構(gòu)耦合分析不僅是可行的,而且具有良好的準(zhǔn)確性.
圖2 受火框架結(jié)構(gòu)(mm)Fig.2 Frame structure exposed to fire(mm)
圖3 梁柱受火面Fig.3 Boundaries exposed to fire of beam and column
在國際標(biāo)準(zhǔn)ISO 834規(guī)定的標(biāo)準(zhǔn)火災(zāi)作用下,框架梁跨中位移(y)隨受火時(shí)間(t)的ANSYS計(jì)算
圖5 防火間梁跨中位移隨受火時(shí)間的變化曲線 Fig.5 Curve between mid-span deflection and exposed-to-fire time of beam in fire compartment
圖4 梁柱節(jié)點(diǎn)單元耦合示意圖Fig.4 Schematic diagram of joint element coup ling
圖6 U型鋼-混凝土梁溫度測試點(diǎn)(mm)Fig.6 Temperature test points of U-section steel concrete beam(mm)
U型鋼-混凝土組合梁試件[11]在國際標(biāo)準(zhǔn) ISO 834規(guī)定的標(biāo)準(zhǔn)火災(zāi)爐中進(jìn)行受火試驗(yàn),其溫度測試點(diǎn)如圖6所示.組合梁試件參數(shù):U型鋼梁寬為150 mm;翼緣寬為50 mm;混凝土板寬為600 mm;混凝土板高為80 mm;U型鋼梁高為300 mm;鋼梁腹板及翼緣板厚均為3 mm.混凝土材料為硅質(zhì)混凝土,厚涂型防火涂料,厚度為6 mm,導(dǎo)熱系數(shù)為0.116 W· (m·K)-1,比熱為1.047 J·(g·K)-1,密度為(400± 20)kg·m-3.按上述方法建模,編制火-結(jié)構(gòu)耦合ANSYS分析程序進(jìn)行計(jì)算,得到的組合梁截面在30,90, 120和150 min時(shí)的溫度場,如圖7所示.
組合梁試件(UL 3)測點(diǎn)1的試驗(yàn)溫度與ANSYS溫度的比較,以及跨中位移-受火時(shí)間曲線,分別如圖8所示.從圖8中可以看出,理論計(jì)算結(jié)果與試驗(yàn)實(shí)測結(jié)果高度一致,特別是在測點(diǎn)1測得的溫度,其受火試驗(yàn)的全過程中都與ANSYS計(jì)算結(jié)果吻合,而且跨中位移-受火時(shí)間曲線也符合很好.這表明火-結(jié)構(gòu)耦合的ANSYS分析模型可用于結(jié)構(gòu)構(gòu)件受火全過程的數(shù)值模擬計(jì)算,計(jì)算結(jié)果具有良好的可靠性.
圖7 U型鋼-混凝土梁截面的溫度場分布(℃)Fig.7 Temperature fields of U-section steel concrete beam(℃)
圖8 受火U型鋼-混凝土梁理論與試驗(yàn)結(jié)果的比較Fig.8 Comparison between simulation results and test results of U-section steel-concrete beam under fire
通過分析和計(jì)算U型鋼-混凝土組合梁試件在國際標(biāo)準(zhǔn)ISO 834規(guī)定的標(biāo)準(zhǔn)火災(zāi)下截面的溫度場和梁的撓度-受火時(shí)間曲線,并將其與試驗(yàn)結(jié)果進(jìn)行比較驗(yàn)證,可以得到如下3點(diǎn)結(jié)論.
(1)在國際標(biāo)準(zhǔn)ISO 834規(guī)定的標(biāo)準(zhǔn)火災(zāi)下,試驗(yàn)梁測點(diǎn)的溫度與ANSYS理論計(jì)算結(jié)果在受火試驗(yàn)全過程高度吻合,試驗(yàn)梁跨中位移-受火時(shí)間曲線也與理論計(jì)算結(jié)果符合良好.表明,模擬分析薄壁型鋼混凝土組合結(jié)構(gòu)的耐火性能是可行的,所建立的火-結(jié)構(gòu)耦合ANSYS分析模型是合適的.
(2)火-結(jié)構(gòu)耦合ANSYS模擬計(jì)算結(jié)果與試驗(yàn)結(jié)果符合良好,應(yīng)用于結(jié)構(gòu)構(gòu)件耐火性能受火全過程的數(shù)值模擬計(jì)算不僅是可行性的,而且具有良好的準(zhǔn)確性和可靠性.
(3)應(yīng)用火-結(jié)構(gòu)耦合ANSYS分析模型,可對不同條件下U型鋼-混凝土組合梁的耐火性能進(jìn)行評判與分析,不僅節(jié)省了大量試驗(yàn)費(fèi)用,而且可為構(gòu)件的耐火試驗(yàn)方案提供依據(jù).
[1]張耀春,毛小勇,曹寶珠.輕鋼-混凝土組合梁的試驗(yàn)研究及非線性有限元分析[J].建筑結(jié)構(gòu)學(xué)報(bào),2003,24(1):26-33.
[2]周天華,何保康,李鑫全,等.帽型冷彎薄壁型鋼-混凝土組合梁的試驗(yàn)研究[J].建筑結(jié)構(gòu),2003,33(1):48-50,14.
[3]沈建華,高軒能,周期源.滑移對薄壁U型鋼-混凝土組合梁撓度的影響[J].華僑大學(xué)學(xué)報(bào):自然科學(xué)版,2009,30 (5):557-562.
[4]李國強(qiáng),吳波,韓林海.結(jié)構(gòu)抗火進(jìn)展與趨勢[J].建筑結(jié)構(gòu)進(jìn)展,2006,8(1):1-13.
[5]毛小勇,肖巖.標(biāo)準(zhǔn)升溫下輕鋼-混凝土組合梁的抗火性能研究[J].湖南大學(xué)學(xué)報(bào):自然科學(xué)版,2006,32(2):64-70.
[6]高軒能,陳明華,王進(jìn).薄壁槽型鋼-混凝土梁在標(biāo)準(zhǔn)火災(zāi)下的溫升特性[J].福州大學(xué)學(xué)報(bào):自然科學(xué)版,2008,36 (3):430-435.
[7]李國強(qiáng),韓林海,樓國彪,等.鋼結(jié)構(gòu)及鋼-混凝土組合結(jié)構(gòu)抗火設(shè)計(jì)[M].北京:中國建筑工業(yè)出版社,2006.
[8]陳明華.薄壁型鋼-混凝土組合梁抗火性能研究[D].泉州:華僑大學(xué),2007.
[9]KRUPPA J,ZHAO B.Fire resistance of composite beams to eurocode 4:Part 1.2[J].Journal of Constructional Steel Research,1995,33(1/2):51-69.
[10]張波,盛合太.ANSYS有限元數(shù)值分析原理與工程應(yīng)用[M].北京:清華大學(xué)出版社,2005.
[11]黃文歡.薄壁型鋼-混凝土組合梁抗火性能試驗(yàn)研究[D].泉州:華僑大學(xué),2008.
(責(zé)任編輯:錢筠英文審校:方德平)
ANSYSSimulation for Fire-Structure Coupling of Thin-Walled U-Section Steel-Concrete Beam
GAO Xuan-neng,HUANGWen-huan,ZHANG Hui-hua
(College of Civil Engineering,Huaqiao University,Quanzhou 362021,China)
To research the fire-resistant behavior of the composite beam in fire,an ANSYSmodel for fire-structure coup ling analysis was established,by the method combining finite element method in space and finite difference method in time,and the temperature fields and the deflection-time curvesof the composite beam s under ISO 834 firewere numerically simulated and analyzed with the thermal-analysis ANSYS program and the nonlinear w hole process analysis ANSYS program.The numerical results show:under ISO standard fire,the temperature of measuring point in tested beam is highly consistent with the ANSYS simulation result,and the curve between mid-span deflection and exposed-to-fire time of the tested beam is in good agreement with the theoretical result.
thin-walled U-section steel;composite beam;fire;fire-structure coup ling;ANSYS simulation
TU 352.5;TU 398.9
A
1000-5013(2011)03-0317-05
2010-09-29
高軒能(1962-),男,教授,主要從事鋼結(jié)構(gòu)穩(wěn)定及工程結(jié)構(gòu)災(zāi)害控制的研究.E-mail:gaoxn@hqu.edu.cn.
國家自然科學(xué)基金資助項(xiàng)目(51008133);福建省科技計(jì)劃重點(diǎn)項(xiàng)目(2006H0030);國務(wù)院僑辦科研基金資助項(xiàng)目(07QZR04);福建省廈門市科技計(jì)劃項(xiàng)目(3502Z20073036)