水下接觸爆炸下鋼筋混凝土板毀傷模式探究

胡錦華

(中鐵十六局集團北京軌道交通工程建設有限公司，北京 101100)

國內外學者圍繞爆炸沖擊下混凝土板的動態響應開展了廣泛的研究，其研究手段主要為現場試驗和數值仿真分析。Wu等[1]通過一系列素混凝土板和不同纖維增強混凝土板的空中抗爆研究，并以普通鋼筋混凝土板作為對照，探究了纖維對混凝土抗爆性能的影響；Zhou等[2]通過空中接觸爆炸實驗，研究了混凝土板的爆炸毀傷全過程，并提出了混凝土板爆炸毀傷塊體尺寸的預測公式；Lavaei等[3]以混凝土板為研究對象，開展了空中爆炸試驗，并分析了爆心距對混凝土板毀傷模式的影響；Li等[4-6]開展了一系列普通鋼筋混凝土板和高強鋼筋混凝土板的空中接觸爆炸試驗，探究了普通混凝土板、纖維增強混凝土板、金屬纖維增強混凝土板的抗爆性能，并提出了對應的預測公式；趙小華等[7]通過拉格朗日-歐拉模型，探究了空中近場爆炸作用下鋼筋混凝土板的毀傷破壞特性；Yun等[8]采用數值分析方法，研究了空中近場爆炸沖擊下鋼板加固混凝土板的毀傷破壞特性；Wang等[9]通過建立的耦合模型，分析了空中近場爆炸荷載作用下鋼筋混凝土板的破壞模式和損傷分布，并通過與現場試驗的對比驗證了仿真計算的可靠性；王高輝等[10]采用數值仿真方法，研究了自由邊界和結構面對水下爆炸沖擊波傳播特性的影響；Yang等[11]通過數值仿真分析，考慮箱涵內水位的影響，分析了地表爆炸對淺埋輸水箱涵的毀傷特性。以上研究主要是圍繞空中爆炸開展的，而相對于空中爆炸，同等藥量下水下爆炸往往能對結構造成更為嚴重的破壞。

圍繞水下爆炸沖擊波的傳播特性和對混凝土結構的毀傷破壞特性，婁浩然等[12]采用離心模型實驗研究了水下爆炸沖擊波傳播特性及氣泡脈動規律，并重現了超重力場下的水下爆炸過程；韓崇剛等[13]基于AUTODYN軟件建立了水下爆炸一維楔形模型，探究了炸藥沖擊波壓力以及氣泡脈動隨時間的變化特性；李惠明等[14]探究了圓柱形裝藥和半球形空穴裝藥水中爆炸過程，對比分析了沖擊波壓力、氣泡直徑等參數；閆秋實等[15]開展了水下近場爆炸沖擊下鋼筋混凝土樁的動力響應和破壞模式研究，并提出了該型鋼筋混凝土樁特定工況下的抗爆安全范圍；孔祥清等[16]通過有限元分析，對比分析了空中和水下近場爆炸作用下爆炸沖擊波的傳播規律，以及鋼筋混凝土板的動力毀傷破壞特性。然而，上述研究主要是圍繞水下非接觸爆炸開展的，而針對水下接觸爆炸作用下鋼筋混凝土板的毀傷破壞模式研究較少。

本文通過現場爆炸試驗以及建立的鋼筋、混凝土、水體、炸藥全耦合模型，對水下接觸爆炸下沖擊波的傳播特性和鋼筋混凝土板的毀傷破壞模式進行了研究，研究結果可為鋼筋混凝土結構的水下防爆設計和研究提供參考。

1 現場試驗

為了探究水下接觸爆炸下鋼筋混凝土板的毀傷模式和破壞特性，并為動態響應過程和毀傷機制的數值分析提供試驗基礎，特進行了水下接觸爆炸現場試驗，試驗所用試件為方形鋼筋混凝土板。

1.1 試驗設置

試驗中所用炸藥類型為巖石乳化炸藥，炸藥質量為50 g,密度1.05 g/cm3，猛度12 mm，Chapman-Jouguet爆速為3 850 m/s，Chapman-Jouguet壓力為3.70 GPa，炸藥放置于鋼筋混凝土板上表面正中心處，即接觸爆炸。炸藥采用電雷管引爆，試驗中電雷管從炸藥頂部垂直插入。

本次試驗所使用的鋼筋混凝土試件尺寸為500 mm×500 mm×60 mm，其中鋼筋混凝土板的配筋為6 mm圓鋼筋，鋼筋間距為100 mm，保護層厚度為20 mm(見圖1)。鋼筋的屈服強度和楊氏模量分別為600 MPa和200 GPa，采用無側限抗壓強度為28.2 MPa的混凝土澆筑。

圖1 鋼筋混凝土板幾何尺寸Fig.1 Geometrical dimensions of reinforced concrete slab

試驗中，采用鋼架對鋼筋混凝土板兩端進行固定，該支架是由4段槽鋼焊接而成，上層卡槽用于試件的安裝固定。安裝時需先將試件卡入槽中，然后在板固定邊的上下面卡入三角鋼，最后將制備好的楔形鋼片卡入三角鋼和槽鋼的間隙中，以提供近似固定的邊界條件。

試驗前將安裝試件后的支撐結構用吊機放入水中，其中板頂面距離水面0.5 m，鋼架的下部為原始湖底，混凝土板的底部距離湖底約為0.2 m，但由于水下接觸爆炸的大部分能量用于毀傷試件，透過板底面傳播至水底的沖擊波較少，同時相對粗糙和疏松的水底將吸收大多數沖擊波能量，因此，可以不用考慮水底反射沖擊波對混凝土板毀傷的影響。

1.2 試驗結果

為更好地觀察鋼筋混凝土試件的破壞模式，給出了鋼筋混凝土板迎爆面和底面的現場試驗結果(見圖2)。

從圖2a中可以看出，水下接觸爆炸沖擊下，鋼筋混凝土迎爆面遭遇了嚴重的沖切毀傷，中部位置的混凝土破碎、掉落，缺失面積占總面積的78%。迎爆面四周殘留的混凝土塊體內充滿裂隙，在內部配筋的支撐下未剝離掉落，內部配筋雖然未出現斷裂破壞，但整體也出現了較大形變。由鋼筋混凝土板底面的破壞結果(見圖2b)可以看出，鋼筋混凝土試件的底面在水下接觸爆炸荷載沖擊下已完全破碎、剝落，可直接觀察到內部的配筋，僅被鋼架固定的兩側有所殘余，鋼筋混凝土板已完全失去強度和剛度。

2 數值仿真分析

2.1 材料本構模型及參數

陶瓷、混凝土和玻璃等材料在剛性炮彈侵徹沖擊和爆炸沖擊等工況中會展現出應變率效應。RHT[17]模型通過引入壓力依賴性、應變速率敏感性、壓縮損傷軟化，并引入了偏應力張量第三不變量對破壞面形狀的影響，能夠較好地描述混凝土材料在侵徹和爆炸沖擊等高加載率荷載作用下的破壞過程。該模型的3個失效面如圖3所示。

圖3 RHT本構模型的3個失效面Fig.3 Three failure surfaces of the RHT constitutive model

其失效面為

(1)

(2)

其損傷定義公式為

(3)

(4)

模型計算中炸藥的爆轟和沖擊波傳播過程通過JWL狀態方程來描述，該模型中壓力、內能和相對體積間可通過公式描述：

(5)

式中：p為爆轟壓力；V為炸藥起爆產物體積與其初始體積的比值；A1、B1、R1、R2、ω為特征參數；E0為單位體積內能；e為爆轟產物的內能。E0=6.0 GJ/m3，密度ρ=1 630 kg/m3，爆速D=6 930 m/s，Pcj=21 GPa，V=1.0，A1=373.77 GPa，B1=3.75 GPa，R1=4.15，R2=0.90，ω=0.35[18]。

通過Polynomial方程描述水體在壓縮和膨脹時的狀態，其公式如下：

壓縮時，即μ>0，其方程為

p=A1μ+A2μ2+A3μ3+(B0+B1μ)ρ0e

(6)

膨脹時，即μ<0，其方程為

p=T1μ+T2μ2+B0ρ0e

(7)

式中：p為水中壓力；μ為壓縮比，μ=ρ/ρ0-1;e為水的內能；ρ0為水密度取為1 g/cm3；A1=T1=2.2×e6kPa，A2=9.54×e6kPa，A3=1.46×e6kPa，B0=B1=0.28，T2=0[18]。

2.2 耦合模型的建立

采用歐拉-拉格朗日耦合數值分析方法對鋼筋混凝土板空中接觸爆炸進行數值建模(見圖4)?？梢钥闯?，在仿真模型中，混凝土材料是使用拉格朗日網格進行描述的，而水體和炸藥則采用歐拉網格進行模擬，且水體四周的切斷邊界被設定為無反射邊界以減小邊界對計算結果的影響。鋼筋是通過Beam單元建模，并通過共節點與混凝土材料耦合，同時假定鋼筋與混凝土之間為理想固結，不發生滑動。這種假定的可行性和準確性已經在文獻中被證實[19]。

圖4 鋼筋混凝土板遭受空中接觸爆炸沖擊的幾何模型Fig.4 Geometric model of reinforced concrete slab subjected to air contact explosion impact

混凝土板的兩端面上下側施加固定約束，與試驗中的邊界條件保持一致。在仿真模型中，炸藥和鋼筋單元的尺寸為5 mm，水體的單元尺寸為10 mm。為保證混凝土與鋼筋的共節點，混凝土板長寬方向的單元尺寸為5 mm，厚度方向的單元尺寸為3 mm，混凝土材料共有200 000個單元。

2.3 動態響應及損傷發展過程

水下接觸爆炸沖擊作用下鋼筋混凝土板的毀傷過程如圖5所示。炸藥在鋼筋混凝土上表面上起爆后，爆轟產物和高壓沖擊波直接造成爆心區域的壓縮毀傷。隨著起爆時間的發展，上表面的毀傷區域面積和深度均不斷增大，最終形成貫穿毀傷區域。此外，爆心下方的混凝土板內配筋出現較大形變，大面積配筋進入塑性區。

圖5 水下接觸爆炸沖擊下鋼筋混凝土板的毀傷過程Fig.5 Damage process of reinforced concrete slab under the impact of underwater contact explosion

從圖5可以看出，炸藥在水中起爆后，爆轟產物以及爆炸產生的水沖擊波將直接作用于板的上表面，與板的迎爆面相互作用，導致板的迎爆面出現圓形毀傷區域(見圖5a)；炸藥起爆0.04 ms后，板的底面中心區域也出現毀傷區域(見圖5b)，這是因為水下爆炸沖擊波傳播至板底面后，因混凝土的波阻抗大于水介質，部分爆炸沖擊波在板底面發生反射，在入射沖擊波和反射沖擊波的聯合作用下，混凝土板底面出現沖切毀傷；隨后爆心處的配筋開始發生形變，板的側面也出現損傷區域(見圖5c)；當t=0.09 ms時，鋼筋混凝土板迎爆面和底面的損傷區域連通，板中部形成了貫穿破壞(見圖5d)；當t=0.6 ms后，鋼筋混凝土板的毀傷區域不再變化，鋼筋混凝土板整體遭遇嚴重的破壞，失去剛度和強度。以上研究表明，水下接觸爆炸將對鋼筋混凝土板造成嚴重的局部和整體毀傷破壞，鋼筋混凝土板水下接觸爆炸數值模擬結果如圖6所示。

圖6 鋼筋混凝土板水下接觸爆炸數值仿真結果Fig.6 Numerical simulation results of underwater contact explosion of reinforced concrete slab

為了分析水下接觸爆炸下鋼筋混凝土板內壓縮波的傳播特性，在鋼筋混凝土板的中部設置了1#、2#、3#共3個測點，測點壓力時程如圖7所示。

圖7 水下接觸爆炸下鋼筋混凝土板內壓力時程Fig.7 Time-history of internal pressure of reinforced concrete slab under underwater contact explosion

從圖7中可以看出，位于爆心正下方的測點1的壓力峰值最大，達到68.49 MPa，而后迅速衰減，并出現多個壓力峰值，這是炸藥爆轟產物及沖擊波在板的迎爆面與底面多次相互作用的結果；而當測點遠離爆心后，其壓力峰值迅速衰減，測點2處的壓力峰值已經降低到21.75 MPa；位于試件邊緣的測點3處的壓力峰值僅為4.38 MPa。此外，相對于測點1，測點2和測點3處的壓力峰值衰減更快，作用時間也更短。

水下接觸爆炸下，位于爆心正下方的混凝土板迎爆面和底面的壓力峰值分別為96.57、53.81 MPa，均遠高于混凝土材料的動抗壓強度，故板迎爆面和底面的毀傷機制主要為壓縮沖切破壞。

3 結論

1)對比數值仿真結果與現場試驗結果可知，建立的數值耦合模型能夠很好地模擬水下接觸爆炸下鋼筋混凝土板的毀傷破壞過程，以及爆炸沖擊波在板內的傳播特性，是一種研究水下爆炸沖擊下鋼筋混凝土結構動態響應的有效工具。

2)通過現場試驗和數值仿真分析發現，水下接觸爆炸作用下，混凝土板整體出現嚴重的沖切損傷，爆心區域板完全破碎、脫落，混凝土板整體變形較大，混凝土板迎爆面和底面的毀傷機制主要為壓縮沖切破壞，因此在鋼筋混凝土板抗爆設計中，應格外關注水下接觸爆炸的安全防護。