负压爆炸载荷作用下固支钢板变形研究

杨锐; 汪泉; 谢守冬; 李瑞; 涂唱畅; 徐小猛; 李孝臣

doi:10.11858/gywlxb.20230685

负压爆炸载荷作用下固支钢板变形研究

doi: 10.11858/gywlxb.20230685

杨锐^1,,
汪泉^1,*, ,,
谢守冬^{2, 3},
李瑞^{1, 4},
涂唱畅¹,
徐小猛¹,
李孝臣¹

1.
安徽理工大学化学工程学院, 安徽淮南　232001
2.
安徽理工大学土木建筑学院, 安徽淮南　232001
3.
宏大爆破工程集团有限责任公司, 广东广州　510000
4.
安徽理工大学煤炭安全精准开采国家地方联合工程研究中心, 安徽淮南　232001

基金项目: 国家自然科学基金（11872002）；精细爆破国家重点实验室开放课题（PBSKL-2022-B-05）；煤炭安全精准开采国家地方联合工程研究中心开放基金（EC2021015）

详细信息

作者简介:
杨　锐（1999－），男，硕士研究生，主要从事爆炸毁伤研究. E-mail：2240036819@qq.com

通讯作者:
汪　泉（1980－），男，博士，教授，博士生导师，主要从事爆炸理论与技术研究. E-mail：wqaust@163.com

中图分类号: O383.1
计量
- 文章访问数: 217
- HTML全文浏览量: 105
- PDF下载量: 87
出版历程
- 收稿日期: 2023-06-21
- 修回日期: 2023-07-11
- 网络出版日期: 2023-09-15
- 刊出日期: 2023-11-07

Deformation of Fixed Support Steel Plate under Explosion Load in Negative Pressure Environment

YANG Rui^1
,,
WANG Quan^{1,*
, ,},
XIE Shoudong^{2, 3},
LI Rui^{1, 4},
TU Changchang¹,
XU Xiaomeng¹,
LI Xiaochen¹

1.
School of Chemical Engineering, Anhui University of Science and Technology, Huainan 232001, Anhui, China
2.
School of Civil Engineering and Architecture, Anhui University of Science and Technology, Huainan 232001, Anhui, China
3.
Hongda Blasting Engineering Group Co., Ltd., Guangzhou 510000, Guangdong, China
4.
Joint National-Local Engineering Research Centre for Safe and Precise Coal Mining, Anhui University of Science and Technology, Huainan 232001, Anhui, China

摘要

摘要: 为研究负压爆炸载荷作用下结构的动态响应，以固支钢板为防护工程的简化单元，开展了负压爆炸实验，探究固支钢板在负压爆炸载荷作用下的变形规律，分析不同负压环境下固支钢板的极限应变和失效条件。采用AUTODYN对负压爆炸载荷作用下固支钢板的动态响应进行数值模拟，通过对比实验结果，验证了数值模拟结果的准确性。结果表明：随着初始环境压力的下降，相同爆距下钢板中心点的最大挠度和最大速度减小；负压爆炸载荷作用下，钢板整体出现塑性大变形，迎爆面形成凹坑，钢板四周在垂直于固支边界指向钢板中心的方向上出现明显的拉伸变形，钢板边缘区的挠度变化基本相同，中心点的最大挠度随着环境压力的下降而减小。通过双向应变假设，确定了钢板的动态极限应变为0.269。建立了负压环境下炸药爆炸冲击波的反射比冲量公式，并对基于刚塑性假设和能量准则提出的失效判据进行检验。研究结果可为负压环境下爆炸空气冲击波威力等效评估、高原环境下目标毁伤评估提供参考。
- 负压 /
- 固支方板 /
- 爆炸载荷 /
- 乳化炸药 /
- 拉伸断裂
Abstract: In order to study the dynamic response of the structure under explosion load in negative pressure environment, the negative pressure explosion experiments were carried out for the fixed supported steel plate, which is as a simplified unit of the protection project. The deformation law, and the ultimate strain and failure conditions of the fixed supported steel plate under different negative pressures were analyzed. The numerical simulation of the dynamic response of the fixed supported steel plate under negative pressure explosion load was carried out by AUTODYN, and the accuracy of the numerical simulation results was verified by comparing the experimental results. The results show that when the initial ambient pressure decreases, for the same burst distance, both the maximum deflection and the maximum velocity at the center point of the steel plate decrease. Under the negative pressure explosion load, the steel plate produces large plastic deformation, the oncoming surface forms a pit, and obvious tensile deformation occurs at the edges perpendicular to the boundary direction. The deflection changes in the edge zone are basically the same, and the maximum deflection at the center point decreases with the decrease of environment pressure. Through the bidirectional strain assumption, the dynamic limit strain of the steel plate is determined to be 0.269. The reflectance specific impulse formula of explosive blast wave under negative pressure environment was established, and the failure criterion based on rigid-plastic hypothesis and energy criterion were examined. The research results can provide a reference for the equivalent evaluation of the shock wave power of explosive air in negative pressure environment and the target damage assessment in plateau environment.
- negative pressure /
- fixed supported square plate /
- explosion load /
- emulsion explosive /
- tensile fracture

HTML全文

近年来，高原环境的军事演习增多，高原环境下的爆炸冲击毁伤问题得到国内外的高度重视。在处理飞机、船舶、车辆等复杂结构在爆炸作用下的毁伤问题时，可将这些复杂结构视为装甲类目标，以其表面的板结构作为研究对象，通过实验、数值模拟和理论分析开展研究^[1–2]。Jacob等^[3]研究了爆距和药量对爆炸载荷作用下固支圆板响应的影响。陈长海等^[4]通过实验分析了方板在近爆载荷作用下的破坏模式，提出了破裂判别条件，从而预测了爆炸载荷作用下钢板是否失效。韩璐等^[5]通过方板毁伤实验建立了不同毁伤模式下的数值模型。Wang等^[6]对聚异氰氨酸酯噁唑烷酮（polyisocyanateoxazodone，POZD）涂层钢板-钢筋混凝土板和普通钢板-钢筋混凝土板进行爆炸实验，分析了药量、POZD厚度、钢板厚度和钢筋混凝土板厚度对POZD涂层钢板-钢筋混凝土板抗爆性能的影响。Remennikov等^[7]采用液体球状炸药对钢板进行近爆实验，探讨了几种防护方案对实验结果的影响，并建立了工程模型。秦业志等^[8]采用RKDG（Runge-Kutta discontinuous Galerkin）方法研究了柱状装药水下近场爆炸作用下单层钢板的响应特性。以往的模拟实验大多针对水下爆炸或常压下爆炸对板结构的毁伤^[9–10]，对于高原、高空等负压环境下爆炸载荷与目标结构变形之间的关系研究较少。当炸药药卷从常压环境转移到负压环境时，为了平衡内外气压，药卷体积膨胀，密度降低，根据爆轰产物的JWL状态方程^[11]，装药密度降低会直接影响爆轰产物膨胀规律，同时，负压环境对冲击波的传播过程亦有影响，致使负压与常压下的毁伤实验结果有所不同。

本研究将对负压爆炸载荷作用下固支钢板的动态响应过程进行实验和数值模拟，以期为负压环境下爆炸冲击波威力等效评估、高原环境下目标毁伤评估等提供参考。

1. 负压爆炸实验

1.1 实验设计

负压爆炸实验在可调真空度的爆炸容器内进行。固支方板作为模拟复杂防护结构的简化单元，采用Q235钢制成，尺寸为1 mm×460 mm×460 mm，四周均匀分布20个直径为16 mm的螺栓孔。为实现方形钢板的边界固定条件，在爆炸容器中心平台上搭建钢制方形支座^[4]，如图1(a)所示。由于存在夹板，因此实验中钢板的有效迎爆面尺寸为300 mm×300 mm。炸药采用乳化炸药，爆热为3610.20 J/g，TNT当量约为乳化炸药质量的0.8629倍^[12]。用绳索将乳化炸药吊置在钢板上方，将容器抽至预期负压后静置2 h，将导爆管雷管插入炸药中心起爆。实验系统如图1(b)所示，实验工况如表1所示，其中p_e为初始环境压力，W为乳化炸药质量。

图 1 实验示意图

Figure 1. Schematic diagram of the experiment

下载: 全尺寸图片幻灯片

表 1 实验工况

Table 1. Experimental conditions

Case	p_e/kPa	W/g	TNT equivalence/g	Explosion distance/mm
1	60	150	130	150
2	80	150	130	150
3	101	150	130	150
4	101	200	170	50

下载: 导出CSV

| 显示表格

1.2 固支钢板在不同负压环境下的变形破坏分析

图2给出了工况1～工况3下固支钢板的变形情况。可以看出，不同负压爆炸载荷作用下固支钢板均未穿孔，整体出现塑性大变形，迎爆面形成凹坑，即毁伤模式为塑性大变形毁伤^[5]。当初始环境压力从101 kPa降到60 kPa时，钢板中心处凹陷程度下降。钢板四周固支端螺栓孔在垂直于固支边界并指向钢板中心的方向上出现明显的拉伸变形，且随着环境压力的下降，钢板四周的拉伸变形逐渐减小。钢板迎爆面对角线处出现明显的塑性铰线，塑性铰线延伸至中心，中心点为凹坑最深处，钢板表面无明显的爆炸产物灼烧现象。

图 2 实验后钢板的形貌

Figure 2. Morphology of the steel plate after the experiment

下载: 全尺寸图片幻灯片

在钢板中线上每隔1 cm布置一个测点，测量钢板中线变形挠度γ。设钢板有效迎爆面的半宽为L，以钢板中心为原点、钢板中线为x轴，得到工况1～工况3中钢板中线的变形轮廓，如图3所示。在10～15 cm区间（边缘区），3种工况下钢板的变形程度相近；在0～10 cm区间（中心区），3种工况下钢板的变形程度不同，最大挠度位于钢板距离炸药中心的最近点（即钢板中心），初始环境压力为101、80、60 kPa时，最大挠度分别为37、33、31 mm，最大挠度随初始环境压力的下降而减小。

图 3 钢板变形轮廓

Figure 3. Deformation contours of steel plates

下载: 全尺寸图片幻灯片

目前，通常用超压-冲量准则^[13]衡量空气爆炸载荷对目标的破坏效应，即认为空气冲击波阵面的最大超压Δp_m与比冲量i的共同作用满足某临界条件时目标破坏。炸药爆炸瞬间，爆炸产物剧烈膨胀并压缩周围空气，形成压力、密度、温度突跃的空气冲击波，其初始强度取决于爆轰波参数、介质的密度及可压缩性等特性。负压环境的典型特征是环境压力和空气密度低于常压。随着爆炸容器内压力的降低，冲击波阵面前未扰动的空气变得稀薄，冲击波峰值超压和比冲量减小^[14]。因此，负压环境通过影响冲击波的传播过程对爆炸载荷起衰减作用。由于爆炸冲击波的比冲量决定实验后钢板的最大挠度^[15]，所以钢板的变形程度随着环境压力的下降而减小。同时，环境压力下降时，爆炸产物膨胀速率下降变慢，空气冲击波速度变快。

工况4增大了乳化炸药的药量并减小了爆距，可视为接触爆炸载荷作用。图4给出了实验后钢板的破坏形貌。钢板中心部位呈花瓣状开裂，3条长度分别为19.20、12.50、15.70 cm的裂纹贯穿整个钢板，将钢板分成3个破片，裂纹延伸至四周固支边界使钢板剪切断裂，有效抗爆面整体剪切失效，从支座上脱落并扭曲成团，毁伤模式为三花瓣状破口毁伤^[5]。固支边界处的剪切断面光滑整齐，垂直于剪切断面方向出现明显的拉伸变形，沿着剪切断面方向出现细微的拉伸变形。相对于初始板厚，中心部位的裂纹边缘呈现明显的减薄现象。

图 4 工况4中钢板的破坏模式

Figure 4. Damage mode of steel plates in case 4

下载: 全尺寸图片幻灯片

1.3 动态极限应变

工况4中，钢板中心部位的裂纹由拉伸断裂所致，边缘出现明显的减薄现象。根据双向应变假设，钢板的动态极限应变可由钢板开裂处的减薄率计算^[16]。在开裂处取极小微元，其厚度为h，体积为hdxdy，设微元在x、y方向的应变分别为ε_x、ε_y，厚度方向由h变成h₁，根据体积不变原理，有

$h{\text{d}}x{\text{d}}y = (1+{\varepsilon _x})(1+{\varepsilon _y}){h_1}{\text{d}}x{\text{d}}y$

(1)

假设微元在x、y方向上的应变是对称的，即ε_x=ε_y，可得

${\varepsilon _{\text{f}}} = {\varepsilon _x} = {\varepsilon _y} = \sqrt {\frac{h}{{{h_1}}}} - 1$

(2)

式中：ε_f为结构的动态极限应变。

测量3个破片在开裂处的钢板厚度，如图4(b)所示，共取8个测点，实验和计算得到的减薄率以及双向极限应变列于表2。考虑到钢板结构的受力和变形比较复杂，测点3和测点8的数据与其他测点的数据相差较大，因此将其舍弃，取其余6个测点的双向极限应变的平均值作为钢板的动态极限应变，得到Q235钢板的动态极限应变为0.269。由于结构的动态极限应变与材料属性、结构形式、加载方式等诸多因素有关，本研究所用薄钢板的尺寸较小，因此得到的双向极限应变0.269是合理的。该结果仅适用于负压爆炸载荷作用下尺寸较小的固支薄钢板中心拉伸断裂情况，对Q235钢结构在类似的爆炸载荷作用下的失效情况具有一定的参考价值。

表 2 开裂处钢板厚度的测量结果

Table 2. Measurement results of steel thickness at crack

Measuring point	Thickness/ mm	Down gauging rate	Bidirectional ultimate strain	Measuring point	Thickness/ mm	Down gauging rate	Bidirectional ultimate strain
1	0.60	0.40	0.291	5	0.64	0.36	0.250
2	0.62	0.38	0.270	6	0.54	0.46	0.361
3	0.52	0.48	0.387	7	0.66	0.34	0.231
4	0.68	0.32	0.213	8	0.70	0.30	0.195

下载: 导出CSV

| 显示表格

1.4 失效条件

近场爆炸时，一般按有效冲量计算爆炸对目标的破坏作用。与单质炸药相比，乳化炸药在负压环境下的做功能力和破坏能力较弱，传统的基于TNT等单质炸药提出的正压冲量公式无法准确预测乳化炸药爆炸冲击波的比冲量。文献[12]给出了乳化炸药在不同负压环境下的入射波正压比冲量公式

${i_{\text{s}}} = 162.5\frac{{\sqrt[{3}]{{{W^2}}}}}{R}\left( \frac{p_{\rm h}}{p_0}\right)^{1/3}$

(3)

式中：装药质量W的单位为kg；R为比例距离，m；p_h和p₀分别为低压环境和标准大气压环境下的环境压力，kPa；i_s为入射波正压比冲量，Pa·s。计算刚性表面的爆炸冲击波时，其药量应按照空中爆炸时药量的2倍^[17]计算，若周围介质可吸收能量，则式(3)中的装药质量应当用2θW表示，其中：θ是考虑爆炸能量耗散在空气中的份额所得出的材料系数，对于本研究中的钢板，θ取1。入射波正压比冲量公式转化为

${i_{\rm {s}}} = 258\frac{{\sqrt[3]{{{W^2}}}}}{R}\left( \frac{p_{\rm {h}}}{p_0}\right)^{1/3}$

(4)

文献[18]给出了冲击波正规反射比冲量i_r与入射波正压比冲量i_s的关系式，即i_r = i_s(1+cos φ₀)，其中：φ₀为爆炸冲击波到钢板的入射角。本实验中，爆心距r=150 mm，钢板半宽L=150 mm，则0° ≤ φ₀ < 45°。因此，不同环境压力下炸药爆炸冲击波的正规反射比冲量为

${i_{\text{r}}} = 516\frac{{\sqrt[3]{{{W^2}}}}}{R}\left( \frac{p_{\rm h}}{p_0}\right)^{1/3}$

(5)

假设结构为理想的刚塑性体，当受到爆炸冲击波作用后，结构产生位移，直至速度为零时，全部动能转化为应变能。根据动量定理，单位体积最大速度v₀ = i_r /ρh，其中ρ为质量密度。

由工况4可知，钢板中心发生拉伸断裂，中心处的应变率^[19]为

$\dot \varepsilon = \frac{{\sqrt {{\varepsilon _{\text{f}}}} {v_0}}}{L}$

(6)

本研究中， L=0.15 m。

设动屈服强度σ_d=ασ₀，其中：σ₀为准静态屈服强度，α为应变率系数。α由Cowper-Symonds关系得到

$\alpha = 1+{\left( \frac {\dot \varepsilon }{D} \right)^{{1 \mathord{\left/ {\vphantom {1 q}} \right. } q}}}$

(7)

式中：D为应变率常数，q为应变率指数。对于Q235钢，q=5，D=40.4 s⁻¹，σ₀=235 MPa。

若结构的最大应变ε_m接近甚至超过极限应变ε_f，则认为结构破裂失效。根据刚塑性假设和能量准则，钢板中心点的单位体积应变能为σ_dε_m。由于近爆时间极短，因此可以认为钢板获得的初始动能全部转化为应变内能

${\sigma _{\text{d}}}{\varepsilon _{\rm {m}}} = \frac{{\rho v_0^2}}{2}$

(8)

文献[4]中定义变量η作为结构的破裂判别参数，η的表达式为

$\eta = \frac{{{\varepsilon _{\text{m}}}}}{{{\varepsilon _{\rm f}}}} = \frac{{\rho v_0^2}}{{2{\sigma _{\text{d}}}{\varepsilon _{\text{f}}}}} = \frac{{i_{\text{r}}^2}}{{2\rho {\sigma _{\text{d}}}{h^2}{\varepsilon _{\text{f}}}}}$

(9)

失效应变采用1.3节中的双向极限应变，即ε_f =0.269。当η <1时，钢板仅发生塑性变形，未发生断裂失效；当η ≥1时，钢板发生断裂失效。

表3给出了工况1～工况4的失效判别条件参数。从表3可以看出：工况1～工况3中，η<1，说明3种负压环境下固支钢板均未穿孔，整体未断裂失效，与实验结果一致；工况4中，η>1，说明钢板断裂失效，与实验结果一致。

表 3 失效判别式的相关参数

Table 3. Relevant parameters for failure discriminant

Case	i_r/(Pa∙s)	v₀/(m∙s⁻¹)	$\dot{\varepsilon }$ /s⁻¹	α	σ_d/MPa	ε_m	η
1	742.1	94.5	326.8	2.52	592.2	0.059	0.219
2	816.8	104.1	359.9	2.55	599.3	0.071	0.264
3	882.8	112.5	389.0	2.57	604.0	0.082	0.305
4	3204.1	408.2	1411.4	3.04	714.4	0.915	3.401

下载: 导出CSV

| 显示表格

通过测量负压爆炸载荷作用下金属板开裂处的厚度，可以计算出金属板的动态极限应变。根据刚塑性假设，由正规反射比冲量公式得出金属板获得的动能，进而根据能量密度准则计算出金属板在爆炸作用下瞬间形变的最大应变。当最大应变超过动态极限应变，即η >1时，金属板断裂失效。当金属板处于近距离爆炸场时，已知金属板的性能参数，可以推算出负压环境下爆炸载荷作用下金属板断裂失效的最小起爆药量。值得注意的是，当选取的金属板厚度较薄，最大应变小于或接近动态极限应变时，金属板也可能断裂失效；当爆距过大时，爆炸载荷可视为均布爆炸载荷，金属板的断裂失效情况可能与局部爆炸载荷作用下的情况有所不同。

2. 数值模拟

2.1 一维空中爆炸及结果映射

采用结果映射Ramap技术，建立一维楔形空气域模型，模拟不同环境压力下TNT爆炸冲击波传播的初始过程，然后将一维结果映射到相同工况的三维模型中继续求解。通过映射可以提高计算效率，得到更精确的结果。一维楔形空气域模型及炸药网格如图5(a)所示。在欧拉空气域中设置半径为26.7 mm的TNT，比例距离为3.0 m/kg^1/3。如图5(b)所示，建立1/4三维模型，包含空气和钢两种物质。空气域的长和宽均为150 mm，高为240 mm，综合考虑网格敏感性和模拟时间成本，将网格尺寸设置为2 mm×2 mm^[20]，边界设为Flow-out透射边界。钢板迎爆面的长和宽均为150 mm，厚度为1 mm，钢板网格尺寸为1 mm×1 mm，在板边缘添加速度为零的约束条件。钢板采用Lagrange算法，板结构和空气域采用流固耦合算法。

图 5 几何模型

Figure 5. Geometrical model

下载: 全尺寸图片幻灯片

2.2 模型材料

TNT采用JWL状态方程描述，相关材料参数取自AUTODYN标准材料模型库。钢板选用Q235钢，采用Shock状态方程描述，其参数取自AUTODYN材料库，Johnson-Cook强度模型参数列于表4，其中：ρ_s为Q235钢的密度，A为初始屈服应力，B为硬化常数，n为硬化指数，C为应变率常数，m为热软化指数，T_m为熔化温度。

表 4 材料参数^[21]

Table 4. Material parameters^[21]

ρ_s/(g·cm⁻³)	A/MPa	B/MPa	n	C	m	T_m/K
7.85	293.8	230.2	0.578	0.0652	0.706	1795

下载: 导出CSV

| 显示表格

空气采用理想气体状态方程描述，即

$p = \left( {K - 1} \right)\rho_{\rm a} e$

(10)

式中：p为气压，K为绝热指数，ρ_a为空气密度， e为空气的初始比内能。不同压力环境下，空气的密度与压力之比相同，即

${\rho _i} = \frac{{{p_i}}}{{{p_j}}}{\rho _j}$

(11)

式中：下标i和j代表不同状态。

根据式(11)，通过改变空气的初始密度，可以实现不同的负压环境。表5列出了不同初始环境压力对应的空气密度。

表 5 不同环境压力对应的空气密度

Table 5. Air densities at different environment pressures

Pressure/kPa	Density/(kg·m⁻³)	Pressure/kPa	Density/(kg·m⁻³)
101	1.225	40	0.484
80	0.967	20	0.242
60	0.725

下载: 导出CSV

| 显示表格

2.3 数值模拟结果分析

图6给出了在不同负压环境下钢板的变形情况。不同负压环境下，钢板均出现明显的塑性大变形，未产生穿孔破损，钢板中心点的挠度最大。当初始环境压力为101、80、60、40、20 kPa时，钢板的最大挠度分别为34.7、32.9、30.6、28.3、27.9 mm；初始环境压力从101 kPa降到80、60、40、20 kPa时，钢板最大挠度分别减小5.3%、11.8%、18.4%和19.6%。随着初始环境压力的下降，冲击波传播过程中空气密度减小，爆炸冲击波的峰值超压和比冲量也减小，致使钢板整体形变程度呈下降趋势。

图 6 不同负压环境下钢板的变形情况

Figure 6. Deformation of steel plate under different negative pressure environments

下载: 全尺寸图片幻灯片

图7为钢板中心点速度曲线。当初始环境压力为101、80、60、40、20 kPa时，钢板中心点最大速度分别为113.7、106.6、83.5、80.3、70.5 m/s，随着初始环境压力下降，钢板中心点最大速度减小。钢板中心点速度在0.66 ms左右降至零以下，即开始反向运动，随后发生振荡，并逐渐衰减，最终停止运动。

图 7 钢板中心点速度曲线

Figure 7. Velocity curves at center point of steel plate

下载: 全尺寸图片幻灯片

数值模拟得到的钢板变形与实验结果一致，证明了数值模型的正确性。当初始环境压力为101、80、60 kPa时，钢板中心点挠度的实验结果与数值模拟结果的相对误差分别为6.2%、0.4%、1.3%，两者符合得较好，但数值模拟结果整体偏小。这是因为实验中冲击波的传播过程还受空气温度和湿度影响，并且钢板螺栓孔处会出现轻微拉伸变形，导致钢板中心挠度的实验结果偏大。另外，数值模型无法完全模拟实验现实状态。考虑到挠度的相对误差小于10%，在可接受的范围之内，因此可以认为数值模型是正确、可靠的。

3. 结　论

通过负压爆炸实验，对负压爆炸载荷作用下固支钢板动态响应的极限应变和失效条件进行了实验研究，同时利用AUTODYN非线性有限元程序进行数值模拟，得到以下主要结论。

(1) 实验结果显示，钢板整体出现塑性大变形，四周出现明显的拉伸变形。随着初始环境压力从101 kPa下降至60 kPa，半径在10～15 cm区间（边缘区）的挠度变化基本相同，而半径在0～10 cm区间（中心区）的挠度减小。

(2) 常压近爆作用下，钢板中心以拉伸断裂为主，裂纹延伸至边缘出现剪切断裂。根据双向应变假设，得到Q235钢板的动态极限应变为0.269。根据乳化炸药在不同负压环境下产生的爆炸冲击波的正规反射比冲量公式，结合基于刚塑性假设和能量准则提出的失效判据，可以有效判断负压爆炸载荷作用下钢板拉伸断裂情况。

(3) 不同环境压力下钢板挠度的数值模拟结果与实验结果的相对误差小于10%，验证了模型的准确性。模拟结果显示，随着初始环境压力从101 kPa下降至20 kPa，冲击波传播过程中空气介质的密度减小，钢板中心点最大挠度减小，中心点最大速度减小。

图 1 实验示意图

Figure 1. Schematic diagram of the experiment

下载: 全尺寸图片幻灯片

图 2 实验后钢板的形貌

Figure 2. Morphology of the steel plate after the experiment

下载: 全尺寸图片幻灯片

图 3 钢板变形轮廓

Figure 3. Deformation contours of steel plates

下载: 全尺寸图片幻灯片

图 4 工况4中钢板的破坏模式

Figure 4. Damage mode of steel plates in case 4

下载: 全尺寸图片幻灯片

图 5 几何模型

Figure 5. Geometrical model

下载: 全尺寸图片幻灯片

图 6 不同负压环境下钢板的变形情况

Figure 6. Deformation of steel plate under different negative pressure environments

下载: 全尺寸图片幻灯片

图 7 钢板中心点速度曲线

Figure 7. Velocity curves at center point of steel plate

下载: 全尺寸图片幻灯片

表 1 实验工况

Table 1. Experimental conditions

Case	p_e/kPa	W/g	TNT equivalence/g	Explosion distance/mm
1	60	150	130	150
2	80	150	130	150
3	101	150	130	150
4	101	200	170	50

下载: 导出CSV

表 2 开裂处钢板厚度的测量结果

Table 2. Measurement results of steel thickness at crack

Measuring point	Thickness/ mm	Down gauging rate	Bidirectional ultimate strain	Measuring point	Thickness/ mm	Down gauging rate	Bidirectional ultimate strain
1	0.60	0.40	0.291	5	0.64	0.36	0.250
2	0.62	0.38	0.270	6	0.54	0.46	0.361
3	0.52	0.48	0.387	7	0.66	0.34	0.231
4	0.68	0.32	0.213	8	0.70	0.30	0.195

下载: 导出CSV

表 3 失效判别式的相关参数

Table 3. Relevant parameters for failure discriminant

Case	i_r/(Pa∙s)	v₀/(m∙s⁻¹)	$\dot{\varepsilon }$ /s⁻¹	α	σ_d/MPa	ε_m	η
1	742.1	94.5	326.8	2.52	592.2	0.059	0.219
2	816.8	104.1	359.9	2.55	599.3	0.071	0.264
3	882.8	112.5	389.0	2.57	604.0	0.082	0.305
4	3204.1	408.2	1411.4	3.04	714.4	0.915	3.401

下载: 导出CSV

表 4 材料参数^[21]

Table 4. Material parameters^[21]

ρ_s/(g·cm⁻³)	A/MPa	B/MPa	n	C	m	T_m/K
7.85	293.8	230.2	0.578	0.0652	0.706	1795

下载: 导出CSV

表 5 不同环境压力对应的空气密度

Table 5. Air densities at different environment pressures

Pressure/kPa	Density/(kg·m⁻³)	Pressure/kPa	Density/(kg·m⁻³)
101	1.225	40	0.484
80	0.967	20	0.242
60	0.725

下载: 导出CSV

参考文献(21)

[1]	TEELING-SMITH R G, NURICK G N. The deformation and tearing of thin circular plates subjected to impulsive loads [J]. International Journal of Impact Engineering, 1991, 11(1): 77–91. doi: 10.1016/0734-743X(91)90032-B
[2]	GUPTA N K, NAGESH. Deformation and tearing of circular plates with varying support conditions under uniform impulsive loads [J]. International Journal of Impact Engineering, 2007, 34(1): 42–59. doi: 10.1016/j.ijimpeng.2006.05.002
[3]	JACOB N, NURICK G N, LANGDON G S. The effect of stand-off distance on the failure of fully clamped circular mild steel plates subjected to blast loads [J]. Engineering Structures, 2007, 29(10): 2723–2736. doi: 10.1016/j.engstruct.2007.01.021
[4]	陈长海, 朱锡, 侯海量, 等. 近距空爆载荷作用下固支方板的变形及破坏模式 [J]. 爆炸与冲击, 2012, 32(4): 368–373. doi: 10.11883/1001-1455(2012)04-0368-08 CHEN C H, ZHU X, HOU H L, et al. Deformation and failure modes of clamped square plates under close-range air blast loads [J]. Explosion and Shock Waves, 2012, 32(4): 368–373. doi: 10.11883/1001-1455(2012)04-0368-08
[5]	韩璐, 袁建飞, 张玉磊, 等. 固支矩形钢板近距爆炸的毁伤特性 [J]. 火炸药学报, 2021, 44(2): 225–232. doi: 10.14077/j.issn.1007-7812.202004023 HAN L, YUAN J F, ZHANG Y L, et al. Damage characteristics of fixed rectangular steel plate under close-in explosion [J]. Chinese Journal of Explosives & Propellants, 2021, 44(2): 225–232. doi: 10.14077/j.issn.1007-7812.202004023
[6]	WANG W, WANG Y P, YANG J C, et al. Investigation on air blast resistance of POZD-coated composite steel plates: experiment and numerical analysis [J]. Composites Part B: Engineering, 2022, 237: 109858. doi: 10.1016/j.compositesb.2022.109858
[7]	REMENNIKOV A, NGO T, MOHOTTI D, et al. Experimental investigation and simplified modeling of response of steel plates subjected to close-in blast loading from spherical liquid explosive charges [J]. International Journal of Impact Engineering, 2017, 101: 78–89. doi: 10.1016/j.ijimpeng.2016.11.013
[8]	秦业志, 王莹, 王志凯, 等. 小当量柱型装药水下近场爆炸固支单层方形钢板毁伤特性研究 [J]. 振动与冲击, 2021, 40(7): 29–36. doi: 10.13465/j.cnki.jvs.2021.07.004 QIN Y Z, WANG Y, WANG Z K, et al. Damage characteristics of fixed single-layer squaresteel plate under near-field underwater explosion of small equivalent column charge [J]. Journal of Vibration and Shock, 2021, 40(7): 29–36. doi: 10.13465/j.cnki.jvs.2021.07.004
[9]	CHEN S, QIN J, DENG S, et al. Experimental and numerical studies on fixed steel sheets subjected to underwater explosion [J]. Materials, 2022, 15(18): 6419. doi: 10.3390/MA15186419
[10]	张斐, 张春辉, 张磊, 等. 多次水下爆炸作用下钢板与焊接钢板冲击损伤特性 [J]. 振动与冲击, 2020, 39(7): 196–201. doi: 10.13465/j.cnki.jvs.2020.07.027 ZHANG F, ZHANG C H, ZHANG L, et al. Impact damage of steel plate and welding steel plate under multiple underwater explosions [J]. Journal of Vibration and Shock, 2020, 39(7): 196–201. doi: 10.13465/j.cnki.jvs.2020.07.027
[11]	赵铮, 陶钢, 杜长星. 爆轰产物JWL状态方程应用研究 [J]. 高压物理学报, 2009, 23(4): 277–282. doi: 10.11858/gywlxb.2009.04.007 ZHAO Z, TAO G, DU C X. Application research on JWL equation of state of detonation products [J]. Chinese Journal of High Pressure Physics, 2009, 23(4): 277–282. doi: 10.11858/gywlxb.2009.04.007
[12]	李孝臣, 汪泉, 谢守冬, 等. 负压条件下球形爆炸容器内乳化炸药冲击波参数研究 [J]. 火炸药学报, 2023, 46(3): 252–259. doi: 10.14077/j.issn.1007-7812.202207001 LI X C, WANG Q, XIE S D, et al. Study of shock wave parameters of emulsified explosives in spherical explosive containers under negative-pressure conditions [J]. Chinese Journal of Explosives & Propellants, 2023, 46(3): 252–259. doi: 10.14077/j.issn.1007-7812.202207001
[13]	冯晓伟, 卢永刚, 李永泽. 飞机目标在爆炸冲击波作用下的毁伤效应评估方法 [J]. 高压物理学报, 2019, 33(4): 045101. doi: 10.11858/gywlxb.20180687 FENG X W, LU Y G, LI Y Z. Damage assessment method of aircraft targets under blast wave [J]. Chinese Journal of High Pressure Physics, 2019, 33(4): 045101. doi: 10.11858/gywlxb.20180687
[14]	李瑞, 李孝臣, 汪泉, 等. 低温和低压环境下炸药爆炸冲击波的传播特性 [J]. 爆炸与冲击, 2023, 43(2): 022301. doi: 10.11883/bzycj-2022-0188 LI R, LI X C, WANG Q, et al. Propagation characteristics of blast wave in diminished ambient temperature and pressure environments [J]. Explosion and Shock Waves, 2023, 43(2): 022301. doi: 10.11883/bzycj-2022-0188
[15]	卢广照, 姜春兰, 毛亮, 等. 薄钢板在CL-20基含铝炸药内爆载荷作用下的变形响应和工程预测 [J]. 兵工学报, 2020, 41(8): 1509–1518. doi: 10.3969/j.issn.1000-1093.2020.08.005 LU G Z, JIANG C L, MAO L, et al. Deformation response and its engineering prediction of steel plate subjected to internal blast loading from CL-20-based aluminized explosive charges [J]. Acta Armamentarii, 2020, 41(8): 1509–1518. doi: 10.3969/j.issn.1000-1093.2020.08.005
[16]	牟金磊, 朱锡, 张振华, 等. 水下爆炸载荷作用下加筋板变形及开裂试验研究 [J]. 振动与冲击, 2008, 27(1): 57–60. doi: 10.13465/j.cnki.jvs.2008.01.009 MU J L, ZHU X, ZHANG Z H, et al. Experimental study on deformation and rupture of stiffened plates subjected to underwater shock [J]. Journal of Vibration and Shock, 2008, 27(1): 57–60. doi: 10.13465/j.cnki.jvs.2008.01.009
[17]	奥尔连科. 爆炸物理学(上册) [M]. 孙承纬, 译. 3版. 北京: 科学出版社, 2011: 463–472.
[18]	孙业斌. 爆炸作用与装药设计 [M]. 北京: 国防工业出版社, 1987: 221–228.
[19]	WEN H M. Deformation and tearing of clamped circular work-hardening plates under impulsive loading [J]. International Journal of Pressure Vessels and Piping, 1998, 75(1): 67–73. doi: 10.1016/S0308-0161(98)00023-4
[20]	胡志乐, 马亮亮, 吴昊, 等. 远距离近地面爆炸空气冲击波计算的网格尺寸优化与验证 [J]. 爆炸与冲击, 2022, 42(11): 114201. doi: 10.11883/bzycj-2021-0499 HU Z L, MA L L, WU H, et al. Optimization and verification of mesh size for air shock wave from large distance and near ground explosion [J]. Explosion and Shock Waves, 2022, 42(11): 114201. doi: 10.11883/bzycj-2021-0499
[21]	郭子涛, 高斌, 郭钊, 等. 基于J-C模型的Q235钢的动态本构关系 [J]. 爆炸与冲击, 2018, 38(4): 804–810. doi: 10.11883/bzycj-2016-0333 GUO Z T, GAO B, GUO Z, et al. Dynamic constitutive relation based on J-C model of Q235 steel [J]. Explosion and Shock Waves, 2018, 38(4): 804–810. doi: 10.11883/bzycj-2016-0333

施引文献

资源附件(0)

访问统计

点击查看大图

图(7) / 表(5)

计量

文章访问数: 217
HTML全文浏览量: 105
PDF下载量: 87

1. 负压爆炸实验
1.1 实验设计
1.2 固支钢板在不同负压环境下的变形破坏分析
1.3 动态极限应变
1.4 失效条件
2. 数值模拟
2.1 一维空中爆炸及结果映射
2.2 模型材料
2.3 数值模拟结果分析
3. 结　论

负压爆炸载荷作用下固支钢板变形研究

doi: 10.11858/gywlxb.20230685

作者简介: 杨 锐（1999－），男，硕士研究生，主要从事爆炸毁伤研究. E-mail：2240036819@qq.com

通讯作者: 汪 泉（1980－），男，博士，教授，博士生导师，主要从事爆炸理论与技术研究. E-mail：wqaust@163.com

计量

出版历程

Deformation of Fixed Support Steel Plate under Explosion Load in Negative Pressure Environment

1. 负压爆炸实验

1.1 实验设计

1.2 固支钢板在不同负压环境下的变形破坏分析

1.3 动态极限应变

1.4 失效条件

2. 数值模拟

2.1 一维空中爆炸及结果映射

2.2 模型材料

2.3 数值模拟结果分析

3. 结 论

计量

出版历程

目录

1. 负压爆炸实验

1.1 实验设计

1.2 固支钢板在不同负压环境下的变形破坏分析

1.3 动态极限应变

1.4 失效条件

2. 数值模拟

2.1 一维空中爆炸及结果映射

2.2 模型材料

2.3 数值模拟结果分析

3. 结 论

作者简介:
杨　锐（1999－），男，硕士研究生，主要从事爆炸毁伤研究. E-mail：2240036819@qq.com

通讯作者:
汪　泉（1980－），男，博士，教授，博士生导师，主要从事爆炸理论与技术研究. E-mail：wqaust@163.com

3. 结　论

3. 结　论