周 清,丁 杰,張新鑫,孫 亮
T形截面梁可與預(yù)制樓板組成疊合式構(gòu)件,因此被廣泛地應(yīng)用于裝配式結(jié)構(gòu)體系中。目前學(xué)者利用LS-DYNA[1-2]對于爆炸荷載作用下矩形截面梁的動力響應(yīng)及破壞機理進行了較為詳細(xì)的研究[3-4],而對T形截面鋼筋混凝土梁研究成果較少。在以往學(xué)者研究的基礎(chǔ)上,筆者利用LS-DYNA研究了配筋率、截面高度、翼緣尺寸等各因素對T形截面鋼筋混凝土梁抗爆性能的影響。
相關(guān)文獻(xiàn)表明[5],在所有的混凝土模型中,MAT _CONCRETE _DAMAGE最能有效地模擬混凝土在高應(yīng)變、大變形下的力學(xué)形態(tài)。另外,數(shù)值模擬結(jié)果對材料的參數(shù)非常敏感。因此,在有限元模型中,準(zhǔn)確定義合理的材料參數(shù)顯得尤為重要。*MAT _ CONCRETE _DAMAGE _REL3模型為*MAT_CONCRETE__DAMAGE的升級版本,在保留了后者優(yōu)點的同時,于模型參數(shù)的確定方面做了簡化,使用此模型模擬混凝土材料,用戶僅需要定義幾個常用的基本參數(shù)。筆者采用該模型模擬C40混凝土,其具體參數(shù)見表1。
表1 C40混凝土模型參數(shù)
采用彈塑性模型*MAT PLASTIC KINEMATIC模擬HRB400鋼筋,其具體參數(shù)見表2。
表2 HRB400鋼筋模型參數(shù)
爆炸荷載作用下,材料的應(yīng)變率可高達(dá)1 000/s,鋼筋與混凝土的強度會有較大程度的提高。在近距離抗爆數(shù)值模擬時通常采用動力增大系數(shù)(DIF)來考慮材料的應(yīng)變率效應(yīng)?;炷翉姸鹊腄IF采用K&C 模型[6],鋼筋強度的 DIF采用 C&S 模型[7]。
混凝土采用Solid 164單元,鋼筋采用Beam 161單元。根據(jù)眾多學(xué)者的數(shù)值模擬經(jīng)驗[8-12],爆炸超壓作用時間很短,可以忽略鋼筋與混凝土之間的相對滑移,故采用共用節(jié)點的方式定義兩者之間的接觸。鋼筋與混凝土有限元網(wǎng)格尺寸均為15 mm。
爆炸荷載作用下,混凝土?xí)l(fā)生損傷甚至破環(huán)。LS-DYNA中,各混凝土材料模型較難直觀地觀察到裂縫及破環(huán),筆者通過添加關(guān)鍵字*MAT_ADD_EROSION來模擬混凝土的破壞情況。
炸藥采用標(biāo)準(zhǔn)TNT,藥量取10 kg(當(dāng)爆炸距離為1 m時,產(chǎn)生的超壓峰值約為5.32 MPa)與20 kg(當(dāng)爆炸距離為1 m時,產(chǎn)生的超壓峰值約為9.77 MPa)兩種。利用 LS-DYNA關(guān)鍵字*LOAD_BLAST__ENHANCED定義任意藥量、任意位置空氣爆炸時標(biāo)準(zhǔn)TNT產(chǎn)生的爆炸荷載,利用關(guān)鍵字*LOAD SEGMEMT _SET定義梁頂面為受爆面[13-15]。利 用 LS - DYNA 關(guān) 鍵 字 *DEFINE_CURVE配合關(guān)鍵字*LOAD _BODY _Y定義梁豎直方向的重力,重力加速度取值9.81 m/s2。
T形截面鋼筋混凝土梁長度為3 m,截面高度為500 mm,截面寬度為200 mm,翼緣寬度為100 mm,翼緣厚度為100 mm,爆炸點位于跨中截面正上方1 m。
圖1 T形截面混凝土梁尺寸及爆炸點位置(單位:mm)
筆者按照以下原則制定各配筋方案:
(1)驗證縱筋配筋率對梁抗爆性能影響時,根據(jù)簡支梁特點,梁主要受力縱筋在下部。所以梁下部采用φ18、φ20、φ22、φ25 四種常用的受力縱筋,上部僅采用2φ14構(gòu)造鋼筋。梁箍筋采用φ8@100。
(2)驗證箍筋配筋率對梁抗爆性能影響時,保持梁下部受力縱筋3φ20不變,梁箍筋采用φ8@100、φ10@100、φ12@100三種常用的箍筋。其具體配筋方案見表3。
表3 縱筋配筋方案
梁整體振動結(jié)果見圖2~圖5及表4。
圖2 藥量 W=10 kg,不同縱筋配筋率的梁時間-豎向位移曲線
圖3 藥量 W=20 kg,不同縱筋配筋率的梁時間-豎向位移曲線
由圖2~圖5及表4分析可以知道:(1)梁在爆炸荷載用下的振動頻率均小于梁固有頻率;(2)增大縱筋配筋率、箍筋配筋率均可增大梁的振動頻率、提高梁的抗爆性能,并且藥量越大作用越明顯;(3)增大箍筋配筋率對梁抗爆性能的提高程度優(yōu)于縱筋配筋率。
圖4 藥量 W=10 kg,不同箍筋配筋率的梁時間-豎向位移曲線
圖5 藥量 W=20 kg,不同箍筋配筋率的梁時間-豎向位移曲線
表4 梁頻率與跨中最大撓度數(shù)據(jù)匯總表
表5為藥量W=20 kg荷載作用下各類鋼筋應(yīng)力匯總表。通過表5分析可以得到以下結(jié)論:(1)跨中翼緣箍筋應(yīng)力超過鋼筋屈服強度;(2)支座翼緣箍筋應(yīng)力沒有達(dá)到屈服強度;(3)跨中縱筋利用率隨著縱筋配筋率增大而減小且沒有超過屈服強度。
表5 鋼筋應(yīng)力匯總表
圖6為藥量W=20 kg,縱筋配筋率、箍筋配筋率最大時的梁破壞圖。通過分析可以得到以下結(jié)論:
(1)縱筋配筋與箍筋配筋率最大時,梁破壞區(qū)主要發(fā)生在跨中梁底與支座處??缰辛旱灼茐妮^為嚴(yán)重,有部分混凝土脫離梁主體。翼緣破壞程度較小。
(2)跨中梁底截面混凝土以彎曲破壞為主,而支座處混凝土以剪切破壞為主。
圖6 縱筋配筋率與箍筋配筋率最大時的梁破壞圖
通過上文研究發(fā)現(xiàn),梁的抗爆性能在配筋率最大時表現(xiàn)最明顯。所以,研究截面高度對梁抗爆性能影響時,配筋方案分別取縱筋面積最大的配筋方案4與箍筋面積最大的配筋方案6。梁截面高度取500 mm、600 mm、700 mm共3種,其在2中配筋方案下的6種截面形式如表6所示。
表6 不同梁高、最大縱筋與箍筋配筋率下,梁的6種截面
圖7、圖8分別為藥量 W=20 kg,縱筋配筋率、箍筋配筋率最大時,不同梁高的梁整體時間-豎向位移曲線。通過分析可以得到以下結(jié)論:(1)梁在爆炸荷載用下的振動頻率均小于梁固有頻率。(2)縱筋、箍筋不變時,采用截面3、截面6的梁振動頻率較采用截面1、截面4明顯增加,剛度明顯增強(見表 7)。
圖7 藥量 W=20 kg、縱筋配筋率最大、不同梁高的梁時間-豎向位移曲線
圖8 藥量 W=20 kg,箍筋配筋率最大、不同梁高的梁時間-豎向位移曲線
表7 梁頻率與跨中最大撓度數(shù)據(jù)匯總表
通過對梁鋼筋應(yīng)力匯總表表8分析可以得到以下結(jié)論:(1)跨中翼緣箍筋應(yīng)力超過鋼筋屈服強度;(2)支座翼緣箍筋應(yīng)力沒有達(dá)到屈服強度;(3)跨中縱筋應(yīng)力沒有達(dá)到屈服強度,鋼筋應(yīng)力利用率隨著梁高增大而減小;(4)支座箍筋均超過屈服強度。
表8 梁鋼筋應(yīng)力匯總表
圖9為藥量 W=20 kg,縱筋配筋率、箍筋配筋率最大時,h=700 mm的梁破壞圖。通過分析可以得到以下結(jié)論:
(1)梁破壞區(qū)主要分布在跨中梁底與支座。與圖6相比,跨中梁底混凝土破壞程度變小。
(2)梁跨中截面以彎曲破壞為主,梁支座截面以剪切破壞為主。
圖9 截面高度 h=700 mm,縱筋與箍筋配筋率最大時的梁破壞圖
與上文所采用的方法相同,研究翼緣寬度對梁抗爆性能影響時,配筋方案分別取縱筋面積最大的配筋方案4與箍筋面積最大的配筋方案6。保持梁高度500 mm不變,梁翼緣截面寬度取100 mm、150 mm、200 mm共3種,其在2中配筋方案下的6種截面形式如表9所示。
表9 不同翼緣截面寬度、最大縱筋與箍筋配筋率下,梁6種截面形式
4.1.1 梁整體振動分析
梁翼緣寬度增大會對梁的抗爆性能產(chǎn)生兩方面的影響:(1)使梁接觸爆炸荷載的面積增大;(2)使梁整體剛度增強。梁抗爆性能應(yīng)為上述兩方面影響的綜合結(jié)果。
圖10、圖11分別為藥量 W=20 kg,縱筋配筋率、箍筋配筋率最大時,不同翼緣寬度梁的整體時間-豎向位移曲線。
圖10 藥量 W=20 kg,縱筋配筋率最大、不同翼緣寬度的梁時間-豎向位移曲線
圖11 藥量 W=20 kg,箍筋配筋率最大、不同翼緣寬度的梁時間-豎向位移曲線
通過分析可以得到以下結(jié)論:(1)梁在爆炸荷載用下的振動頻率均小于梁固有頻率;(2)縱筋、箍筋配筋率不變時,梁豎向撓度與振動頻率隨著翼緣寬度的增大而減?。ㄒ姳?0)。
以上說明,梁翼緣寬度越大,梁振動頻率越小,破壞程度越大。
4.1.2 鋼筋應(yīng)力分析
通過對梁鋼筋應(yīng)力匯總表11分析可以得到以下結(jié)論:(1)跨中翼緣箍筋應(yīng)力超過鋼筋屈服強度;(2)支座翼緣箍筋應(yīng)力隨著梁翼緣寬度的增大而增大;(3)梁跨中縱筋利用率隨著翼緣寬度增大而增大;(4)梁支座箍筋均超過屈服強度。
表10 梁頻率與跨中最大撓度數(shù)據(jù)匯總表
表11 梁鋼筋應(yīng)力匯總表
4.1.3 混凝土破壞程度分析
圖12表示藥量 W=20 kg,縱筋配筋率、箍筋配筋率最大時,翼緣寬度s=200 mm的梁破壞圖,通過分析可以得到以下結(jié)論:
(1)縱筋配筋率最大時,梁發(fā)生嚴(yán)重破壞??缰辛旱谆炷链竺娣e脫落,支座處出現(xiàn)大量裂縫。梁上部翼緣破壞嚴(yán)重,翼緣沿其與腹板的交線幾乎發(fā)生折斷現(xiàn)象。
(2)箍筋配筋率最大時,梁同樣發(fā)生嚴(yán)重破壞,但梁破壞程度較縱筋配筋率最大時的情況小。
圖12 翼緣寬度 s=200 mm,縱筋與箍筋配筋率最大時的梁破壞圖
與上文所采用的方法相同,研究翼緣厚度對梁抗爆性能影響時,配筋方案分別取縱筋面積最大的配筋方案4與箍筋面積最大的配筋方案6。梁翼緣厚度取100 mm、150 mm、200 mm共3種,其在2中配筋方案下的6種截面形式如表12所示。
表12 不同翼緣厚度、最大縱筋與箍筋配筋率下,梁的6種截面形式
4.2.1 梁整體振動分析
圖13、圖14分別為藥量 W=20 kg,縱筋配筋率、箍筋配筋率最大時,不同翼緣厚度的梁整體時間-豎向位移曲線。通過分析可以得到以下結(jié)論:(1)縱筋、箍筋配筋率不變時,梁振動頻率隨著翼緣厚度的增大而增大;(2)縱筋配筋率不變時,采用方案1的梁跨中撓度約為20 mm,而采用方案3梁跨中截面撓度約為10 mm,并且梁頻率明顯增大。箍筋配筋率不變時,采用方案4的梁跨中撓度約為10 mm,而采用方案6梁跨中截面撓度約為7 mm(見表13)。
圖13 藥量 W=20 kg、縱筋配筋率最大、不同翼緣厚度的梁整體時間-豎向位移曲線
圖14 藥量 W=20 kg,箍筋配筋率最大、不同翼緣厚度的梁整體時間-豎向位移曲線
表13 梁頻率與跨中最大撓度數(shù)據(jù)匯總表
以上說明,在保持最大縱筋配筋率、箍筋配筋率不變的條件下,增大翼緣厚度可以有效地提高梁抗爆性能。
4.2.2 鋼筋應(yīng)力分析
通過對梁鋼筋應(yīng)力匯總表14分析可以得到以下結(jié)論:
(1)跨中翼緣箍筋應(yīng)力超過鋼筋屈服強度。隨著翼緣厚度增大鋼筋應(yīng)力減小。
(2)支座翼緣箍筋應(yīng)力超過屈服強度,隨著梁翼緣厚度的增大變化不大,在550 MPa上下浮動。
(3)跨中縱筋利用率隨著翼緣厚度增大而增大。
(4)支座箍筋應(yīng)力均超過屈服強度。各工況下箍筋應(yīng)力變化不大,約為600 MPa。
表14 梁鋼筋應(yīng)力匯總表
4.2.3 混凝土破壞程度分析
圖15為藥量 W=20 kg,縱筋配筋率、箍筋配筋率最大時,翼緣厚度h1=200 mm的梁破壞圖與混凝土塑性應(yīng)變分布圖。通過分析可以得到以下結(jié)論:梁破壞區(qū)主要分布在跨中梁底、梁側(cè)與支座翼緣處??缰幸砭壊糠只炷涟l(fā)生破壞。與增加翼緣寬度相反,增加翼緣厚度可以降低梁頂面混凝土破壞程度。
圖15 翼緣厚度 h1=200 mm、縱筋與箍筋配筋率最大時的梁破壞圖
利用LS-DYNA軟件,對影響簡支T形截面鋼筋混凝土梁抗爆性能的因素進行了詳細(xì)的分析。通過研究發(fā)現(xiàn):
(1)增大縱筋配筋率與箍筋配筋率均能提高梁的抗爆性能。相同條件下,增大箍筋配筋率較增大縱筋配筋率效果更好。
(2)增大梁的截面高度會較大程度地減小爆炸荷載作用下梁跨中截面的撓度,使梁振動頻率變大、剛度增加,有效地提高梁的抗爆性能。
(3)增大翼緣寬度雖然一定程度上增加了梁的整體剛度,但同時也增加了爆炸荷載作用面積。綜合分析,增大翼緣寬度不利于提高梁的抗爆性能。增大翼緣厚度可以有效地減小爆炸荷載作用下梁跨中截面的豎向撓度,增加梁的振動頻率,提高梁的抗爆性能。同時,增加翼緣的厚度可以使得梁頂截面混凝土破壞程度降低。
(4)跨中梁底縱筋應(yīng)力一般均小于鋼筋屈服強度,而支座箍筋應(yīng)力均超過鋼筋屈服強度且在各種工況下應(yīng)力值較為穩(wěn)定。
(5)“增加箍筋配筋率+增加梁截面高度+增加翼緣厚度”為提高T形截面鋼筋混凝土簡支梁抗爆性能的最優(yōu)組合形式。
[1] 劉成清,何 斌,陳 弛,等.ANSYS/LS-DYNA工程結(jié)構(gòu)抗震、抗撞擊與連續(xù)倒塌分析[M].北京:中國建筑工業(yè)出版社,2014.
[2] 唐長剛.LS-DYNA有限元分析及仿真[M].北京:電子工業(yè)出版社,2014.
[3] 李猛深,李 杰,李 宏,等.爆炸荷載下鋼筋混凝土梁的變形和破壞[J].爆炸與沖擊,2015,35(2):177-183.
[4] 都 浩,鄧 艸凡,杜榮強.爆炸荷載作用下鋼筋混凝土梁動力響應(yīng)的數(shù)值分析[J].山東科技大學(xué)學(xué)報(自然科學(xué)版),2010,29(6):50-54.
[5] 李忠獻(xiàn),師燕超.建筑結(jié)構(gòu)抗爆分析理論[M].北京:科學(xué)出版社,2015.
[6] Bischoff P H,Perry S H.Compressive behavior of concrete at high stain rate[J].Material and Structure,1991,24(144):425-450.
[7] Malvar L.Review of static and dynamic properties of steel reinforcing bars[J].ACI Material Journal,1998,95(5):609-616.
[8] 田 力,付協(xié)偉.爆炸荷載作用下高層鋼筋混凝土框架結(jié)構(gòu)的連續(xù)倒塌機理研究[J].科學(xué)技術(shù)與工程,2016,16(3):232-238.
[9] 師燕超,李忠獻(xiàn),郝 洪.爆炸荷載作用下鋼筋混凝土框架的連續(xù)倒塌分析[J].解放軍理工大學(xué)學(xué)報,2007,8(6):652-658.
[10] 周繼凱,莊天文,慕建磊,等.爆炸荷載下鋼筋混凝土框架連續(xù)倒坍特性及影響因素研究[J].混凝土,2015,5(4):5-13.
[11] 周 清,牛寧寧,劉海員,等.剪跨比對爆炸荷載作用下鋼筋混凝土梁抗爆性能的影響[J].水利與建筑工程學(xué)報,2017,15(1):219-224.
[12] 謝 欣,趙 文,關(guān)永平.地鐵車站STS工法鋼管混凝土橫向連接正截面受彎數(shù)值模擬[J].水利與建筑工程學(xué)報,2014,12(6):74-77.
[13] 宋 力,丁 慧,王建平.高強箍筋約束高強混凝土方柱軸壓性能數(shù)值試驗研究[J].水利與建筑工程學(xué)報,2015,13(2):57-61.
[14] 龔順風(fēng),鄧 歡,朱升波,等.近爆作用下鋼筋混凝土板動態(tài)破壞的數(shù)值模擬研究[J].振動與沖擊,2012,31(2):20-24.
[15] 李忠獻(xiàn),師燕超,史祥生.爆炸荷載作用下鋼筋混凝土板破壞評定方法[J].建筑結(jié)構(gòu)學(xué)報,2009,30(6):60-66.