焦晓龙，赵鹏铎，姚养无，张 磊，李旭东，池 海

（1. 海军研究院，北京 100161；2. 中北大学机电工程学院，山西 太原 030051）

舰艇作为掌握海洋守卫主动权的利器，只有具备较强生命力时，才能有机会发挥舰载武器的威力，才能说具备出色的作战能力。随着半穿甲武器的迅速发展，学者们对炸药在舰艇空间有限的封闭舱室内爆炸的相关问题进行了研究。侯海量等[1]进行了典型单舱室结构内爆模型实验，分析了内爆载荷和舱室板架结构失效模式，指出角隅部位的汇聚冲击波强度远大于壁面反射冲击波、主要失效模式为沿角隅撕裂并发生大挠度外翻。另外，樊壮卿等[2]利用流固耦合算法，对该实验进行了进一步研究，不仅验证了毁伤过程的一致性，还得到了单舱室结构的失效特征和主要破坏载荷。随着对内爆下单舱室结构毁伤效应研究的不断深入，学者们对内爆下多舱室结构的毁伤情况也进行了探索。孔祥韶等[3-4]、严波等[5]对内爆炸环境下舷侧防护结构的响应进行了研究：孔祥韶等[3-4]通过数值模拟和实验研究，分析了高速破片和冲击波耦合作用下舷侧防护结构的冲击响应；严波等[5]利用数值模拟，分析了隔板间距疏密对破坏形式的影响。李营等[6]开展了3 个并连舱室结构在舱内爆炸作用下的毁伤特性实验，测量了爆炸破片形态尺寸、破片速度和冲击波载荷等，分析了塑性变形、毁伤模式等结构毁伤特点。姚术健[7]基于固支方板的运动方程，先通过量纲分析法，给出了适用于爆炸舱（爆炸发生所在舱室）的无量纲毁伤数，再根据实验结果和数值模拟，指出单层多箱室结构主要有十字形和非十字形两种破坏模式，最后给出了破坏模式快速预测方法，该方法考虑了结构的尺度效应和应变率的影响且适用于工程应用。

目前，对内爆下多舱室结构毁伤效应的研究多集中于爆炸舱，而涉及邻舱（邻近舱室）的还很少。实际上，半穿甲反舰武器能对多个舱室甚至舱段造成破坏，内爆载荷对邻舱的设备和人员的威胁非常大[8]。所以，需要对内爆载荷作用下包括邻舱在内的整个多舱室结构毁伤效应进行研究，本文中将探讨在不同药量TNT 内爆下大尺寸多舱室结构的毁伤效应。

1 数值计算模型和方法

采用AUTODYN 软件的多物质流固耦合算法，能够模拟计算舱内爆炸过程中载荷与结构间的相互作用，可以准确模拟爆轰波的产生、反射以及汇聚等载荷特性，并能准确反映结构的响应。

1.1 有限元模型

多舱室结构模型由边长2 m 的立方体型单舱室按3×3×3 堆叠而成，图1 为完整模型和1/8 模型。

图1 多舱室结构模型和1/8 模型Fig.1 Whole model and 1/8 model of multi-cabin structure

有限元模型包括舱室结构有限元模型、空气域有限元模型和炸药有限元模型3 种。模型舱壁和甲的厚度均为6 mm，远小于舱壁结构尺寸，因此舱室结构模型采用4 节点壳单元，空气和炸药有限元模型均采用8 节点体单元。考虑结构会发生大变形，为保证变形后的结构依然能受到爆炸载荷作用，多舱室模型完全浸没在空气域中。结构具有对称性，为简化计算和节省时间，采用1/8 结构模型，如图2（a）所示，整体有限元模型共有约196 万个单元。

图2 多舱室结构有限元模型Fig.2 Finite element model of multi-cabin structure

1.2 模型材料和参数

与舱内爆炸载荷相比，甲板和舱壁为薄板，在舱内爆炸载荷作用下会发生塑性大变形，甚至产生破口、连接部位和开口部位的撕裂。考虑舱壁的失效，而所使用钢材的屈服极限较高，本文中设置失效应变为0.3[9]。

空气采用理想气体状态方程[10]描述，该方程可以准确反映实际气体的物理行为。状态方程为：

式中：p 为压强，γ 为绝热指数，ρ 为密度，e 为内能。状态方程参数分别为：γ=1.4，ρ=1.225 kg/m3, e=206.8 kJ/kg。

TNT 炸药采用Jones-Wilkins-Lee（JWL）状态方程描述，该方程可以精确计算炸药爆轰驱动金属板壳结构的加速过程，被广泛应用于工程计算中[11]。状态方程为：

式中：A、B、R1、R2、ω 为实验拟合参数，E 为爆轰产物的体积内能，v=ρ0/ρt为爆轰产物的相对比容，ρ0为炸药初始密度，ρt为爆轰后密度。

2 数值模拟方法的验证

2.1 实验验证

采用实验数据和经验公式相结合的方法，对本文的数值模拟方法进行验证。根据连赟猛[12]的实验，建立有限元模型（见图2（b））进行验证。其中，实验装置为密闭长方体形钢筋混凝土结构，如图3（a）所示。装置壁厚150 mm（由实验中测得的最大载荷强度1.18 MPa），可忽略实验装置壁面的动态变形，即将壁面视为刚体，模拟计算时可用相同厚度的钢材代替。测点分布[12]如图3（b）所示。

图3 实验装置和测点布置[12]Fig.3 Experimental device and pressure gauge arrangement[12]

文献[12]中给出了壁压数据。其中，测点P1、P2 正对爆源（见图3（b）），这两个测点的压力曲线的首峰为第1 次正反射冲击波峰，用正反射超压的计算公式[13]求解验证：

式中：Δpr为反射压力，Δpf为入射压力，p0为环境压力，γ 为绝热指数。绝热指数γ 的取值随入射超压强度有所改变，依据文献[13]，本文中取γ=1.4。

入射压力为爆炸产生的初始冲击波超压，为：

式中：Δpf为冲击波峰值超压，MPa；Z=R/W1/3为相对距离，R 为爆距，m，W 为爆炸当量，kg。

该经验公式以大量实验为基础，可信度高，为我国国防工程设计规范中的空爆冲击波超压计算公式[14]。

2.2 验证结果

数值模拟的各测点压力峰值见表1，并利用经验公式加以佐证，提高本次研究的置信度。其中，相对于实测结果，数值模拟计算结果的平均误差为5.37%，最大误差为8.92%，计算超压值与实测值接近。可以认为，本文中采用的数值模拟方法能够较好地模拟舱内爆炸时的壁面压力特征，可靠性较高。

表1 实测压力峰值和数值模拟计算结果的对比Table 1 Comparison of measured peak pressures and numerical simulation results

3 结构毁伤等级划分

多舱室结构的毁伤情况复杂，为方便表述，进行以下分类：（1）依据相对爆炸舱的位置，将邻舱分为共面邻舱、共边界邻舱和共点邻舱；（2）依据变形失效类型，将邻舱舱壁及边界进行分类编号，如图4 所示，具体描述见表2。

图4 结构分类Fig.4 Structural classification

表2 结构分类Table 2 Structural classification

分析计算结果，发现多舱室结构的毁伤集中体现在各个舱壁的变形失效上，主要包括舱壁挠曲大变形、舱壁中心冲切破坏、舱壁沿边界发生翻转撕裂、舱壁在边界处剪切失效等，这些变形破坏随着爆炸载荷的增强相继发生。为了体现内爆载荷对多舱室结构的毁伤特点，结合以往经验划分毁伤等级，见表3。各种毁伤等级对应的典型毁伤情况，如图5 所示。

表3 毁伤等级的描述Table 3 Description of damage grade

图5 各种毁伤等级对应的典型毁伤情况Fig.5 Typical damage situations corresponding to various damage grades

上述毁伤情况符合舱内爆炸载荷作用下典型舱室结构的毁伤模式。为了支撑上述结论，对主要内爆冲击载荷之中的准静态压力采用经验公式验证。李德聪等[15]对一些常用的计算公式进行了总结，并认为当m/V＜0.5 时，劳氏规范公式的计算结果与其他公式相差较大。结合本文情况，选用Carlson 公式进行验证：

式中：pqs为准静态压力，MPa；m 为TNT 药量，kg；V 为结构容积，m3。

这里，验证了准静态压力作用明显的几个TNT 药量状况，见表4。由于准静态压力在某个范围内浮动，采用确定的值，结果难免有误差。经对比，计算结果平均相对误差5.995%，最大误差9.43%，结果在可接受的精度范围内。表4中，计算结果普遍低于经验公式的，可能是因计算中未考虑炸药后燃烧效应，这需进一步研究。可以认为，本文中采用的数值模拟计算方法能够较好模拟舱内爆炸，具有进一步研究的可行性。

表4 数值模拟结果与准静态压力经验公式结果的对比Table 4 Comparison between numerical simulation results and quasi-static pressure empirical formula results

4 量纲分析和快速毁伤预测

4.1 量纲分析

式（8）中的5 个变量均无量纲，但各变量间关系不确定。为了进一步探索结构变形规律，采用指数乘积的形式表示：

用式（9）对计算结果进行分析，发现δ/H 和l/H 均与m/V 存在较好的线性关系，如图6 所示，可令式（9）中α=-1、γ=1。由文献[7]，δ/H 与其分析得到的无量纲量存在明显线性关系，可认为在某个确定结构中，炸药释放总能量Q 和H 的比值Q/H 与δ/H 线性相关。对同种炸药，单位质量炸药释放的能量为定值，即m 与Q 线性相关，因而m/H 与δ/H 线性相关。δ 和l 是独立的，则m/H 与l/H 也线性相关，则可令β=1。定义内爆下多舱室结构毁伤效应无量纲毁伤数为Bc，则：

图6 最大挠度和裂缝长度与TNT 药量的关系Fig.6 Relationship of maximum deflection and tearing length with TNT charge

4.2 快速毁伤预测方法及验证

经上述分析，采用舱壁最大挠度δ 和固定边界裂缝长度l 评价多舱室结构的毁伤情况，确定毁伤等级。拟合得到毁伤等级预测曲线（见图7）和对应公式：

由于5 级毁伤表现为三面交接处反向撕裂并伴随着边界扭曲，不适于用δ 或l 表征，因此图7 中未显示。

图7 内爆下多舱室结构毁伤等级预测曲线Fig.7 Prediction curve of damage grade of multi-cabin structure under internal blast

图7 中各线段不连续的原因包括系统误差和偶然误差，主要有：（1）各类舱壁所受载荷类型有差异；（2）各类舱壁的受载方向不同，在图7 中体现为各线段斜率不同；（3）有某药量爆炸时，舱壁处于弹性阶段，未发生塑性变形，导致毁伤等级偏小；（4）表征量无法准确取得。

利用文献[23]中的实验简述快速毁伤预测方法，并作检验。如图8 所示，该文献的装置由相同结构的小舱室组成，每个小舱室结构尺寸为0.15 m×0.125 m×0.1 m，纵向板厚3 mm，结构材料屈服强度σ 为464 MPa，TNT 药量为0.05 kg。

图8 实验装置[23]Fig.8 Experimental device[23]

根据实验情况，利用式（10）求得Bc为83.9，代入式（11）第3 个式子，得到毁伤等级为3.06，对应的破损现象为爆炸舱舱壁沿边界被撕裂失效后飞出(见表3)，且飞出舱壁与后续内爆载荷耦合作用于A1 类舱壁，使其发生轻微破损。文献[23]中实验后装置内部破损情况如图9 所示，有3 个爆炸舱舱壁破损严重，且每个均有3 条边界完全撕裂并发生大挠度外翻，装置的其他结构变形不明显。

图9 装置内部失效状况[23]Fig.9 Device internal failure[23]

与实验结果相比，快速预测结果偏大，原因有：（1）快速预测方法基于数值模拟计算得出，其单舱室尺寸为2 m×2 m×2 m，而该实验的单舱室尺寸为0.15 m×0.125 m×0.1 m；（2）数值模拟中，舱壁间以理想的连接方式处理，未考虑焊接的影响；（3）数值模拟中未考虑药孔的影响，实验装置有多处药孔（见图9）；（4）数值模拟中爆点位置在结构中心，而实验中存在偏差。对以上原因，可以通过大量实验结果来修正、完善。总之，采用本文预测方法，可对内爆下多舱室结构的毁伤情况进行快速预测。

5 结 论

运用数值模拟，研究了内爆载荷作用下多舱室结构的毁伤特性。运用ANSYS/AUTODYN 非线性有限元软件中的多物质流固耦合算法，进行了大量不同TNT 药量的数值模拟计算。并依据文献[12]中的实验数据、置信度较高的经验公式，验证了数值模拟计算方法的可行性。然后，对计算结果运用量纲分析，推导了关于球形TNT 结构中心内爆下典型多舱室结构的无量纲毁伤数，并给出了基于本文模拟计算结果的毁伤等级预测公式和快速毁伤预测方法。主要结论如下。

（1）内爆载荷下多舱室结构的毁伤特征可用基础结构的不同失效形式描述，如舱壁主要包括舱壁中心冲切失效、舱壁撕裂失效、舱壁在边界处剪切失效及三面连接处反向撕裂失效。

（2）内爆载荷下舱壁挠曲变形产生的最大挠度和撕裂失效产生的裂缝长度均与m/V 存在明显线性关系；内爆载荷下多舱室结构毁伤程度与m/V 密切相关，并可用关于m/V 的方程表示。

（3）基于基础结构失效模式，结合无量纲分析的快速预测方法，能够实现内爆下多舱室结构的快速毁伤预测，可用于内爆下多舱室结构的毁伤特性研究，为舰船毁伤研究提供参考。

（4）还发现内爆载荷泄出后，泄口处会产生负压区，可能使破片飞回，对结构造成二次破坏，对此还需进一步研究。