Rock burst prediction based on data preprocessing and improved sparrow algorithm

在当今世界的大环境下，无论是民用领域还是军事领域，爆炸冲击的防护都变得越来越重要。传统防护爆炸冲击的方法是提高单层金属防护介质的厚度，但是单纯增加防护介质的厚度不仅会增加成本，还会增加整个防护结构的重量，降低装备的机动性。因此，研究由多层材料所组成的防护结构具有非常重要的意义。聚脲弹性体作为一种新兴的聚合物材料，具有良好的抗爆炸冲击性能，被广泛应用在建筑和舰艇等大型防护结构中。美国海军陆战队研究部门通过研究后发现，在军车底部喷涂聚脲可以在很大程度上减少爆炸冲击对人员和车辆的破坏^[1]。宋彬等^[2]对比了聚脲弹性体夹层、无夹层和橡胶夹层3种防爆罐在相同质量炸药爆炸载荷作用下的抗爆性能，发现聚脲弹性体夹层防爆罐的抗爆性能优于其他两种防爆罐。甘云丹等^[3]对比分析了单层钢板与聚脲涂覆钢板在水中爆炸下的抗爆性能，发现聚脲弹性体可以明显提高钢板的抗爆能力。Bahei-El-Din等^[4]发现，在用纤维复合材料为夹层、泡沫为芯层的传统三明治结构中加入聚脲层可以提高其防护性能。Amini等^[5]通过分析单位厚度的破坏能量后得出：在钢板的背面涂覆聚脲时，可以较明显地提高钢板的抗爆性能。赵鹏铎等^[6]对不同涂覆方式下单钢板和箱体结构的抗爆性能进行了较详细的试验研究，结果表明：在等面密度条件下，钢板迎爆面涂覆聚脲并不能提高其抗爆性能；而在等钢板厚度条件下，在单层钢板和箱体结构上涂覆聚脲能够提高其抗爆性能，且聚脲涂覆在钢板背面时其抗爆性能更好。Tekalur等^[7]对由乙烯酯玻璃纤维、聚脲等材料组成的复合结构进行了抗爆性能研究，结果表明相较于单一的乙烯酯玻璃纤维板，聚脲涂层可以明显提高其抗爆性能，且采用聚脲作为夹心层的三明治结构具有最好的抗爆性能。LeBlanc等^[8]研究了聚脲层的厚度和位置对聚脲、E-Glass/Epoxy复合结构防护性能的影响，发现相对于聚脲涂覆在受冲击测，聚脲涂覆在背侧更能大幅度提高防护结构的抗爆性能。Ackland等^[9]采用局部爆炸加载方法研究了相同面密度聚脲涂覆钢板复合结构的抗爆性能，结果表明单一钢板结构的抗爆性能优于聚脲涂覆结构的抗爆性能。大多数学者的研究显示，在等厚度钢板上涂覆聚脲层可以在一定程度上提高钢板的抗冲击性能，但是定量地研究爆心距、炸药质量以及涂覆聚脲厚度对聚脲涂覆钢板结构抗爆性能的响应规律还比较少。

本研究利用LS-DYNA有限元软件分析聚脲涂覆位置对钢板抗爆性能的影响，并结合量纲分析方法讨论爆心距、炸药质量、涂覆聚脲厚度对钢板抗爆性能的响应规律，以期为聚脲弹性体涂覆结构抗爆性能的评估与防护结构设计提供参考。

1.   理论分析

聚脲涂覆钢板在爆炸载荷作用下的响应涉及多个物理量，可通过量纲分析方法得到各物理量之间的关系。图1为球形TNT炸药与聚脲涂覆钢板的相互作用示意。

图  1  球形TNT与聚脲涂覆钢板

Figure  1.  Schematic of spherical TNT and polyurea coated steel plate

下载: 全尺寸图片 幻灯片

在没有发生断裂的情况下，钢板在爆炸载荷作用下的最大位移可作为其抗爆性能的重要指标之一，而决定最大位移ω的物理量包含以下几个方面。

（1）球形炸药的物理量参数，包括：炸药量Q，装药密度ρ_e、单位质量炸药所释放的化学能E_e、爆轰产物的膨胀指数γ_e。

（2）钢板的物理量参数，包括：密度ρ_s、弹性模量E_s、泊松比v_s、屈服极限σ_s、厚度D、钢板尺寸l。

（3）聚脲的物理量参数，包括：密度ρ_p、弹性模量E_p、泊松比v_p、拉伸强度σ_p、厚度d、尺寸l_p。

（4）钢板与聚脲之间的粘结力，包括：拉伸强度σ、剪切强度τ。

（5）空间几何参数：球形装药质心到靶板的垂直距离h。

基于以上各物理量参数，钢板最大位移ω与各物理量之间的关系函数可写为

取装药密度ρ_e、单位质量炸药所释放的化学能E_e、钢板厚度D作为基本量，（1）式可以化为

如果保持钢板、聚脲和炸药的材料性质以及钢板的几何尺寸不变，即

(1)式就能变为更简单的无量纲关系式

当保持聚脲的几何尺寸及厚度、钢板与聚脲之间的粘结力、炸高不变时，(4)式化为简单的无量纲关系

当保持聚脲的尺寸及厚度、钢板与聚脲之间的粘结力、炸药量不变时，(4)式可变为简单的无量纲关系，即

当保持聚脲涂层的尺寸l_p、钢板与聚脲之间的粘结力、炸药量、炸高不变时，(4)式可化为

上面几个式子表明，当保持钢板与聚脲之间的粘结力不变时，钢板的最大位移只和炸高、炸药质量、聚脲涂层厚度有关，具体的无量纲关系式可通过试验或数值模拟确定。下面将用一系列数值仿真来确定相关函数关系。

2.   有限元计算模型

2.1   仿真建模及验证

为了验证所用模拟方法的有效性，基于LS-DYNA有限元软件对文献[6]中爆炸载荷作用下聚脲涂覆钢板的试验结果进行仿真分析。试验中聚脲涂覆底材为边长400 mm的方形Q235A钢板，板厚1.2 mm，钢板迎爆面涂覆6 mm厚聚脲层。聚脲涂覆钢板结构固定在两个夹板之间，40 g柱形炸药在50 mm爆心距下对聚脲涂覆钢板的加载区域为250 mm×250 mm，更详细的试验方案见文献[6]。

在数值模拟建模中，由于爆炸载荷和聚脲涂覆钢板结构具有对称性，所以只建立1/4模型，并且在模型的边界处施加固支约束和对称约束。试验中聚脲涂覆钢板的变形区域为250 mm×250 mm，为了节省计算时间，分析时只考虑聚脲涂覆钢板的实际变形区域。钢板、聚脲层用Lagrange网格离散，聚脲层采8节点实体单元SOLID164，而钢板采用壳单元的建模方式。聚脲涂覆钢板与炸药的有限元模型及尺寸如图2所示。

图  2  聚脲涂覆钢板结构与炸药模型

Figure  2.  Schematic of polyurea-coated structure and explosive

下载: 全尺寸图片 幻灯片

Q235A钢采用动力学塑性本构模型（*MAT_PLASTIC_KINEMATIC）。由于聚脲是一种超黏弹性体，具有很明显的应变率效应，故采用6参数的Mooney-Rivlin模型^[10]，具体的材料参数列于表1。

表  1  聚脲弹性体材料参数^[6]

Table  1.  Parameters of polyurea elastomer material^[6]

Density/(kg·m–3)	Tensile strength/MPa	Elongation at break/%	Tear strength/(kN·m–1)	Shore hardness (HA)
1020	24	400	85	85–95

下载: 导出CSV 

| 显示表格

用关键字*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_TIEBREAK^[11]定义钢板与聚脲层之间的黏结力，设置两种材料之间的黏结强度。计算中当钢板与聚脲界面上的应力大于它们之间的黏结强度时，钢板与聚脲发生分层破坏。

图3给出文献[6]中工况3的靶板变形情况，图4为相应的数值模拟位移云图。可见，中心处最大位移的试验结果与仿真结果吻合得非常好，验证了本研究所采用模型的有效性。

图  3  试验中聚脲涂覆钢板的变形^[6]

Figure  3.  Deformation of polyurea coated steel plate in test^[6]

下载: 全尺寸图片 幻灯片

图  4  数值模拟中聚脲涂覆钢板的变形

Figure  4.  Deformation of polyurea coated steel plate in simulation

下载: 全尺寸图片 幻灯片

2.2   计算模型

首先对聚脲与钢板组成的具有相同面密度的单一钢板结构、正面涂层、背面涂层和正背面均有涂层4种结构在爆炸载荷作用下的动态响应进行分析，并从钢板的最大位移、塑性应变两个角度分析不同结构的抗爆性能，具体结构设计方案见表2。复合结构的尺寸均为800 mm×800 mm（方形板），钢板厚度为3 mm，聚脲涂层的厚度为4 mm，球形炸药位于方形板中心位置的正上方，炸高为200 mm。

表  2  结构设计方案

Table  2.  Structural design schemes

Condition number	Structural diagram	Steel plate thickness/mm	Polyurea thickness/mm	TNT dose/kg	Detonation distance/mm
A1	Contrast condition	3	0	1	200
A2	Single steel plate	3.52	0	1	200
A3	Front coating	3	4	1	200
A4	Back coating	3	4	1	200
A5	Double sided coating	3	2+2	1	200
Note: Q235A steel plate; Polyurea.

下载: 导出CSV 

| 显示表格

钢板为Q235A钢，采用动力学塑性本构模型（*MAT_PLASTIC_KINEMATIC），强化方式为随动强化，并且考虑材料的应变率效应，其屈服应力放大系数的计算公式为

式中：σ_Y为屈服应力，σ_d为动态屈服应力，${\dot \varepsilon }$为应变率，P和C为材料应变率效应常数。Q235A钢的具体参数见表3。

表  3  Q235A钢材料参数^[12]

Table  3.  Parameters of Q235A steel material^[12]

Density/(kg·m–3)	Young’s modulus/GPa	Poisson’s ratio	σY/MPa	Tangent modulus/GPa	C/s–1	P
7850	210	0.3	240	2	105	5

下载: 导出CSV 

| 显示表格

2.3   复合结构抗爆性分析

上述5种防护结构在相同当量爆炸载荷作用下的变形情况如图5所示。球形炸药位于防护结构中心位置的正上方，所以钢板中心区域的位移最大。为了更好地对比钢板的最大位移和塑性变形情况，提取钢板中心区域单元的位移、塑性应变随时间的变化，如图6、图7所示。

图  5  不同工况中钢板的位移云图

Figure  5.  Displacement cloud chart of steel plate in different working conditions

下载: 全尺寸图片 幻灯片

图  6  钢板最大变形处的位移时程曲线

Figure  6.  Displacement-time curves at maximum deflections of steel plates

下载: 全尺寸图片 幻灯片

图  7  钢板最大变形处的塑性应变时程曲线

Figure  7.  Plastic strain-time curves at maximum deflections of steel plates

下载: 全尺寸图片 幻灯片

由图6、图7可知，纯钢板、无涂覆层的A1工况中的钢板变形大于其余工况，且A3>A5>A4>A2， A3、A5、A4、A2工况中钢板的最大变形分别比A1减小了7.0、7.7、7.9、14.0 mm。比较钢板中心区域的最大塑性应变得到，A1>A3>A5>A4>A2，且A3、A5、A4、A2工况中钢板的最大塑性应变分别比A1减小了11.2%、12.4%、15.5%、16.9%。这说明在等厚度钢板上涂覆一定厚度的聚脲材料能够提高钢板的抗爆性能，并且聚脲涂覆在钢板背爆面的抗爆效果优于其他两种方式的涂覆。背面涂层结构受到爆炸冲击时，聚脲层的变形、破裂等弥散了较多冲击能量，从而减小钢板的变形；而当聚脲涂覆在迎爆面时，聚脲层首先受到爆炸冲击，聚脲层与钢板发生大面积的分层现象，且聚脲层发生撕裂脱落，从而在一定程度上减弱了爆炸冲击对钢板的作用：这两种涂覆方式都对爆炸冲击能量有一定的耗散作用，只是两者相比，背面涂覆方式更好。在等面密度条件下（A2～A5），纯钢板结构比聚脲涂覆钢板复合结构的抗爆性能更好，说明聚脲更应该涂覆在已有结构表面而不是用于设计等面密度的复合结构。

3.   聚脲涂覆钢板最大位移的影响参数分析

由2.3节可知，在等厚度钢板上涂覆一定厚度的聚脲层能够提高钢板的抗爆性能，且聚脲涂覆在钢板背面的方式更好。本节以A4工况为例，用量纲方法分析爆心距、炸药质量、聚脲层厚度对聚脲涂覆钢板变形的影响。

3.1   爆心距对聚脲涂覆钢板最大位移的影响

保持炸药质量、聚脲层厚度不变，改变球形炸药质心距钢板中心位置的距离(0.2～0.8 m)。首先对得到的钢板中心区域最大位移数据进行无量纲化，再用指数函数对无量纲数据进行拟合，得到如图8所示的无量纲化关系曲线。

图  8  爆心距对聚脲涂覆钢板最大位移的影响

Figure  8.  Relationship between the detonation distance and the maximum displacement of polyurea coated steel plate

下载: 全尺寸图片 幻灯片

从图8中可以看出，在保持其他变量不变的条件下，随着炸药质心与钢板之间距离的增加，钢板的最大位移近似呈指数形式递减，拟合关系式为$\omega/D=2.55+40.3\,{\rm {e}}^{-0.005h/D}$。由球形炸药在无限空气介质中爆炸的超压、比冲量经验公式可知，当保持炸药质量不变时，随着爆心距的增加，结构受到的超压峰值和冲量下降，爆炸载荷传递到结构上的能量也随之减少，因此聚脲涂覆钢板的最大变形量随爆心距的增加而近似呈指数形式递减。

3.2   炸药质量对聚脲涂覆钢板最大位移的影响

同样以A4工况为例，保持涂层聚脲厚度和炸高不变，改变炸药的质量，分析炸药质量对聚脲涂覆钢板最大位移的影响规律。和3.1节一样，首先对得到的聚脲涂覆钢板最大位移数据进行无量纲化，再用线性函数拟合无量纲数据，得到如图9所示的关系曲线。

图  9  炸药质量对聚脲涂覆钢板最大位移的影响

Figure  9.  Relationship between the mass of explosive and the maximum displacement of polyurea coated steel plate

下载: 全尺寸图片 幻灯片

从图9中可以看出，在保持其他变量不变的情况下，随着球形炸药质量的增加，钢板中心区域的最大位移近似呈线性增加趋势，拟合线性关系式为$\omega/D=2.33 \sqrt[{^{^3}}\!\!\!\!]{\dfrac{Q}{{\rho}_{\rm {e}}D^{3}}}-32.75$。同样从球形炸药在无限空气介质中爆炸的超压、比冲量经验公式可知，当保持爆心距不变时，随着炸药质量的增加，爆炸载荷作用在结构上的超压和比冲量增大，传递到复合结构上的能量也随之增加，因此聚脲涂覆钢板的最大位移随炸药质量的增加近似呈线性增长趋势。

3.3   聚脲厚度对聚脲涂覆钢板最大位移的影响

保持炸药质量、炸高不变，改变涂覆在钢板背面的聚脲层的厚度。同样对该条件下得到的所有钢板的最大位移数据进行无量纲化，用线性函数对无量纲数据进行拟合，得到如图10所示的关系曲线。

图  10  聚脲层厚度对聚脲涂覆钢板最大位移的影响

Figure  10.  Relationship between the thickness of polyurea layer and the maximum displacement of polyurea coated steel plate

下载: 全尺寸图片 幻灯片

从图10中可以看出，在保持其他变量不变的情况下，随着涂覆在钢板背面聚脲层厚度的增加，钢板的最大位移近似呈线性减小趋势，拟合关系式为$\omega/D=33.98-1.372(d/D)$。钢板的最大变形随着涂覆聚脲层厚度的增加而减小的主要原因是：增加聚脲层厚度在一定程度上增加了结构的刚度，从而整体上提高了结构的抗爆性能；将聚脲涂覆在背面时，由于聚脲材料的黏弹性特性，复合结构在受到爆炸冲击后，一部分以破裂碎片动能的形式带走大量冲击能量，一部分能量被聚脲层弥散和吸收，且聚脲层越厚，其破裂的程度越严重，破片的数量越多，以动能形式耗散的能量越大，所以聚脲层越厚，其吸收和弥散的能量越多。

4.   结　论

通过数值仿真分析了聚脲涂覆方式对钢板抗爆性的影响，结合量纲分析方法讨论了爆心距、炸药质量、涂覆聚脲厚度对钢板中心最大位移的影响，得到以下结论：

（1）在等厚度钢板上涂覆聚脲能够提高钢板整体的抗爆性能，且聚脲涂覆在钢板背爆面的效果优于其他两种方式，而在等面密度条件下，纯钢板结构比聚脲涂覆钢板复合结构的抗爆性能更好；

（2）在保持其他变量不变的情况下，随着爆心距的增大，聚脲涂覆钢板的最大位移近似呈指数形式递减；

（3）改变炸药质量，聚脲涂覆钢板的最大位移近似随着炸药质量的增加而呈线性增加；

（4）改变涂覆在钢板背爆面的聚脲厚度，其他条件保持不变，则聚脲涂覆钢板的最大位移随着涂层厚度的增加近似呈线性递减趋势。