程 帅，师莹菊，殷文骏，刘文祥，唐仕英，张德志

（西北核技术研究院强动载与效应实验室，陕西 西安 710024）

泡沫铝具有轻质、可加工性强、成本低等优点；且承受动态载荷时，能够在固定的应力下发生较大的塑性变形，被认为是一种理想的缓冲、吸能材料，已经广泛应用于各类抗爆结构的设计。任新见等[1]采用泡沫铝夹心结构，提高了圆柱形排爆管的性能并减少了质量。刘新让等[2]通过数值模拟发现在相同的质量和爆炸当量下，泡沫铝夹芯筒的整体变形小于实心钢筒。Goel 等[3]和Santosa 等[4]通过实验和数值模拟研究了泡沫铝提升结构抗爆性能的机理，并指出泡沫铝的厚度是影响抗爆性能的重要因素。张培文等[5]通过数值模拟分析了钢板厚度对泡沫铝夹芯结构抗爆性能的影响，发现泡沫铝能够有效降低爆炸载荷下结构的变形。

然而，总结成功应用泡沫铝提高结构抗爆能力的工作不难发现，这些案例中存在爆炸比距离偏大、爆炸载荷强度不高、泡沫材料厚度使被防护结构尺寸显著增大等问题。随着爆炸比距离减小，泡沫防护结构的变形、破坏模式将发生变化[6]。低强度冲击波与泡沫材料作用时，泡沫材料只在弹性和平台应力区域工作，可通过变形吸收冲击波能量；但高强度冲击波会引起泡沫材料的致密化，导致被防护结构发生更严重的破坏[7]。Skews 等[8]发现强冲击波作用下泡沫材料中形成的速度低于空气冲击波的压缩波可能是导致变形增大的重要原因。Li 等[9]使用一维弹簧-质量模型分析了泡沫材料导致冲击波强化的条件；Tan 等[10-11]使用一维冲击模型分析了泡沫材料对冲击载荷的增强作用，并指出惯性和局部压缩现象是导致载荷增强的主要原因。Lopatnikov 等[12]通过数值模拟和霍普金森杆实验研究了强冲击载荷下泡沫金属的响应。大量的研究结果表明，当冲击波强度足够、压缩波可以到达被保护结构时，波的反射就会引起载荷增强现象，导致结构发生更严重的破坏[13-14]。

综上，采用泡沫材料可以减少结构质量，提高结构的抗爆能力，但同时必须确定泡沫材料能够起防护作用的极限载荷强度，否则可能导致结构发生更严重的破坏。为在不改变钢筒内径、壁厚等特征尺寸的情况下，提高承受内部爆炸载荷钢筒的抗爆能力，本文中通过数值模拟分析泡沫铝内衬对钢筒变形的影响，总结泡沫铝内衬影响钢筒变形的主要机制，以期为钢筒抗爆能力提升设计提供参考。

1 实 验

如图1 所示，本文实验以内径为100 mm、壁厚为6 和12 mm 的钢筒为加固对象，在钢筒内部添加泡沫铝内衬，对比钢筒在内部爆炸载荷下的变形，研究泡沫铝内衬对抗爆能力的影响。根据秦学军等[15]的研究结果，承受内部爆炸载荷的细长柱形钢筒，当钢筒长径比超过6 时，可以忽略边界对爆心环面变形的影响。因此，实验中选取的钢筒长度为600 mm。钢筒材料为Q345R，是压力容器、抗爆结构设计的常用材料。

实验时用卡箍将钢筒固定在图2 所示的刚性平台上，再将泡沫铝内衬和炸药装入钢筒中心。其中泡沫铝内衬的外直径为100 mm，密度为0.48 g/cm3，孔隙率约为82%，其静态条件下测得的应力应变曲线如图3 所示。由于内部爆炸载荷作用下，钢筒的变形主要集中在爆心环面2 倍半径距离之内[16]，因此实验中选用泡沫铝的长度为200 mm。使用基于多普勒干涉原理的光子多普勒测速仪[17-18]（PDV，photonic Doppler velocimetry）测量钢筒爆心环面外表面的响应历程。通过调节结构使光纤探头垂直对准到钢筒爆心环面的外表面。

共进行了5 次对比实验，实验的工况设置如表1 所示。爆炸当量为10 g TNT 的实验中使用壁厚为6 mm 的钢筒。在钢筒内部增加10 mm厚的泡沫铝内衬后，钢筒爆心环面结构响应的位移峰值小幅度增加；2 次实验后钢筒爆心环面均无径向膨胀，可认为钢筒变形较小、处于弹性状态。爆炸当量约为180 g TNT 时，实验后钢筒存在明显的径向膨胀：在没有泡沫铝内衬的情况下，测得爆心环面径向膨胀5.3 mm，计算对应的塑性变形约为8.4%；泡沫铝内衬厚度为5 mm 时，爆心环面径向膨胀达到7.8 mm，对应的塑性变形约为12.4%；泡沫铝内衬厚度为15 mm 时，实验后钢筒破裂，如图4 所示。综上，由爆炸当量为180 g TNT 的对比实验可知，泡沫铝内衬厚度增加后，不但没有起到防护作用，反而导致结构发生了更严重的破坏。

图2 实验用钢筒和泡沫铝内衬Fig. 2 The steel cylinder and aluminum foam lining used in the experiment

表1 实验设置和钢筒外表面的变形情况Table 1 The experimental setup and deformation of the outer surfaces of the steel cylinders

图3 泡沫铝静态应力应变曲线Fig. 3 Static stress-strain curve of aluminum foam

2 数值模拟

为深入分析实验数据，找出泡沫铝内衬起防护/破坏作用的机理和规律，使用商业软件建立如图5所示的二维轴对称模型。模型主要包括钢筒、泡沫铝内衬、炸药和空气4 部分。其中钢筒和泡沫铝内衬采用Lagrange 网格，网格尺寸为2 mm×2 mm。采用双线性本构模型描述钢的力学行为，材料密度为7 830 kg/m3，弹性模量为200 GPa，泊松比为0.3，切线模量为7.7 GPa。考虑到动态加载条件下，由于应变率效应钢的屈服强度会提高，模型中设置钢的屈服强度为700 MPa[19]。对泡沫铝内衬采用CRUSHABLE_ FOAM 模型，模型中采用如图3所示的应力应变曲线。

采用欧拉网格计算爆炸流场，其中对空气采用理想气体状态方程，初始密度为1.29 kg/m3，初始比内能为2.068×105J/kg。对炸药采用JWL 状态方程：

图5 用于数值模拟的二维轴对称模型Fig. 5 The two-dimensional axisymmetric model used in numerical simulation

模型的边界条件设置如图5 所示，下边界为轴对称边界，左边界为对称边界，上边界和右边界为自由流出边界。通过流固耦合方法计算钢筒和泡沫铝的动态响应，输出钢筒爆心环面变形的计算结果，图6对比了数值模拟与实验的变形历程曲线，表1 对比了爆炸当量为10 g TNT 时钢筒弹性变形的应变峰值和爆炸当量为180 g TNT 时钢筒的残余应变。对于爆炸当量为10 g TNT 的2 次实验，钢筒处于弹性变形阶段，数值模拟得到的挠度曲线的第1 个峰值和后期振动趋势与实验数据基本一致。对于爆炸当量约为180 g TNT 的3 次实验，钢筒中未添加泡沫铝的实验仅有实验后测得的残余应变数据[16]，数值模拟得到的钢筒变形与实验数据基本一致。添加厚度为5 mm 的泡沫铝内衬时，数值模拟得到的钢筒外表面变形历程与实验数据基本一致。添加厚度为15 mm 的泡沫铝内衬时，实验中钢筒破碎，通过数值模拟预计的钢筒外壁位移达到13 mm；同时计算的钢筒外壁塑性变形达到19.6%，内壁塑性变形达到25.7%，超过Q345R 材料的断裂延伸率[20]，因而钢筒破碎。图6 中还比较了10 g TNT 实验中泡沫铝塑性变形分布的实验和数值模拟结果，两者基本吻合。综上，数值模拟结果与实验数据具有较好的一致性，计算模型基本可靠。

3 讨 论

本部分将结合实验结果，讨论弹性变形和塑性变形2 种情况下，泡沫铝对钢筒变形的影响机制，解释实验中泡沫铝内衬导致钢筒变形增大的原因。

对于弹性变形情况，结合表1 中爆炸当量为10 g TNT 的实验，计算了在壁厚为6 mm 的实验钢筒内不添加内衬、添加厚度为10 mm 的泡沫铝内衬、添加厚度为18 mm 的泡沫铝内衬时，爆心环面的内壁载荷和外壁位移曲线，结果如图7 所示。添加泡沫铝内衬后，钢筒内壁受到载荷峰值显著降低，但载荷的作用时间显著拉长。随着泡沫铝厚度增大，载荷到时逐渐提前，原因是波在固体介质中传播速度更快；但受泡沫铝塑性变形过程的影响，载荷峰值的到时越来越晚。另一方面，随着泡沫铝内衬厚度增大，钢筒外壁径向膨胀先增大后减小，说明泡沫铝内衬厚度足够时，才能对钢筒起到保护作用；相反，如果泡沫铝内衬厚度不足，还可能导致钢筒发生更严重的变形。图8 对比了爆炸当量为6～10 g TNT 时，无泡沫铝内衬、泡沫铝内衬厚度分别为10 和18 mm 时爆心环面外壁的径向膨胀位移峰值，进一步证明了上述泡沫铝内衬厚度对加固效果的影响规律。

图7 泡沫铝内衬对钢筒内壁载荷历程的影响(10 g TNT)Fig. 7 Influences of aluminum foam linings on the pressure loading on the inner surfaces of steel cylinders (10 g TNT)

图8 泡沫铝内衬对钢筒外壁弹性径向膨胀的影响Fig. 8 Influences of aluminum foam linings on the elastic radial expansion of the outer walls of steel cylinders

为研究塑性变形条件下，泡沫铝内衬对钢筒变形的影响，仍然以壁厚6 mm 的实验钢筒为对象，比较爆炸当量为8～140 g TNT，钢筒内添加10、18 mm 厚的泡沫铝内衬和不添加内衬时的钢筒塑性变形的计算结果，如图9 所示。爆炸当量超过10 g TNT 后，钢筒开始进入塑性变形状态；随着爆炸当量的增大，泡沫铝对钢筒变形的影响分为3 种模式。

图9 泡沫铝内衬对钢筒塑性变形的影响Fig. 9 Influences of aluminum foam linings on plastic deformation of steel cylinders

模式1 发生在钢筒刚刚开始发生塑性变形时。如图9 所示，爆炸当量为14 g TNT 时，添加10 mm 厚泡沫铝内衬的钢筒塑性变形增大，而添加18 mm 厚泡沫铝内衬的钢筒未发生塑性变形。这说明泡沫铝内衬厚度足够时，泡沫铝内衬可通过塑性变形吸收能量，起到理想的防护效果。

随着爆炸当量增大，泡沫铝内衬对钢筒变形的影响变为模式2：此时爆炸载荷的强度超过了泡沫铝内衬变形吸能的能力，泡沫铝内衬不再起到防护效果，反而导致钢筒塑性变形增大。由图10(a)可见，爆炸当量为30 g TNT 时，添加厚度为10 和18 mm 的泡沫铝内衬后，受泡沫铝内衬塑性变形的影响，钢筒内壁应力曲线的峰值到时延迟，应力峰值提高，且载荷的持续时间延长、载荷冲量大幅增加，因此钢筒塑性变形增大。模式2 中，可将泡沫铝内衬视为冲击波传播的可压缩介质，由于冲击波在固体中传播衰减弱于在空气中，因此钢筒变形增大。

继续增大爆炸当量，泡沫铝内衬对钢筒变形的影响进入模式3：泡沫铝厚度越大，钢筒的塑性变形越小。由图10(b)可见，爆炸当量为100 g TNT 时，泡沫铝内衬仅影响钢筒内壁载荷峰值到时，对载荷峰值和峰值后的衰减历程影响较小。但在特征相近的载荷下，泡沫铝内衬越厚，结构整体的质量就越大，因而钢筒的塑性变形减小。

图10 泡沫铝内衬对钢筒内壁载荷历程的影响Fig. 10 Influences of aluminum foam linings on pressure loading on the inner surfaces of steel cylinders

4 结 论

为在不改变钢筒尺寸的前提下提高钢筒的抗爆性能，通过实验和数值模拟分析了泡沫铝内衬对承受内部爆炸载荷钢筒变形的影响规律，并发现只有当添加足够厚度的泡沫铝内衬时才能减小钢筒变形；如果泡沫铝内衬厚度不足，反而会导致钢筒发生更大的变形，甚至破裂。根据进一步的数值计算结果可以发现，对于固定尺寸的含泡沫铝内衬钢筒，随着爆炸当量和载荷强度提高，泡沫铝内衬对钢筒内壁载荷特征的影响机理不同，泡沫铝内衬对钢筒变形的影响可能包含3 种模式：

模式1，泡沫铝通过塑性变形吸收能量，减小钢筒变形；

模式2，爆炸载荷强度超过泡沫铝变形吸能能力，此时泡沫铝内衬将导致钢筒内壁载荷峰值提高、持续时间延长，钢筒变形增大或破坏；

模式3，泡沫铝内衬对载荷特征的影响较小，但可增大结构质量，减小钢筒塑性变形。