落錘沖擊下鋼筋混凝土梁動態(tài)響應(yīng)數(shù)值模擬分析*

李青山,路國運(yùn)

(1. 山西瀟河建筑產(chǎn)業(yè)有限公司, 山西 太原 030000; 2. 太原理工大學(xué), 山西 太原 030024)

0 引言

鋼筋混凝土結(jié)構(gòu)在服役過程中往往會遭受各種不可預(yù)知的沖擊荷載,如汽車碰撞、落石沖擊、爆炸、水流沖擊等。對于有特殊要求的建筑結(jié)構(gòu),在進(jìn)行結(jié)構(gòu)設(shè)計時,其抗沖擊性能需加以考慮。研究并深入了解鋼筋在沖擊過程中的作用,對于提高結(jié)構(gòu)抗沖擊性能、優(yōu)化結(jié)構(gòu)設(shè)計具有重要工程意義。在鋼筋混凝土抗沖擊性能測試中,許多學(xué)者采用落錘沖擊試驗方法對錘頭形狀[1]、沖擊速度和能量[2]、試件慣性效應(yīng)[3]、配筋率[4]等進(jìn)行分析。結(jié)果表明,可采用任意形狀錘頭對鋼筋混凝土梁進(jìn)行沖擊試驗且結(jié)構(gòu)慣性力占沖擊力的2/3以上[5],因而許多研究中采用支反力作為評估鋼筋混凝土梁在破壞過程中受到的真實彎曲荷載。在沖擊荷載下,配筋率對鋼筋混凝土梁斷裂性能和鋼筋在裂縫擴(kuò)展過程中的限裂作用具有一定影響[6],且試件最大承載力由配筋率決定[7]。試件制備過程中,不同學(xué)者利用鋼筋和混凝土材料性能具有一定差異,且落錘沖擊試驗中不同技術(shù)人員、試驗設(shè)備、場地、數(shù)據(jù)采集和加載條件的控制存在誤差導(dǎo)致試驗結(jié)果具有不同程度離散性。隨著鋼筋混凝土結(jié)構(gòu)抗沖擊性能不斷深入研究,學(xué)者們對特定試驗環(huán)境和試件進(jìn)行試驗后得到的試驗結(jié)果均具有局限性。采用試驗方法可直觀獲得鋼筋混凝土結(jié)構(gòu)抗沖擊性能參數(shù),但過高的試驗成本和誤差致使其可靠性降低,試驗結(jié)果得不到有效推廣。

隨著計算機(jī)軟件開發(fā)和有限元技術(shù)在工程領(lǐng)域的廣泛應(yīng)用,數(shù)值仿真計算為工程結(jié)構(gòu)受力分析提供了有效參考依據(jù)。李敏等[8]利用ABAQUS模擬了地震荷載作用下內(nèi)加載速率對鋼筋混凝土梁力學(xué)行為的影響。Ahmed[9]采用ABAQUS對梁在沖擊荷載作用下的動力特性進(jìn)行了研究,并與已發(fā)表的試驗結(jié)果進(jìn)行對比后對有限元模型結(jié)構(gòu)響應(yīng)進(jìn)行了分析。Meng等[10]探討了鋼筋強(qiáng)度、混凝土強(qiáng)度等級和鋼筋保護(hù)層厚度對梁抗沖擊性能的影響,數(shù)值模擬結(jié)果表明:隨著鋼筋強(qiáng)度提高,梁承載力提高,但增加幅度會隨著鋼筋強(qiáng)度的提高而逐漸減小。Jiang等[11]模擬了鋼筋混凝土梁在動態(tài)沖擊下的結(jié)構(gòu)響應(yīng)。有限元建模存在本構(gòu)模型不同、網(wǎng)格劃分精細(xì)程度不一等差別,與試驗結(jié)果進(jìn)行對比后可驗證其有效性。本文基于ABAQUS有限元軟件,模擬鋼筋混凝土梁在沖擊荷載下的動態(tài)響應(yīng)。通過與試驗結(jié)果對比,驗證了模型有效性?；谠撃Ｐ?討論了配筋率一定的情況下不同黏結(jié)面積鋼筋混凝土梁承載力、抗變形能力和破壞模式。

1 有限元模型

曾翔等[12]利用落錘對無腹筋的鋼筋混凝土梁進(jìn)行三點(diǎn)彎曲試驗,探討了落錘重量、沖擊速度和沖擊能量對鋼筋混凝土梁抗沖擊行為的影響。本文對該試驗進(jìn)行有限元建模,所用模型尺寸與參考文獻(xiàn)一致,如圖1所示[13]?；炷量箟簭?qiáng)度為45.8MPa,試驗用混凝土梁長3 000mm,跨度為 1 860mm, 截面尺寸為150mm×310mm。在梁底部配置縱筋且不考慮箍筋配置情況?？v向鋼筋為3根φ16二級螺紋鋼筋(HRB335),保護(hù)層厚度為25mm。在對試件進(jìn)行澆筑前,抽取3根鋼筋進(jìn)行單軸拉伸試驗。由試驗結(jié)果可知鋼筋屈服強(qiáng)度為394MPa,極限強(qiáng)度為525MPa,彈性模量為1.94×105MPa。

圖1 梁尺寸及配筋截面

表1 混凝土模型參數(shù)[13]

圖2 有限元模型

鋼筋采用三折線模型,有限元模型中網(wǎng)格單元為三維8節(jié)點(diǎn)縮減積分單元(C3D8R),網(wǎng)格尺寸為5mm×5mm×5mm。由試驗結(jié)果可知,不同沖擊能量下鋼筋混凝土梁表面出現(xiàn)裂紋后混凝土與鋼筋間的黏結(jié)被破壞,部分混凝土出現(xiàn)崩落且鋼筋產(chǎn)生大變形甚至達(dá)到塑性變形。裂紋以外混凝土和鋼筋間的黏結(jié)完好,鋼筋有效保持了混凝土試件在沖擊破壞后的完整性。因此,有限元模型中可默認(rèn)沖擊荷載下鋼筋與混凝土間的完美黏結(jié),不考慮鋼筋拔出效應(yīng)。在有限元模型中鋼筋以Embed形式嵌入混凝土基體中。試驗表明[14],高應(yīng)變率作用下建筑鋼筋彈性模量保持不變,鋼筋HPB235彈性模量為2.1×105MPa。鋼筋屈服強(qiáng)度和極限強(qiáng)度隨應(yīng)變率增大而提高,變化規(guī)律可分別用式(1),(2)表示:

(1)

(2)

2 有限元模型驗證

將模擬得到的破壞模式與試驗結(jié)果進(jìn)行對比,如圖3所示。試驗中通過高速攝影記錄了鋼筋混凝土梁在不同沖擊速度下變形、損傷及破壞狀況。圖3a表示沖擊速度為5.71m/s時梁表面最終損傷破壞情況,為對試驗拍攝照片處理后得到的裂紋分布圖,圖3b為相對應(yīng)沖擊速度下梁拉伸損傷云圖。由圖3可知,配筋率為1.51%的鋼筋混凝土梁在沖擊荷載下出現(xiàn)剪切裂紋,混凝土梁上表面受沖擊區(qū)域發(fā)生壓潰破壞。此外在支座處出現(xiàn)豎直裂紋,這是由于沖擊初始時刻剪跨比為3.47的鋼筋混凝土梁主要發(fā)生剪切破壞;隨著剪切裂紋產(chǎn)生,混凝土破壞且內(nèi)部能量得到釋放,多個部分混凝土逐漸斷裂。由于鋼筋與混凝土間黏結(jié)良好,試件保持了較好完整性。此時跨度以外混凝土梁與已破壞部分自重導(dǎo)致支座處彎矩迅速增大,彎曲裂紋萌生。有限元模型中梁體跨中附近混凝土損傷嚴(yán)重,且支座豎直方向彎曲裂紋明顯。由于底部鋼筋彈性模量遠(yuǎn)大于混凝土,彎曲裂紋擴(kuò)展至鋼筋附近后停止繼續(xù)擴(kuò)展,有限元模型中左支座附近混凝土損傷情況與試驗結(jié)果吻合。因此,采用該模型可有效模擬鋼凝土梁遭受沖擊形態(tài)后的破壞形態(tài)。

圖3 梁損傷破壞狀況

由圖4可知,沖擊速度為5.71m/s時沖擊力時程曲線與數(shù)值模擬結(jié)果吻合較好。試驗中隨著錘頭與試件上表面的接觸,安裝在錘頭與配重之間的壓力傳感器信號迅速增大至最大值。試驗結(jié)果中沖擊力信號觸發(fā)后約0.25ms后沖擊力達(dá)最大值296kN;數(shù)值模擬沖擊力峰值距觸發(fā)時間為0.21ms,峰值為298.62kN。試驗中沖擊力信號在4～12ms區(qū)間內(nèi)震蕩明顯,這是由于混凝土中裂紋不斷擴(kuò)展導(dǎo)致試件內(nèi)部發(fā)生破壞,錘頭與梁體之間反復(fù)碰撞接觸和分離,沖擊力信號呈震蕩逐漸減小趨勢。有限元模型中混凝土單元中損傷單元數(shù)量逐漸加劇且被刪除。試驗中混凝土盡管發(fā)生破壞并存在部分混凝土碎塊崩落現(xiàn)象,但梁體內(nèi)部混凝土不斷擠壓并未完全離開梁體。試驗與有限元模擬自身的差異性致使試驗結(jié)果與模擬結(jié)果無法完全吻合,這對于評估鋼筋混凝土梁抗沖擊性能影響較小,因此本文對沖擊力峰值后的曲線軟化部分暫不做具體分析。對比可知,本文建立的數(shù)值模型可有效模擬文獻(xiàn)中鋼筋混凝土梁破壞過程及沖擊力響應(yīng)。

圖4 沖擊力時程曲線

3 黏結(jié)面積對梁承載力和破壞模式的影響

配筋率是工程建設(shè)中對于鋼筋混凝土結(jié)構(gòu)常用力學(xué)性能參數(shù),一般指面積配筋率。鋼筋混凝土結(jié)構(gòu)中縱向受力鋼筋面積與結(jié)構(gòu)橫截面積有效面積之比,即受拉鋼筋面積與主梁面積之比,計算公式為:p=As/A,其中As為縱向受力鋼筋截面面積,A為主梁面積。試驗中配筋率為1.51%。

基于本文模型,通過調(diào)整單根鋼筋直徑及結(jié)構(gòu)內(nèi)部鋼筋數(shù)量得到配筋率一定時鋼筋與混凝土黏結(jié)面積不同的試件。設(shè)置鋼筋直徑分別為8,16,27.72mm,分別對應(yīng)鋼筋數(shù)量為12,3,1根,其黏結(jié)面積由S=nπdL得到,其中n為鋼筋數(shù)量,d為鋼筋直徑,L為鋼筋長度(鋼筋長度均為3 000mm)。對應(yīng)黏結(jié)面積分別為0.904,0.452,0.261m2,分別標(biāo)記為S1,S2,S3?？梢娕浣盥室欢〞r鋼筋面積減小,鋼筋數(shù)量增多可有效提高鋼筋與混凝土間黏結(jié)面積。

不同黏結(jié)面積鋼筋混凝土梁在沖擊速度為5.71m/s時的沖擊力時程曲線如圖5所示。模擬結(jié)果表明,沖擊速度一定時,不同黏結(jié)面積試件沖擊力變化趨勢相同,最大沖擊力隨著黏結(jié)面積減小而增大,且沖擊力達(dá)到峰值時間隨黏結(jié)面積的增大而減小。黏結(jié)面積為0.452,0.261m2時,沖擊力峰值較接近?？梢婐そY(jié)面積減小但鋼筋直徑增大可提高鋼筋混凝土梁抗沖擊性能。

圖5 不同黏結(jié)面積試件沖擊力時程曲線

不同黏結(jié)面積鋼筋混凝土梁在5.71m/s沖擊速度下的跨中撓度時程曲線如圖6所示。隨著黏結(jié)面積減小,鋼筋混凝土梁跨中撓度增長速度減緩且最大值減小。這是由于配筋率一定時鋼筋與混凝土黏結(jié)面積減小的同時鋼筋直徑增大。相同彎矩下直徑較大鋼筋抗變形能力更強(qiáng)。因此,同一沖擊速度下,相同配筋率鋼筋混凝土梁抗變形能力隨黏結(jié)面積減小而降低。在工程實際設(shè)計中,對于已確定配筋率的構(gòu)件,調(diào)整鋼筋直徑和數(shù)量可增大鋼筋與混凝土黏結(jié)面積,但由于單根鋼筋直徑減小過多,試件抗變形能力減弱,沖擊荷載下梁體變形更大,更容易遭到破壞。

圖6 不同黏結(jié)面積試件跨中撓度時程曲線

準(zhǔn)靜態(tài)加載過程中,施加在試件上的荷載全部用于試件破壞,因而支座反力與外荷載大小相同、方向相反。沖擊荷載下試件破壞過程與靜態(tài)破壞過程存在差異。落錘沖擊過程中,已有試驗結(jié)果表明一部分沖擊力用于試件加速運(yùn)動,另一部分沖擊力用于試件破壞。與靜態(tài)荷載下試件的破壞過程不同,沖擊力并不是試件破壞時受到的真實彎曲荷載。在動態(tài)三點(diǎn)彎曲過程中,將支反力作為試件破壞過程中受到的真實彎曲荷載比沖擊力更有說服性。因此,提取支座反力進(jìn)行對比,如圖7所示。由圖7可知,黏結(jié)面積越大,支反力峰值越大。

圖7 不同黏結(jié)面積試件支反力時程曲線

沖擊過程中,各組分耗能無法通過試驗方法獲得,也就無法明確混凝土和鋼筋在試件破壞過程中的能量耗散情況。而鋼筋混凝土結(jié)構(gòu)中不同組分耗能占比對于合理設(shè)計并優(yōu)化結(jié)構(gòu)具有重要工程意義。數(shù)值模型中鋼筋材料為彈塑性模型且不發(fā)生損傷,混凝土采用損傷塑性模型且混凝土單元會出現(xiàn)損傷。在模型中可通過提取混凝土損傷耗能時程曲線來對比分析鋼筋在沖擊過程中的作用。如圖8所示,不同鋼筋混凝土試件在落錘沖擊過程中混凝土損傷耗能呈相似增長趨勢。在沖擊初始前期迅速增大;當(dāng)應(yīng)力傳遞至鋼筋附近時,混凝土損傷減弱;隨著鋼筋變形加劇,混凝土內(nèi)部裂紋繼續(xù)擴(kuò)展,損傷耗能再次快速增長直至最大值,此時試件完全破壞。對比不同黏結(jié)面積鋼筋混凝土試件損傷耗能時程曲線可知,鋼筋直徑減小、黏結(jié)面積增大,導(dǎo)致混凝土應(yīng)力傳播范圍增大,在短暫沖擊時間內(nèi)大量混凝土單元發(fā)生變形。而應(yīng)力分散減弱了內(nèi)部局部單元損傷情況,導(dǎo)致整體損傷耗能減小。

圖8 不同黏結(jié)面積試件損傷耗能時程曲線

4 結(jié)語

本文利用ABAQUS建立鋼筋混凝土梁落錘沖擊有限元模型,通過與文獻(xiàn)中試驗結(jié)果對比后驗證了模型有效性。利用該模型模擬相同沖擊速度下不同鋼筋與混凝土黏結(jié)面積鋼筋混凝土梁動態(tài)結(jié)構(gòu)響應(yīng),詳細(xì)討論了黏結(jié)面積對結(jié)構(gòu)抗沖擊性能、變形能力和混凝土損傷耗能情況。得到以下結(jié)論。

1)配筋率一定時,減小鋼筋直徑、提高鋼筋數(shù)量可促使鋼筋與混凝土間黏結(jié)面積增大。黏結(jié)面積的增大導(dǎo)致沖擊力峰值減小,試件抗變形能力減弱。

2)支反力可作為評價沖擊過程中試件受到的真實彎曲荷載,且支反力隨著黏結(jié)面積減小而降低,試件破壞程度減弱。

3)鋼筋與混凝土間黏結(jié)面積的增大導(dǎo)致試件內(nèi)部應(yīng)力傳播更為廣泛,混凝土損傷單元減少,試件損傷耗能減弱,保護(hù)了試件完整性。