钢化玻璃冲击波毁伤效应测试结果分析

钟巍; 寿列枫; 何增; 李伟昌; 雷鸣; 刘俊; 田宙; 浦锡锋; 王仲琦

doi:10.11883/bzycj-2017-0070

钢化玻璃冲击波毁伤效应测试结果分析

doi: 10.11883/bzycj-2017-0070

1.
西北核技术研究所, 陕西西安 710024
2.
湘潭大学材料科学与工程学院, 湖南湘潭 411105
3.
北京理工大学爆炸科学与技术国家重点实验室, 北京 100081

基金项目:

国家自然科学基金项目 91330205

详细信息

作者简介:
钟巍(1986-), 男, 博士研究生, 助理研究员, zhongwei2016@pku.edu.cn

中图分类号: O383.2
计量
- 文章访问数: 5564
- HTML全文浏览量: 1411
- PDF下载量: 152
- 被引次数: 12
出版历程
- 收稿日期: 2017-03-02
- 修回日期: 2017-07-12
- 刊出日期: 2018-09-25

On experimental blast parameters for damage effect of monolithic tempered glass subjected to blast loading

1.
Northwest Institute of Nuclear Technology, Xi'an 710024, Shaanxi, China
2.
Faculty of Material Science and Engineering, Xiangtan University, Xiangtan 411105, Hunan, China
3.
State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China

摘要

摘要: 为研究爆炸冲击波对钢化玻璃的毁伤阈值，开展了钢化玻璃冲击波毁伤效应实验。对每一发实验进行了爆炸参数测试，获得了冲击波超压随时间变化历程的实验数据。通过对实验数据的处理与分析，得到了爆炸冲击波基本参数，包括正反射超压、冲击波到时、正压作用时间、正反射冲量、负反射超压、负压作用时间和负反射冲量等。将实验结果与CONWEP计算结果进行了比较，比较结果表明二者误差很小，证明冲击波测试结果准确可靠。对冲击波正负反射参数进行了比较，结果表明冲击波正反射参数明显高于冲击波负反射参数。
- 爆炸冲击波 /
- 钢化玻璃 /
- 毁伤阈值 /
- 爆炸参数 /
- CONWEP
Abstract: Damaging effects of blast waves on tempered glasses were experimentally investigated to explore the blast damage threshold of the tempered glasses. The blast shock wave parameters were tested in each experiment, and the experimental overpressure-time histories were obtained. By analyzing the experimental data, the blast shock wave parameters were obtained, such as normally reflected overpressure, arrival time of the shock wave, duration of the normally reflected overpressure, normally reflected impulse, negative reflected overpressure, duration of the negative reflected overpressure, and negative reflected impulse. Compared with the results computed by the CONWEP program, the relative errors were small, which verified that the tested results were accurate and reliable. Meanwhile, the normally reflected parameters are greater than the corresponding negative reflected ones.
- blast shock wave /
- tempered glass /
- damage threshold /
- blast parameter /
- CONWEP

HTML全文

现代鱼雷可通过多种发射方式攻击敌对目标, 如飞机空投或舰载发射等, 但不管采用何种发射方式, 都要经过空中、入水、水下航行至命中目标的全弹道过程。入水是一个短暂、变化激烈而又复杂的力学过程, 在入水的整个过程中, 鱼雷将经历撞水、浸水、带空泡航行、全浸湿后转入受控运动等4个阶段。入水过程中雷体会激起周围流体介质的运动, 反过来, 流体介质对结构又施加各种反作用力。入水冲击力有2个分量相当重要:一个是轴向力或阻力, 会引起减加速度, 可能使雷头变形, 破坏鱼雷内部的仪表设备; 另一个是法向力或升力, 会形成力矩, 产生横向角速度, 影响弹道, 并可能使雷体皱折或断裂。因此研究空投鱼雷入水冲击力对鱼雷结构的总体设计以及水下弹道设计等方面具有重要的工程意义^[1]。

本文中拟基于鱼雷入水动力学模型、动态非线性有限原理论和耦合算法^[2], 计算计算鱼雷入水时受到的作用力, 建立空投鱼雷入水冲击有限元模型, 对不同工况进行数值模拟, 以期为预测鱼雷入水冲击载荷提供参考。

1. 鱼雷入水受力分析

作用在鱼雷上的力可以分成2类:鱼雷重力和流体反作用力, 流体反作用力中包括浮力和流体动力^[3]。重力和浮力的大小和方向都是明确的, 容易确定的。流体动力不仅取决于鱼雷的外形、流体的特性, 也取决于鱼雷入水运动状态参数(如速度、加速度、姿态角、姿态角速度)等。可以用适合于惯性参考系的牛顿运动定律来描述鱼雷斜入水时纵平面内的受力及运动。鱼雷入水时运动和几何参数如图 1所示。

图 1 鱼雷入水运动简图

Figure 1. Schematic of torpedo entry into water

下载: 全尺寸图片幻灯片

浮力及其力矩为:

$\begin{array}{c} F_{\mathrm{b}}=\rho g \int_{-L_{2}}^{L_{1}} A_{\mathrm{b}}(l) \mathrm{d} l \end{array}$

(1)

$\begin{array}{c} M_{\mathrm{b}}=F_{\mathrm{b}} \int_{L}^{L_{1}} \sin (\phi) L(l) \mathrm{d} l \end{array}$

(2)

式中:A_b为沾湿部分有效截面积。

轴向黏性阻力、法向黏性阻力及其对重心力矩为:

$\begin{aligned} F_{\mathrm{d}, \mathrm{a}} & =\frac{1}{2} \rho C_{\mathrm{d}, \mathrm{a}} A_{\mathrm{e}} v_{\mathrm{a}}\left|v_{\mathrm{a}}\right| \end{aligned}$

(3)

$\begin{aligned} F_{\mathrm{d}, \mathrm{n}} & =\rho C_{\mathrm{d}, \mathrm{n}} \int_{L}^{L_{1}} r v_{\mathrm{n}}\left|v_{\mathrm{n}}\right| \mathrm{d} l \end{aligned}$

(4)

$\begin{aligned} M_{\mathrm{d}, \mathrm{n}} & =\rho C_{\mathrm{d}, \mathrm{n}} \int_{l}^{L_{1}} r l v_{\mathrm{n}}\left|v_{\mathrm{n}}\right| \mathrm{d} l \end{aligned}$

(5)

式中:v_a、v_n分别为轴向、法向速度, 轴向黏性阻力系数C_{d, a}=0.3, 法向粘性阻力系数C_{d, n}=1.1, A_e为轴向有效面积, r为轴向截面半径。

轴向、法向附加力及其对重心的力矩为:

$\begin{array}{c} F_{\mathrm{m}, \mathrm{a}}=\frac{\mathrm{d} m_{\mathrm{a}}}{\mathrm{d} v} \frac{\mathrm{d} v}{\mathrm{~d} t}+m_{\mathrm{a}} A_{\mathrm{a}}\left(\frac{L}{L_{1}+L_{2}}\right) \end{array}$

(6)

$\begin{array}{c} F_{\mathrm{m}, \mathrm{n}}=\int_{-L_{2}}^{L_{1}}\left(v_{\mathrm{n}} \frac{\mathrm{d} m}{\mathrm{~d} t}+m \frac{\mathrm{d} v_{\mathrm{n}}}{\mathrm{d} t}\right) \mathrm{d} l \end{array}$

(7)

$\begin{array}{c} M_{\mathrm{m}, \mathrm{n}}=\int_{-L_{2}}^{L_{1}} l \mathrm{~d} F_{\mathrm{m}, \mathrm{n}} \end{array}$

(8)

式中:m_a为轴向附加质量, m为单位长度的横截面在浸深h处的附加质量, 其计算方法可由文献[4]确定。根据以上受力分析, 建立鱼雷入水过程中动力学方程:

$\left\{\begin{array}{l} M X^{\prime \prime}=F_{\mathrm{m}, \mathrm{n}} \cos \phi-F_{\mathrm{m}, \mathrm{a}} \sin \phi-F_{\mathrm{m}, \mathrm{a}} \sin \phi+F_{\mathrm{m}, \mathrm{n}} \cos \phi \\ M Z^{\prime \prime}=F_{\mathrm{m}, \mathrm{n}} \sin \phi+F_{\mathrm{m}, \mathrm{a}} \cos \phi+F_{\mathrm{m}, \mathrm{a}} \sin \phi+F_{\mathrm{m}, \mathrm{n}} \sin \phi-M g \\ I_{\phi}^{\prime \prime}=M_{\mathrm{b}}+M_{\mathrm{m}, \mathrm{n}}+M_{\mathrm{d}, \mathrm{n}} \end{array}\right.$

(9)

式中:M为鱼雷总质量, I为对重心的转动惯量。将上式转化成矩阵形式:

$\left(\begin{array}{ccc} M+m_{11} & m_{12} & m_{13} \\ m_{21} & M+m_{22} & m_{23} \\ m_{31} & m_{32} & M+m_{33} \end{array}\right)\left(\begin{array}{c} X^{\prime \prime} \\ Z^{\prime \prime} \\ \phi^{\prime \prime} \end{array}\right)=\left(\begin{array}{c} F_{x} \\ F_{z} \\ F_{\phi} \end{array}\right)$

(10)

式中:m_ij(i, j=1, 2, 3)对应于加速度的质量项。式(10)是一个二阶常微分方程组, 将其降阶转化为一阶方程组, 采用变步长Runge-Kutta法求解, 对于各积分项采用变步长辛普森法进行计算。取MK46鱼雷为计算目标, 入水角为45°, 入水速度为32m/s, 攻角为零, 通过使用Matlab进行数值计算, 计算结果如表 1所示。

表 1 鱼雷入水时力学参数

Table 1. Mechnical parameters for torpedo entry into water

t/ms	F_{m, a}/kN	F_{m, n}/kN	M_{m, n}/(kN·m)
0.34	28.5	19.7	121.6
1.00	10.1	6.3	38.9
2.00	2.1	1.7	10.5
3.00	3.6	2.3	14.2
4.00	2.1	2.0	12.3
5.00	1.8	1.2	7.4
6.00	1.6	1.1	6.8

下载: 导出CSV

| 显示表格

由表 1可以看出, 鱼雷在入水的瞬间就会受到很大的轴向以及法向冲击载荷, 从而也产生了很大的法向弯矩, 这都会对鱼雷的结果以及姿态产生很大的影响。但持续时间很短, 随着不断入水, 轴向力和法向力迅速衰减, 趋于某一常值。

2. 数值模型

以MK46鱼雷为算例, 依据所研究问题的侧重点, 对鱼雷和流体进行了必要的简化。由于主要研究对象为鱼雷入水冲击时的流体动力情况, 不考虑载荷下鱼雷壳体会不会发生塑性变形以及结构破坏, 故将鱼雷作为刚体建模, 且不考虑分段。对于鱼雷入水, 由于其入水速度高, 粘性仅仅表现在附着在耦合面边界层很薄的一部分, 对数值计算的结果影响很小, 在本文中不予考虑^[5-7]。

鱼雷结构采用刚体假设, 泊松比μ=0.3, 质量m=235kg, 空气和水选用可压缩理想流体本构材料关系, 水压力由多项式状态方程描述。

对于压缩状态:

$p=a_{1} \mu+a_{2} \mu^{2}+a_{3} \mu^{3}+\left(b_{0}+b_{1} \mu\right) \rho_{0} e$

(11)

对于拉伸状态:

$p=a_{1} \mu+\left(b_{0}+b_{1} \mu\right) \rho_{0} e$

(12)

空气的压力由Gamma状态方程来描述:

$p=(\gamma-1) \rho_{1} e$

(13)

水密度为ρ₀=1 000kg/m³, 体积模量K=2.2GPa, 质量内能e=83.95kJ/kg; 空气密度ρ₁=1.293kg/m³, 质量热容比γ=1.4, 质量内能e=212.65kJ/kg

雷体结构采用拉格朗日四节点四边形壳建模, 共有4 800个节点, 4 798个单元, 有限元模型如图 2所示。水和空气选用欧拉六面体实体单元, 考虑到入水问题的具体情况, 在近水面处网格较密, 远离水面采用等差网格, 逐渐稀疏。鱼雷的整个外表面为流固耦合面, 采用一般耦合算法。鱼雷头部距液面初始高度0.2m。鱼雷入水(45°入水角)模型如图 3所示。

图 2 鱼雷有限元网格模型

Figure 2. Finite element model of torpedo

下载: 全尺寸图片幻灯片

图 3 鱼雷入水有限元模型

Figure 3. Finite element model for torpedo entry into water

下载: 全尺寸图片幻灯片

3. 鱼雷入水冲击载荷分析

3.1 鱼雷入水载荷分析

选取初始入水速度为32m/s, 入水角度为45°, 初始步长为1μm, 总计算时间为0.06s。将结果文件导入MSC.Dytran计算后可以得到该工况下鱼雷入水的分析结果。如图 4所示为鱼雷入水过程中耦合面受到的作用力时间历程曲线。

图 4 鱼雷入水耦合面作用力时间历程曲线

Figure 4. Coupling force-time curve for torpedo entry into water

下载: 全尺寸图片幻灯片

由图 4中可以看出, 鱼雷在入水时, 会遭受极大的冲击载荷, 随着鱼雷不断入水, 冲击载荷迅速的衰减下去, 但还是会不断受到冲击力的作用, 3ms后冲击过程基本结束, 作用力趋于稳定。作用力峰值达到35.397kN, 时间约为0.32ms, 由于作用时间很短, 不会影响鱼雷入水的姿态以及运动。计算得到的不同时刻作用力大小也与相同工况下Matlab计算出的结果(表 1)比较接近, 从而从进一步验证了数值模型的可靠性。

图 5所示为入水1.8ms时鱼雷外壳以及水所受到的压力云图。根据鱼雷在入水过程中受到的压力变化情况可得出, 入水初始时接触面的压力最大, 峰值为25.8MPa, 向边缘逐渐减小。对于平头鱼雷入水情况, 其撞水瞬间是一种碰撞现象, 撞水初期可以假设其为一平板撞击可压缩水面。基于Von Karman一元碰撞理论, 刚性平板撞击压缩水面的撞击压力峰值可由下式估算^[8]:

图 5 入水时鱼雷及水的压力云图

Figure 5. Pressure contour for torpedo entry into water

下载: 全尺寸图片幻灯片

$p_{\max }=\rho c v_{\perp} \sin \varphi$

(14)

式中:ρ为水的密度, 入水角度φ=45°, 水的声速c=1 480m/s, 垂直入水速度v_⊥=22.6m/s。计算得最大理论压力峰值为33.45MPa, 较计算结果高。这主要是由于没有考虑到空气被压入水中的空气垫效应而减小了压力。对于一般雷头壳体, 选用的材料为硬铝合金, 其屈服应力为330MPa, 故壳体所受压力值不可能使壳体发生塑性变形或破坏。

3.2 入水速度对鱼雷入水的影响

为了进一步分析入水过程, 以入水角度为45°的工况下分别给鱼雷加以从10~40m/s范围内变化的轴向入水速度进行计算。根据数值模拟结果, 将不同速度下的鱼雷头部所受到的压力峰值进行曲线拟合, 如图 6所示。

图 6 不同速度下鱼雷头部所受到的压力峰值

Figure 6. The pressure peak of torpedo head at different water-entry velocities

下载: 全尺寸图片幻灯片

所得拟合二次曲线函数:

$F=223 x^{2}+3523 x+12480 \quad 10 \mathrm{~m} / \mathrm{s} \leqslant x \leqslant 40 \mathrm{~m} / \mathrm{s}$

(15)

将不同速度下的鱼雷入水弹道曲线进行比对分析, 如图 7所示, 可以得出, 当冲击力作用在鱼雷头部时, 由于不通过重心, 使得鱼雷产生俯仰角速度, 对弹道产生严重影响。随着鱼雷入水速度的增大, 弹道偏离现象逐渐减小, 说明入水速度越高, 越有利于消除法向力引起的弯矩对鱼雷带来的影响。

图 7 不同入水速度鱼雷弹道曲线图

Figure 7. Torpedo trajectory curves at different water-entry velocities

下载: 全尺寸图片幻灯片

3.3 入水角度对鱼雷入水的影响

取入水速度为40m/s, 改变入水角度, 得到不同角度下的冲击力峰值变化曲线, 如表 2所示。

表 2 不同角度入水时鱼雷所受最大冲击力

Table 2. The maximum impact force at different water-entry angles

ϕ/(°)	F_{m, a}/kN	F_{m, n}/kN
30	23.9	21.8
45	30.9	19.2
60	36.5	11.7
75	38.5	6.5
90	68.0	0

下载: 导出CSV

| 显示表格

由表 2可知, 入水轴向冲击力随着入水角度的增大而增大, 尤其在垂直入水时轴向力急剧增大, 这说明垂直入水是鱼雷入水最不安全的情况, 会对鱼雷的纵向结构强度提出了很高的要求, 应该尽量避免。

而法向力随着入水角度的增大而减小, 这说明当鱼雷以小角度入水时应该注意避免法向力太大而引起的忽扑现象。同时根据计算结果分析不同角度入水弹道轨迹, 发现随着入水角度的变小, 鱼雷运动轨迹偏离初始角度的情况变的更为严重, 这也同样说明法向作用力及其所引起的弯矩在小角度入水情况下大于大角度入水时的现象。

3.4 不同头型对雷体入水的影响

分别选取了截头锥形壳体, 尖头壳体以及半圆头壳体这3种头部壳体形状的结构, 以入水速度30m/s对鱼雷垂直入水工况进行计算, 其有限元模型如图 8所示。

图 8 不同头型鱼雷有限元模型

Figure 8. Finite element models for torpedo with different headtypes

下载: 全尺寸图片幻灯片

3种不同头型鱼雷垂直入水工况的入水耦合面作用力时间历程曲线如图 9所示。

图 9 不同头型入水耦合面作用力时程曲线

Figure 9. Coupling force-time curve for water-entry torpedoes with different head shapes

下载: 全尺寸图片幻灯片

由图 9可以看出, 鱼雷在开始入水时都会受到很大的冲击载荷, 但截头锥形的头部形状受到的压力峰值最大, 压力作用的时间也是最长的, 其次是圆头雷, 最小的是尖头雷。这说明鱼雷入水时, 如果减小头部与水的接触面积, 则可以有效地减小鱼雷所受到的冲击载荷。

4. 结论

分析了鱼雷入水时受到的流体动力, 应用MATLAB编程求解了鱼雷入水过程的流体动力参数, 并通过使用MSC.Dytran仿真模拟了空投鱼雷入水过程。通过建立鱼雷入水三维有限元模型, 对鱼雷入水过程雷体所遭受的冲击力及其在不同入水速度和角度下冲击压力峰值、不同头部形状冲击压力峰值特点、鱼雷入水弹道做了分析与讨论。得出如下结论:

(1) 空投鱼雷入水时会受到很高的冲击载荷, 存在明显的冲击载荷峰值, 对雷头部结构将造成极大的影响, 随着鱼雷不断进入水中, 冲击载荷逐渐趋于常值。因此, 在设计鱼雷壳体的时候有必要强化其前端面附近的强度设计或在头部加装缓冲头帽, 提高其耐冲击的性能, 以确保其入水时的安全性;

(2) 随着鱼雷入水速度的增加, 作用在鱼雷壳体上的冲击载荷也不断升高, 这对鱼雷雷体的强度提出了很高的要求, 所以应尽可能减小空投鱼雷在入水时的速度, 如在空投时雷身附有降落伞;

(3) 随着鱼雷入水角度的增加, 作用在鱼雷上的法向作用力不断减小, 但轴向力不断增大, 从而由于法向弯矩引起的忽扑现象逐渐消失。垂直入水时轴向作用力达到最高值, 这是鱼雷入水最危险的工况, 应该尽量避免;

(4) 对于垂直入水工况, 不同头型鱼雷的入水载荷各不相同, 减小鱼雷与水的接触面积会明显减小入水载荷峰值的大小。

图 1 PCB压力传感器的安装

Figure 1. Installation of PCB pressure sensor

下载: 全尺寸图片幻灯片

图 2 理想爆炸超压时间历程曲线及其三角形简化模型

Figure 2. The ideal explosion shock wave overpressure-time curve and the triangle simplification model

下载: 全尺寸图片幻灯片

图 3 爆炸冲击波传播示意图

Figure 3. A diagram of the blast shock propagation

下载: 全尺寸图片幻灯片

图 4 第1轮实验正反射超压时间历程曲线与相应的CONWEP程序计算结果的比较

Figure 4. Comparison of the normally reflected overpressure-time curves obtained in the first round experimentswith the corresponding ones computed by CONWEP

下载: 全尺寸图片幻灯片

图 5 第2轮实验正反射超压时间历程曲线与相应的CONWEP程序计算结果的比较

Figure 5. Comparison of the normally reflected overpressure-time curves obtained in the second round experimentswith the corresponding ones computed by CONWEP

下载: 全尺寸图片幻灯片

图 6 第3轮实验正反射超压时间历程曲线与相应的CONWEP程序计算结果的比较

Figure 6. Comparison of the normally reflected overpressure-time curves obtained in the third round experimentswith the corresponding ones computed by CONWEP

下载: 全尺寸图片幻灯片

图 7 正反射超压随比例爆心距的变化

Figure 7. Normally reflected overpressure varying with scaled distance of blast

下载: 全尺寸图片幻灯片

图 8 比例冲击波到时随比例爆心距的变化

Figure 8. Scaled shock wave arrival time varying with scaled distance of blast

下载: 全尺寸图片幻灯片

图 9 比例正压作用时间随比例爆心距的变化

Figure 9. Scaled duration of normally reflected overpressure varying with scaled distance of blast

下载: 全尺寸图片幻灯片

图 10 比例正反射冲量随比例爆心距的变化

Figure 10. Scaled normally reflected impulse varyingwith scaled distance of blast

下载: 全尺寸图片幻灯片

图 11 正负反射超压随比例爆心距的变化

Figure 11. Normally and negative reflected overpressuresvarying with scaled distance of blast

下载: 全尺寸图片幻灯片

图 12 比例正、负反射超压作用时间随比例爆心距的变化

Figure 12. Scaled duration of normally and negative reflected overpressures varying with scaled distance of blast

下载: 全尺寸图片幻灯片

图 13 比例正、负反射冲量随比例爆心距的变化

Figure 13. Scaled normally and negative reflected impulsesvarying with scaled distance of blast

下载: 全尺寸图片幻灯片

表 1 实验中使用的TNT装药和玻璃的基本信息

Table 1. Basic information of the TNT charge and glass in each experiment

实验序号	W/kg	σ_c/MPa	l/mm	w/mm	h/mm	备注
1-1-1	5	69	1 600	1 300	10	不带夹层
1-2-2	5	69	1 600	1 300	10	不带夹层
1-3-3	5	69	1 600	1 300	10	不带夹层
1-4-4	5	69	1 600	1 300	10	不带夹层
1-5-5	5	69	1 600	1 300	10	不带夹层
1-6-6	5	69	1 600	1 300	10	不带夹层
1-7-7	5	69	1 600	1 300	10	不带夹层
1-8-8	5	69	1 600	1 300	10	不带夹层
1-9-9	5	69	1 600	1 300	6	不带夹层
1-10-10	5	69	1 600	1 300	6	不带夹层
1-11-11	5	69	1 600	1 300	6	不带夹层
1-12-12	5	69	1 600	1 300	6	不带夹层
1-13-13	5	69	1 600	1 300	6	不带夹层
1-14-14	5	69	1 600	1 300	10	不带夹层
2-1-15	10	69	1 600	1 300	15	不带夹层
2-2-16	10	69	1 600	1 300	15	不带夹层
2-3-17	10	69	1 600	1 300	15	不带夹层
2-4-18	5	69	1 600	1 300	6	不带夹层
2-5-19	5	69	1 600	1 300	6	不带夹层
2-6-20	5	69	1 600	1 300	6	不带夹层
2-7-21	5	69	1 600	1 300	10	不带夹层
2-8-22	5	69	1 600	1 300	10	不带夹层
2-9-23	10	69	1 600	1 300	10	不带夹层
2-10-24	10	69	1 600	1 300	10	不带夹层
2-11-25	10	69	1 600	1 300	10	不带夹层
2-12-26	10	69	1 600	1 300	10	不带夹层
2-13-27	10	69	1 600	1 300	10	不带夹层
3-1-28	10	69	1 600	1 300	6+1.52PVB+6	带PVB夹层
3-2-29	10	69	1 600	1 300	6+1.52PVB+6	带PVB夹层
3-3-30	10	69	1 600	1 300	6+1.52PVB+6	带PVB夹层
3-4-31	10	69	1 600	1 300	6+1.52PVB+6	带PVB夹层
3-5-32	10	69	1 600	1 300	6+1.52PVB+6	带PVB夹层

下载: 导出CSV

表 2 位置O、A、B、C处的冲击波参数计算结果

Table 2. Calculation results of shock wave parameters for locations O, A, B and C

OP/m	测点	R/m	p_s/kPa		I_s/(kPa·ms)
OP/m	测点	R/m	5 kg	10 kg	5 kg	10 kg
4.000	O	4.000	415.30		335.10
	A	4.075	394.60		327.90
	B或C	4.052	400.80		330.10
10.000	O	10.000	52.74	80.88	120.20	194.70
	A	10.030	52.46	80.40	119.80	194.10
	B或C	10.021	52.55	80.56	120.00	194.30
15.000	O	15.000	28.02	39.57	77.99	125.60
	A	15.020	27.96	39.49	77.88	125.40
	B或C	15.014	27.98	39.51	77.91	125.50
20.000	O	20.000	18.82	25.86	57.59	92.55
	A	20.015	18.81	25.83	57.55	92.47
	B或C	20.011	18.81	25.84	57.56	92.49
25.000	O	25.000	13.92	19.03	45.57	73.17
	A	25.012	13.91	19.01	45.55	73.13
	B或C	25.008	13.91	19.02	45.56	73.14

下载: 导出CSV

表 3 第1轮实验正反射冲击波参数测试结果和CONWEP程序计算结果

Table 3. Normally reflected blast parameters obtained in the first round experimentsand the corresponding ones computed by CONWEP

实验序号	W/kg	R/m	(R·W^-1/3)/(m·kg^-1/3)	方法	p_s⁺/kPa	t_a/ms	t_d⁺/ms	I_s⁺/(kPa·ms)
1-14-14	5	4	2.339	实验	429.00	4.22	2.87	615.61
1-14-14	5	4	2.339	CONWEP	415.30	4.42	3.78	335.10
1-4-4	5	8	4.678	实验	115.28	13.56	4.31	248.43
1-4-4	5	8	4.678	CONWEP	115.30	12.72	6.30	231.20
1-9-9	5	8	4.678	实验	119.51	13.31	4.39	262.32
1-9-9	5	8	4.678	CONWEP	115.30	12.72	6.30	231.20
1-6-6	5	10	5.848	实验	70.51	18.00	5.19	182.97
1-6-6	5	10	5.848	CONWEP	74.70	17.84	6.86	180.70
1-7-7	5	11	6.433	实验	59.02	20.07	5.68	167.62
1-7-7	5	11	6.433	CONWEP	63.08	20.47	7.09	162.80
1-8-8	5	11.5	6.725	实验	43.83	21.13	3.91	85.69
1-8-8	5	11.5	6.725	CONWEP	58.47	21.80	7.19	155.20
1-8-8修正	5	11.5	6.725	实验	57.12	21.91	5.69	162.51
1-8-8修正	5	11.5	6.725	CONWEP	58.47	21.80	7.19	155.20
1-1-1	5	12	7.018	实验	57.47	23.50	5.27	151.43
1-1-1	5	12	7.018	CONWEP	54.45	23.14	7.29	148.20
1-5-5	5	12	7.018	实验	55.22	23.75	5.71	157.64
1-5-5	5	12	7.018	CONWEP	54.45	23.14	7.29	148.20
1-3-3	5	15	8.772	实验	36.65	31.29	6.11	111.96
1-3-3	5	15	8.772	CONWEP	38.29	31.25	7.84	116.60
1-10-10	5	15	8.772	实验	39.27	31.70	6.00	117.86
1-10-10	5	15	8.772	CONWEP	38.29	31.25	7.85	116.60
1-13-13	5	16	9.357	实验	36.41	33.90	6.52	118.66
1-13-13	5	16	9.357	CONWEP	34.77	33.98	8.01	108.90
1-2-2	5	17	9.942	实验	27.62	36.43	3.87	53.44
1-2-2	5	17	9.942	CONWEP	23.50	38.58	7.17	68.33
1-12-12	5	17	9.942	实验	25.31	36.92	5.29	66.94
1-12-12	5	17	9.942	CONWEP	23.50	38.58	7.17	68.33
1-11-11	5	20	11.696	实验	19.10	46.58	5.80	55.39
1-11-11	5	20	11.696	CONWEP	18.82	46.90	7.53	57.59

下载: 导出CSV

表 4 第2轮实验正反射冲击波参数测试结果和CONWEP程序计算结果

Table 4. Normally reflected blast parameters obtained in the second round experimentsand the corresponding ones computed by CONWEP

实验序号	W/kg	R/m	(R·W^-1/3)/(m·kg^-1/3)	方法	p_s⁺/kPa	t_a/ms	t_d⁺/ms	I_s⁺/(kPa·ms)
2-1-15	10	6	2.785	实验	262.02	7.36	4.39	575.13
2-1-15	10	6	2.785	CONWEP	260.40	7.56	5.53	345.00
2-2-16	10	10	4.642	实验
2-2-16	10	10	4.642	CONWEP	80.88	17.34	6.99	194.70
2-3-17	10	14	6.498	实验	52.25	26.90	6.41	167.46
2-3-17	10	14	6.498	CONWEP	61.99	26.17	8.96	202.90
2-9-23	10	15	6.962	实验	48.32	29.36	7.81	188.69
2-9-23	10	15	6.962	CONWEP	55.17	28.83	9.17	188.30
2-13-27	10	16	7.427	实验
2-13-27	10	16	7.427	CONWEP	49.64	31.52	9.37	175.60
2-12-26	10	17	7.891	实验
2-12-26	10	17	7.891	CONWEP	45.07	34.22	9.55	164.50
2-8-22	5	15	8.772	实验	34.07	31.81	6.40	109.02
2-8-22	5	15	8.772	CONWEP	38.29	31.25	7.85	116.60
2-11-25	10	20	9.283	实验	29.97	41.68	6.91	103.55
2-11-25	10	20	9.283	CONWEP	35.18	42.38	10.07	138.30
2-4-18	5	17	9.942	实验	23.68	37.37	5.94	70.33
2-4-18	5	17	9.942	CONWEP	23.50	38.58	7.17	68.33
2-10-24	10	25	11.604	实验	24.02	55.34	8.84	106.17
2-10-24	10	25	11.604	CONWEP	25.53	56.15	10.82	109.20
2-5-19	5	20	11.696	实验	17.00	47.29	6.52	55.42
2-5-19	5	20	11.696	CONWEP	18.82	46.90	7.53	57.59
2-7-21	5	20	11.696	实验	20.80	44.29	5.87	61.05
2-7-21	5	20	11.696	CONWEP	18.82	46.90	7.53	57.59
2-6-20	5	25	14.620	实验	14.68	58.61	6.18	45.36
2-6-20	5	25	14.620	CONWEP	13.92	60.86	8.03	45.57

下载: 导出CSV

表 5 第3轮实验正反射冲击波参数测试结果和CONWEP程序计算结果

Table 5. Normally reflected blast parameters obtained in the third round experimentsand the corresponding ones computed by CONWEP

实验序号	W/kg	R/m	(R·W^-1/3)/(m·kg^-1/3)	方法	p_s⁺/kPa	t_a/ms	t_d⁺/ms	I_s⁺/(kPa·ms)
3-1-28	10	5	2.32	实验	394.40	4.50	2.20	433.84
3-1-28	10	5	2.32	CONWEP	424.40	5.48	4.72	426.10
3-2-29	10	8	3.71	实验	139.40	12.10	4.50	313.65
3-2-29	10	8	3.71	CONWEP	129.70	12.24	6.41	249.20
3-3-30	10	10	4.60	实验	106.50	15.10	5.70	303.53
3-3-30	10	10	4.60	CONWEP	117.20	15.82	7.91	293.90
3-4-31	10	15	6.96	实验	54.60	28.70	6.30	171.99
3-4-31	10	15	6.96	CONWEP	55.17	28.83	9.17	188.30
3-5-32	10	18	8.35	实验	41.20	37.20	7.30	135.78
3-5-32	10	18	8.35	CONWEP	41.24	36.93	9.73	154.80

下载: 导出CSV

表 6 第1轮实验负反射冲击波参数测试结果

Table 6. Negative reflected blast parameters obtained in the first round experiments

实验序号	W/kg	R/m	(R·W^-1/3)/(m·kg^-1/3)	p_s^-/kPa	t_d^-/ms	I^-/(kPa·ms)
1-1-1	5	12	7.018	16.10	3.12	25.07
1-2-2	5	17	9.942	14.29	2.43	17.36
1-3-3	5	15	8.772	7.47	2.44	9.12
1-4-4	5	8	4.678	20.31	3.60	36.52
1-5-5	5	12	7.018	13.62	3.57	24.34
1-6-6	5	10	5.848	18.70	4.89	45.74
1-7-7	5	11	6.433	14.94	3.13	23.36
1-8-8	5	11.5	6.725	14.87	2.43	18.07
1-8-8修正	5	11.5	6.725	14.28	2.43	36.13
1-9-9	5	8	4.678	27.75	3.76	52.22
1-10-10	5	15	8.772	11.96	3.62	21.67
1-11-11	5	20	11.696
1-12-12	5	17	9.942
1-13-13	5	16	9.357	9.37	3.37	15.78
1-14-14	5	4	2.339

下载: 导出CSV

表 7 第2轮实验负反射冲击波参数测试结果

Table 7. Negative reflected blast parameters obtained in the second round experiments

实验序号	W/kg	R/m	(R·W^-1/3)/(m·kg^-1/3)	p_s^-/kPa	t_d^-/ms	I^-/(kPa·ms)
2-1-15	10	6	2.785	50.14	4.41	110.56
2-2-16	10	10	4.642
2-3-17	10	14	6.498	9.59	4.32	20.71
2-4-18	5	17	9.942	5.55	4.07	11.29
2-5-19	5	20	11.696	2.96	4.55	6.73
2-6-20	5	25	14.620	5.66	5.49	15.54
2-7-21	5	20	11.696	5.92	3.22	9.53
2-8-22	5	15	8.772	7.66	4.81	18.42
2-9-23	10	15	6.962	7.77	4.60	17.87
2-10-24	10	25	11.604	5.79	6.43	18.61
2-11-25	10	20	9.283	7.22	6.84	24.69
2-12-26	10	17	7.891
2-13-27	10	16	7.427

下载: 导出CSV

表 8 第3轮实验负反射冲击波参数测试结果

Table 8. Negative reflected blast parameters obtained in the third round experiments

实验序号	W/kg	R/m	(R·W^-1/3)/(m·kg^-1/3)	p_s^-/kPa	t_d^-/ms	I^-/(kPa·ms)
3-1-28	10	5	2.32
3-2-29	10	8	3.71	27.40	4.10	56.17
3-3-30	10	10	4.60	15.40	4.40	33.88
3-4-31	10	15	6.96	12.10	5.10	30.86
3-5-32	10	18	8.35	8.60	6.80	29.24

下载: 导出CSV

参考文献(14)

[1]	李峰.玻璃幕墙结构在爆炸荷载作用下的动力响应分析[D].西安: 长安大学, 2009: 1-7. http://cdmd.cnki.com.cn/Article/CDMD-11941-2009210690.htm
[2]	BRUCK T. The economic consequences of new global terrorism[J]. Economic Bulletin, 2002, 39(10):327-332. doi: 10.1007/s10160-002-0168-8
[3]	CORLEY W G, MLAKAR P F, SOZEN M A, et al. The Oklahoma city bombing:summary and recommendations for multihazard mitigation[J]. Journal of Performance of Constructed Facilities, 1998, 12(3):100-112. doi: 10.1061/(ASCE)0887-3828(1998)12:3(100)
[4]	MLAKAR P F, CORLEY W G, SOZEN M A, et al. The Oklahoma city bombing:analysis of blast damage to the murrah building[J]. Journal of Performance of Constructed Facilities, 1998, 12(3):113-119. doi: 10.1061/(ASCE)0887-3828(1998)12:3(113)
[5]	NORVILLE H S, HARVILL N, CONRATH E J, et al. Glass-related injuries in Oklahoma city bombing[J]. Journal of Performance of Constructed Facilities, 1999, 13(2):50-56. doi: 10.1061/(ASCE)0887-3828(1999)13:2(50)
[6]	ZHANG Xihong, HAO Hong, WANG Zhongqi. Experimental investigation of monolithic tempered glass fragment characteristics subjected to blast loads[J]. Engineering Structures, 2014, 75:259-275. doi: 10.1016/j.engstruct.2014.06.014
[7]	ZHANG Xihong, HAO Hong, MA Guowei. Parametric study of laminated glass window response to blast loads[J]. Engineering Structures, 2013, 56:1707-1717. doi: 10.1016/j.engstruct.2013.08.007
[8]	IDRISS J S, LOWAK M J. Empirical evaluation of glazing systems in response to blast loads[C]//Proceedings of the 2014 Structure Congress. Boston, USA, 2014: 258-269. doi: 10.1061/9780784413357.024
[9]	中国建筑科学研究院.JGJ102-2003玻璃幕墙工程技术规范[S].北京: 中国建筑工业出版社, 2003: 60-68.
[10]	常规武器防护设计原理(TM5-855-1)[M].方秦, 吴平安, 张育林, 译.南京: 工程兵工程学院, 1997: 1-20.
[11]	姜鹏飞, 唐德高, 伍俊.装配式组合防爆墙抗高速破片侵彻性能研究[J].振动与冲击, 2008, 27(7):131-135. doi: 10.3969/j.issn.1000-3835.2008.07.030 JIANG Pengfei, TANG Degao, WU Jun. Analysis on penetration characteristics of assembled blast-resistance wall impacted by fragments with high velocities[J]. Journal of Vibration and Shock, 2008, 27(7):131-135. doi: 10.3969/j.issn.1000-3835.2008.07.030
[12]	丁娜娜, 严东晋, 徐源.悬臂式防爆墙消波作用数值分析及其应用[J].地下空间与工程学报, 2010, 6(3):645-649. http://d.old.wanfangdata.com.cn/Periodical/dxkj201003037 DING Nana, YAN Dongjin, XU Yuan. Numerical analysis and application of shielding effect of cantilever blast wall[J]. Chinese Journal of Underground Space and Engineering, 2010, 6(3):645-649. http://d.old.wanfangdata.com.cn/Periodical/dxkj201003037
[13]	夏志成, 王静, 施翠英.内爆炸作用下抗爆单元防护隔墙厚度计算[J].弹箭与制导学报, 2009, 29(1):144-147. doi: 10.3969/j.issn.1673-9728.2009.01.041 XIA Zhicheng, WANG Jing, SHI Cuiying. Calculation of thickness of blast wall of antiknock cell under internal explosion[J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2009, 29(1):144-147. doi: 10.3969/j.issn.1673-9728.2009.01.041
[14]	倪晋峰.单层球面网壳基于CONWEP外爆响应分析及试验设计[D].哈尔滨: 哈尔滨工业大学, 2012: 13-19. http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=D238759

施引文献

期刊类型引用(7)

1.	李昂. 爆炸载荷下玻璃破坏的研究. 信息记录材料. 2022(02): 12-14 . 百度学术
2.	谢东升，孙滔，史卓鹏，智生龙，史淋升，李海涛. 输电线路爆破除冰动态特性. 爆炸与冲击. 2021(06): 151-159 . 本站查看
3.	张文岚，杜礼明，章芝霖，李梓豪. 爆炸荷载作用下高速列车车窗玻璃的动态响应. 大连交通大学学报. 2021(05): 7-12 . 百度学术
4.	钟巍，田宙，寿列枫. 基于线性代数的大规模快速量纲分析算法及其在爆炸与冲击工程研究中的应用. 计算数学. 2020(02): 170-195 . 百度学术
5.	王澍霏，钟巍，王智环，姚成宝，贾雷明，田宙. 爆炸荷载作用下单层钢化玻璃破坏状态的理论预测方法. 现代应用物理. 2020(03): 87-93 . 百度学术
6.	钟巍，寿列枫，王仲琦，田宙，刘俊，何增，李伟昌. PVB夹层钢化玻璃冲击波毁伤效应实验研究. 北京理工大学学报. 2019(06): 565-570 . 百度学术
7.	钟巍，寿列枫，田宙，何增，李伟昌，刘俊. 10 mm厚无夹层钢化玻璃冲击波毁伤效应实验研究. 科学技术与工程. 2019(34): 48-55 . 百度学术

其他类型引用(5)

资源附件(0)

访问统计

图(13) / 表(8)

计量

文章访问数: 5564
HTML全文浏览量: 1411
PDF下载量: 152
被引次数: 12

1. 鱼雷入水受力分析
2. 数值模型
3. 鱼雷入水冲击载荷分析
3.1 鱼雷入水载荷分析
3.2 入水速度对鱼雷入水的影响
3.3 入水角度对鱼雷入水的影响
3.4 不同头型对雷体入水的影响
4. 结论

钢化玻璃冲击波毁伤效应测试结果分析

doi: 10.11883/bzycj-2017-0070

作者简介: 钟巍(1986-), 男, 博士研究生, 助理研究员, zhongwei2016@pku.edu.cn

计量

出版历程

On experimental blast parameters for damage effect of monolithic tempered glass subjected to blast loading

1. 鱼雷入水受力分析

2. 数值模型

3. 鱼雷入水冲击载荷分析

3.1 鱼雷入水载荷分析

3.2 入水速度对鱼雷入水的影响

3.3 入水角度对鱼雷入水的影响

3.4 不同头型对雷体入水的影响

4. 结论

期刊类型引用(7)

其他类型引用(5)

计量

出版历程

目录

1. 鱼雷入水受力分析

2. 数值模型

3. 鱼雷入水冲击载荷分析

3.1 鱼雷入水载荷分析

3.2 入水速度对鱼雷入水的影响

3.3 入水角度对鱼雷入水的影响

3.4 不同头型对雷体入水的影响

4. 结论

作者简介:
钟巍(1986-), 男, 博士研究生, 助理研究员, zhongwei2016@pku.edu.cn