近爆作用下鋼筋混凝土板動態破壞的數值模擬研究

龔順風，鄧 歡，朱升波，金偉良

（浙江大學 結構工程研究所，杭州 310058）

1 數值模型

在數值模擬中，爆炸沖擊波首先在如圖1所示的1D楔形空氣中傳播，當沖擊波波陣面到達距離爆心0.5 m處，即相當于接近鋼筋混凝土板上表面處，將其映射到3D模型中計算其與鋼筋混凝土板的相互作用。板的凈尺寸為1.0 m ×2.0 m ×0.1 m，鋼筋直徑為8 mm，長邊方向鋼筋間距為100 mm，短邊方向鋼筋間距為120 mm，雙層雙向布置，如圖2所示。炸藥置于距板頂面中心0.5 m高度處。鋼筋混凝土板的兩個長邊自由，短邊采用簡支約束，以此模擬單向板的工作狀態。鋼筋和混凝土均采用Lagrange網格進行建模，其中混凝土單元數為464 208，鋼筋為9 872?？諝獠捎肊uler-FCT算法，邊界面設為物質流出。在爆炸荷載作用下，由于作用的時間非常短，通常可以假定鋼筋和混凝土之間的粘結是完好的，同時混凝土單元的大變形會導致網格的扭曲，為解決網格的大變形問題，模型中采用侵蝕算法（Erosion）處理，這種算法將失效的混凝土單元自動刪除，避免了網格扭曲造成的計算精度下降、計算步長變小等問題。

圖1 TNT與空氣的一維楔形模型Fig.1 Geometry of 1D wedge with axial symmetry

圖2 鋼筋布置Fig.2 Rebar arrangement

2 材料模型

2.1 空氣和炸藥的狀態方程

通常將空氣假定為理想氣體，采用線性多項式狀態方程來描述爆炸后空氣的性能。根據Gama準則，空氣的狀態方程可表示為［9－10］：

式中：pa為氣體壓力；γ為比熱比，理想氣體取值為1.4；初始時刻的空氣密度ρa=1.29 kg/m3；E為空氣單位體積的內能，其初始值E0=2.5×105J/m3。對于不同炸藥爆炸的數值模擬，一般采用經驗的JWL狀態方程來描述炸藥爆轟過程壓力和內能及相對體積的關系［9，10］：

式中：p為爆炸產生的壓力，V為單位體積裝藥產生的爆轟產物的體積，即當壓力為p時的體積與初始體積的比值，ψ是炸藥的比內能；A、B、R1、R2和ω是由相關爆炸實驗確定的材料參數。對于TNT炸藥，初始密度ρ0=1 630 kg/m3，單位體積的初始內能 ψ0=6.0 ×103MJ/m3，其他相關參數見表1。

表1 炸藥的材料參數Tab.1 Material parameters of explosive

2.2 鋼筋

鋼筋采用J－C模型［11］，該模型適合模擬材料在大應變、高應變率和高溫情況下的力學性能。在爆炸荷載作用下，熱傳導的時間遠大于爆轟持續的時間，且由于變形熱導致的溫升并不大，材料出現強化繼而接近失效時，應力減小的幅度比較有限，所以溫度軟化效應并不明顯。爆炸沖擊荷載作用下混凝土結構中鋼筋的應變率可以達到 50 s－1～100 s－1［11－12］，圖 5 給出鋼筋在不同應變率下的應力應變關系曲線，該材料本構模型可表示為：

式中：εp為等效塑性應變為等效塑性應變率為參考塑性應變率，取=1.0 s－1。A表示材料的屈服強度是對材料強化段的描述，C為應變率敏感參數。鋼筋的失效應變取為0.1，其他相關參數見表2。

測試關注網絡參數N,節點參數Rs,Rc,以及覆蓋百分比p對無線傳感器生命期的影響,即關注在上述參數取值不同的情況下,節點所求得的覆蓋集情況以及算法的性能表現。

表2 鋼筋的材料參數Tab.2 Material parameters of reinforcing bar

2.3 混凝土

混凝土采用AUTODYN軟件材料庫中的CONC－35 MPA材料模型，該模型包括P－α狀態方程，RHT［12］強度模型和損傷模型以及侵蝕算法。RHT強度模型通過引進三個不同的強度面可以較好地表示混凝土材料具有屈服強度、最大強度及殘余強度的特性，能較為合理地描述混凝土從彈性到失效的整個過程，是當前比較適合用來模擬混凝土材料在爆炸沖擊荷載作用下動力特性的模型。

2.3.1 三個強度面

最大強度面定義為靜水壓力、偏應力面上羅德角和應變率的函數，可表示為：

其中，受壓子午線為：

式中：ψ=rt/rc，rt為靜水壓力軸到拉伸子午線的距離，rc為靜水壓力軸到受壓子午線的距離，羅德角θ是第二和第三偏應力不變量的函數，可以由下式計算得到：

在RHT模型里，引入屈服面主要考慮材料應變的硬化效應，將最大強度面縮放一定比例來得到：

式中：比例系數Fe為屈服強度和最大強度的比值，是靜水壓力的線性函數；Fcap（p）為彈性極限的帽子函數，用于限制靜水壓力作用下的彈性偏應力。

為了描述壓碎后混凝土材料的剩余強度，RHT模型單獨引入了殘余強度面，為壓力的指數函數，其定義為：

式中：Yr為殘余強度，B、M為殘余強度面參數，缺省值分別為1.6 和 0.61。

圖6為相應的三個強度面示意圖，對于在屈服強度面和最大強度面、最大強度面和殘余強度面之間的加載面，可以通過相互間的線性插值得到。

圖3 有限元數值模型Fig.3 Cross-section view of numerical model

圖4 沖擊波圖Fig.4 Pressure contours

圖5 不同應變率下鋼筋的應力應變關系Fig.5 Constitutive model of rebar at different strain rate

圖6 最大強度、屈服強度和殘余強度面Fig.6 Maximum strength，yield strength and residual strength surfaces

2.3.2 損傷定義及侵蝕算法

RHT模型認為損傷是塑性應變的累積，損傷變量D定義為：

式中：Δεp為等效塑性應變增量，Δμp為塑性體積應變增量，D1、D2為損傷參數。損傷度D隨塑性應變的累積而增長。材料完全失效后，損傷變量D=1.0，材料只能承受有限壓力，不能承受任何拉應力。

混凝土材料的應變率敏感性是引起動荷載作用下混凝土材料力學性能顯著區別于其準靜態荷載作用下的主要因素?；炷猎跍熟o態荷載作用下的峰值拉應變為0.000 2（為峰值壓應變的1/10），考慮到軟化階段，混凝土完全失效時的應變可認為是0.000 6。而在爆炸荷載作用下，混凝土的應變率為10 s－1～100 s－1，相應的強度動力放大因子可達到5.0甚至更大，但是斷裂應變的增強比強度的放大要小，故混凝土的失效拉應變取為0.001。

2.3.3 采用修正后參數對應力軟化特性的影響

文獻［13］對幾種廣泛使用的混凝土材料模型進行了全面的分析評價，并指出應用于AUTODYN軟件中的RHT模型缺省參數存在較大的問題。為了得到材料軟化的實際特性，對參數B、M、D1、D2和最小失效拉應變εminf進行修正。對于殘余強度面的可調參數B和M建議取值0.7和0.8，εminf取0.000 8代替原來的缺省值 0.01，系數D1推薦取值 0.015，D2=1.0 保持不變。混凝土的密度ρ=2.55×103kg/m3，抗壓強度為35.0 MPa，抗拉強度取為 3.5 MPa，剪切模量G為1.67 ×104MPa。

為了驗證RHT模型的實際特性，進行了數值模擬試驗用以評價該材料模型在不同荷載作用下的應力－應變關系。數值試驗采用單元體法，荷載由具有一定速率的位移控制，作用于節點上，施加必要的約束以限制剛體運動，如圖7所示。數值試驗顯示使用缺省參數時，RHT模型表現出在拉壓作用下應變軟化范圍出現不合實際的遞減過程，如圖8（a）所示的單軸拉伸試驗，完全破壞狀態下應變的數量級為0.01，遠遠大于一般物理實驗實際觀察的拉伸破壞應變，同樣的單軸受壓問題可見圖8（b）。采用修正后的參數，混凝土材料能得到比缺省參數較為合理的力學性狀。

圖7 RHT模型的單軸試驗模型Fig.7 Typical model setup of a uniaxial compression test

3 數值模擬

3.1 爆炸荷載

炸藥在空氣中爆炸，是一個能量高速釋放的過程，高溫高壓的爆炸產物高速膨脹，強烈壓縮周圍空氣，形成壓力很高的空氣沖擊波。當爆炸面產生的壓應力波傳播到介質與空氣的分界面時，沖擊波將發生反射，反射壓力通常為入射超壓的2～8倍。1969年，美國軍方出版的TM5－1300［14］結構抗爆作用規范給出的大范圍內爆炸荷載參數已被先前諸多應用所驗證，如圖9所示，為各比例距離峰值超壓的變化情況，其中Pso為TNT炸藥在自由空氣中爆炸時的沖擊波峰值超壓，Pr為反射沖擊波的峰值超壓。計算所得的鋼筋混凝土板面中心的壓力時程曲線如圖10所示，爆炸壓力在開始階段迅速達到峰值，并隨時間推移在1 ms內超壓降低為0。反射波的峰值超壓比TM5－1300值在自由空氣中測得的要大，但與TM5－1300反射超壓相比則要小，主要是因為數值模擬中鋼筋混凝土板的剛度與實際測量中的剛性體有一定偏差。

3.2 鋼筋混凝土板動力響應與破壞模式分析

在爆炸作用下板跨中最大撓度是衡量其破壞程度的一個重要指標。由圖11和圖12可知，最大速度出現在爆炸沖擊波結束時，而后結構由于慣性作用將繼續向下運動，另外，隨著炸藥量的增加，板跨中的撓度隨之增大。當炸藥量為1 kg時，板在20 ms時停止運動，最終位移為17 mm。當炸藥量為2 kg時，板的豎向位移為61 mm。5 kg TNT炸藥作用下，板跨中在0.77 ms達到最大速度 16.7 m/s，隨后速度迅速衰減，在20 ms時位移為143 mm。

圖8 默認參數與新參數的RHT模型應力應變曲線比較Fig.8 Comparison of stress-strain curves generated by RHT model with default and newly proposed model parameter values

圖9 TNT球形裝藥海平面處爆炸荷載參數Fig.9 Pressure and reflected pressure for a spherical TNT detonation at sea level

圖10 鋼筋混凝土板面處的反射壓力時程曲線Fig.10 Reflected pressure on the surface of the RC slab

圖11 不同炸藥量下板中心點豎向速度時程曲線Fig.11 Vertical velocity of middle point in slab under different blast loadings

圖12 不同炸藥量下板中心點豎向位移時程曲線Fig.12 Vertical displacement of middle point in slab under different blast loadings

鋼筋混凝土板在不同炸藥量作用下的破壞模式如圖13所示。由圖可知，在不同的爆炸荷載作用下，壓縮應力波均在板的迎爆面產生損傷，并傳播至板的背爆面形成強拉伸波，造成背爆面混凝土的層裂和崩塌。隨著炸藥量的增大，鋼筋混凝土板的破壞越嚴重，并表現出不同的破壞模式。在1 kg TNT炸藥作用時，板正面縱向中部有一條寬度為100 mm左右的混凝土壓碎帶，背面縱向中部有一條寬度約5 mm的受拉破壞裂縫，板主要發生彎曲破壞；當炸藥量為2 kg時，混凝土板受拉區破壞較為嚴重，板背面縱向中部出現多條受拉破壞的裂縫，表現出明顯的雙向彎曲。在5 kg TNT炸藥作用時，鋼筋混凝土板破壞嚴重，板整體表現為彎曲破壞，中部橫向裂縫增多，混凝土板被分割成許多塊體，部分鋼筋屈服但沒發生斷裂；當炸藥量為10 kg時，板的破壞主要集中在中央，呈現彎曲沖切耦合破壞，板中央部分縱向受拉鋼筋發生斷裂，在迎爆面與背爆面之間的板中心區域發生上下貫通破壞，形成了一個直徑大于1.0 m的爆坑，有混凝土碎塊飛出。

Fig.13 不同炸藥量作用下鋼筋混凝土板的損傷Fig.13 Damage distribution of RC slab under different explosive charge

圖14 2kg Comp B炸藥作用下鋼筋混凝土板的損傷Fig.14 Damage distribution of RC slab in 2 kg Comp B charge case

3.3 數值模擬與試驗結果的比較

選用文獻［15］鋼筋混凝土板試件進行建模，板的凈尺寸為1.3 m ×1.0 m ×0.12 m，鋼筋直徑10 mm，縱向受力鋼筋間距為100 mm，分布鋼筋間距為120 mm，雙層雙向布置。2 kg Comp B炸藥置于距板頂面中心0.6 m高度處。在爆炸荷載作用下，壓縮應力波在板的迎爆面產生損傷，并傳播至板的背爆面形成強拉伸波，造成背爆面混凝土的剝落和層裂，破壞過程如圖14所示。鋼筋混凝土板背面縱向中部較橫向中部先出現受拉破壞裂縫，在2 ms時，橫向裂縫發展，最終寬度約為5 mm，板發生明顯的雙向彎曲破壞。在板中心區域，出現環形裂縫，具有雙向彎曲的角部裂縫特征。鋼筋混凝土板受拉區破壞嚴重，板背面中部出現多條橫向受拉破壞的裂縫，板正面和背面的破壞特征與試驗得到的圖15和圖16［15］吻合較好。

圖15 正面混凝土壓縮裂縫Fig.15 Compression crack of RC slab upper surface

圖16 背面的受拉破壞Fig.16 Tensile damage of RC slab back surface

4 結論

本文采用AUTODYN軟件建立了炸藥、空氣和鋼筋混凝土板的三維數值模型，考慮了應變率對鋼筋和混凝土材料動力本構特性的影響以及炸藥－空氣－結構之間的流固耦合相互作用，分析了不同炸藥量作用下鋼筋混凝土板的損傷機理和破壞特征，結果表明：

（1）鋼筋混凝土板面承受的爆炸荷載與自由場沖擊波大小和板的剛度密切相關。入射超壓越大，則反射越強，升壓時間短，衰減迅速。

（2）建立的流固耦合數值模擬方法可以合理地展現鋼筋混凝土板從混凝土開裂、碎片形成、部分鋼筋屈服斷裂到板局部震塌的動態演變過程。

（3）隨著炸藥量的加大，鋼筋混凝土板的破壞逐漸由整體彎曲破壞轉變為局部的沖切破壞。

［1］Gong S F，Lu Y，Jin W L.Simulation of air blast load and its effect on RC structures［J］.Transactions of Tianjin University，2006，12（suppl）：165－170.

［2］Yuan L，Gong S F，Jin W L.Spallation mechanism of RC slabs under contact detonation［J］.Transaction of Tianjin University，2008，14（6）：464－469.

［3］龔順風，金偉良，何 勇.內部爆炸荷載作用下鋼筋混凝土板的動力響應研究［J］.振動工程學報，2008，21（5）：516－520.

［4］閻 石，張 亮，王 丹.鋼筋混凝土板在爆炸荷載作用下的破壞模式分析［J］.沈陽建筑大學學報，2005，21（3）：177－180.

［5］張 舵，盧芳云，王瑞峰.鋼筋混凝土板在爆炸作用下的破壞研究［J］.彈道學報，2008，20（2）：13－16.

［6］ Gong S F，Lu Y，Tu Z G，et al.Validation study on numerical simulation of RC response to close-in blast with a fully coupled model［J］. StructuralEngineering and Mechanics，2009，32（2）：283 －300.

［7］ Zhou X Q，Kuznetsov V A，Hao H，et al.Numerical prediction of concrete slab response to blast loading［J］.International Journal of Impact Engineering，2008，35（10）：1186－1200.

［8］Xu K，Lu Y.Numerical simulation study of spallation in reinforced concrete plates subjected to blast loading［J］.Computers＆ Structures，2005，84（6）：431－438.

［9］ AUTODYN theory manual，revision 6.1 ［M］.California：Century Dynamics，San Ramon，2006.

［10］ AUTODYN User Manual，Revision 6.1 ［M］.California：Century Dynamics，San Ramon，2006.

［11］Johnson G R，Cook W H.A constitutive model and data for metals subjected to large strains，high strain rates and high temperatures［C］.Proceedings of the seventh international symposium on Ballistics.Hague：1983，541－547.

［12］ Riedel W，Thoma K，Hiermaier S，et al.Penetration of reinforced concrete by BETA－B－500 numerical analysis using a new macroscopic concrete model for hydrocodes［C］.Proceeding of the Ninth International Symposium on Interaction of the Effects of Munitions with Structures，1999，315－322.

［13］Tu Z G，Lu Y.Evaluation of typical concrete material models used in hydrocodes for high dynamic response simulations［J］.International Journal of Impact Engineering，2009，36（1）：132－146.

［14］ TM5－1300.Structures to resist the effects of accidental Explosions［S］.US Army，USA，1990.

［15］孫文彬.鋼筋混凝土對邊簡支板的爆炸試驗［J］.力學與實踐，2008，30（4）：58－60.